HCBravoLab / metagenomeSeq

Statistical analysis for sparse high-throughput sequencing
66 stars 20 forks source link

Model subject ID in differential analysis using fitFeatureModel #80

Open K-Pimentel opened 4 years ago

K-Pimentel commented 4 years ago

I would like to include subject ID in the design matrix of a fitFeatureModel as I am analysing data from a crossover study and it is thus useful to control for subject ID. As far as I can see the fitFeatureModel does not (yet?) offer the same flexibility as the model.matrix as accepted by fitZig()

I receive an error message when I attempt the analysis: Can't analyze currently

Are more complex study designs supported by the fitFeatureModel() function or is this planned for the future? Thank you.