Question: Combined p-values for DE genes from RNA-seq meta-analysis
0
gravatar for ks
2.9 years ago by
ks30
ks30 wrote:

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 fisher meta analysis • 1.2k views
ADD COMMENTlink modified 2.9 years ago by chris86250 • written 2.9 years ago by ks30

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 REPLYlink written 2.9 years ago by chris86250

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

ADD REPLYlink written 2.9 years ago by ks30
0
gravatar for chris86
2.9 years ago by
chris86250
United Kingdom, London
chris86250 wrote:

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 COMMENTlink modified 2.9 years ago • written 2.9 years ago by chris86250
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1998 users visited in the last hour