Skip to content

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

@cfljam

Description

@cfljam
  • 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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions