MVP.HE.Vg.Ve {rMVP} | R Documentation |
Build date: Feb 2, 2017 Last update: Feb 2, 2019
MVP.HE.Vg.Ve(y, X, K)
y |
phenotype |
X |
genotype |
K |
kinship matrix |
vg, ve, and delta
Translated from C++(GEMMA, Xiang Zhou) to R by: Haohao Zhang
phePath <- system.file("extdata", "07_other", "mvp.phe", package = "rMVP") phenotype <- read.table(phePath, header=TRUE) print(dim(phenotype)) genoPath <- system.file("extdata", "06_mvp-impute", "mvp.imp.geno.desc", package = "rMVP") genotype <- attach.big.matrix(genoPath) print(dim(genotype)) K <- MVP.K.VanRaden(genotype) vc <- MVP.HE.Vg.Ve(y=phenotype[,2], X=matrix(1, nrow(phenotype)), K=K) print(vc)