grmPhenoEval function

function to take a data.frame coming from framePhenoRec and GRM and analyze them with individuals as a random effect with a GRM covariance matrix

grmPhenoEval(phenoDF, grm)

Arguments

phenoDF	A data.frame of phenotypic observations. See framePhenoRec for details
grm	A genomic relationship matrix

Value

Named real vector of the BLUPs of all individuals in phenoDF (names are the individual ids), with appropriate weights by error variance of the observation

Details

Given all the phenotypic records calculate the GEBV for each individual using all its records

Examples

phenoDF <- framePhenoRec(records, bsp)
#> Error in framePhenoRec(records, bsp): object 'bsp' not found
grm <- makeGRM(records, bsp, SP)
#> Loading required package: sommer
#> Loading required package: Matrix
#> Loading required package: MASS
#> 
#> Attaching package: ‘MASS’
#> The following object is masked from ‘package:dplyr’:
#> 
#>     select
#> Loading required package: lattice
#> Loading required package: crayon
#> Error in makeGRM(records, bsp, SP): object 'records' not found
grmBLUPs <- grmPhenoEval(phenoDF, grm)
#> Error in factor(phenoDF$id, levels = rownames(grm)): object 'phenoDF' not found