Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Incorrect phenotype labelling of Beta Matrix for multiple phenotypes in K Models #5

Open
cfljam opened this issue Feb 13, 2025 · 0 comments

Comments

@cfljam
Copy link
Contributor

cfljam commented Feb 13, 2025

  • When running mulitple phenotypes, the beta matrix heading is always phenotype 1
  • doesnt seem to be case w Q models
  • reprex with test data
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)

Check the outputs

> 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
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant