Combined p-values for DE genes from RNA-seq meta-analysis
1
0
Entering edit mode
8.5 years ago
ks ▴ 30

Hello, I have raw p-values from edgeR for two comparisons of interest: control vs genotype1 and control vs genotype2. I'd like to see if there is an aggregate effect of genotype on DE.

The Internet has a number of answers on how to do a Fisher's combined test, and some resources on how to do this in R (such as the metaRNASeq package). But edgeR (as well as DESeq2) give two-sided p-values, while the R functions for combined tests requires a one-sided p-value. Does anyone know how to proceed? I'm also not sure if this is a valid way to get an aggregate result, since these two tests have the same control samples and are not independent. Does anyone have thoughts on whether it would be better to run the analysis in edgeR again, treating genotype1 and genotype2 as one group?

RNA-Seq meta analysis fisher • 2.9k views
ADD COMMENT
0
Entering edit mode

You haven't explained which studies are separate? Are control vs genotype1 and control vs genotype2 two completely separate experiments with new controls and new genotype samples?

ADD REPLY
0
Entering edit mode

No, they are from the same experiment -- the samples sequenced were control, genotype1, and genotype2.

ADD REPLY
0
Entering edit mode
8.5 years ago
chris86 ▴ 400

You can use limma for RNA-seq analysis now it has been updated. You can see in Aarons answer on this page he describes how to do contrasts for a A+B versus C comparison which is what you seem to want to do. I think this is what you are looking for - at a guess.

con.2 <- makeContrasts(A - (B+C)/2, levels=design)

https://support.bioconductor.org/p/69764/

ADD COMMENT

Login before adding your answer.

Traffic: 1801 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6