library(mpQTL)
## form a phenotype matrix with additional column
mppheno2 <- cbind(mppheno,phenotype2=sample(mppheno[,1]))
## Fit K model
result_mix2 <- map.QTL(phenotypes = mppheno2,
genotypes = mphapdose,
ploidy = 4,
map = mpmap,
K = T,
binsize = 1)
> head(result_mix2$phenotype1$beta[[1]])
phenotype1
161.7023575
69 2.3290261
35 0.3849862
50 1.7214127
9 0.9017646
42 0.9288540
> head(result_mix2$phenotype2$beta[[1]])
phenotype1
161.7023575
69 -1.1577997
35 0.5484196
50 -0.7968072
9 -0.4091906
42 -1.1837361