Entering edit mode
8.2 years ago
star
▴
350
I have RNAseq data of patient and control samples with different age and gender. I would like to know whether these co-variate affect gene expression levels. If they do not have any effect remove them from model and create simpler model for differential expression analysis using edgeR and DESeq2. Would you please let me know if there is a function in these packages to find the answer? Actually, I would like to see details of models in edgeR and DESeq2 and then find which co-varate affect gene expression and which one do not have affect?
Coincidentally a colleague told me today about http://bioconductor.org/packages/release/bioc/html/pcaExplorer.html, which would allow you to analyze (the effect of) covariates on your data. I haven't tested it, but perhaps it's something worth trying for your question.