Question: PWF from GOseq to adjust q-value in Cuffdiff
gravatar for trakhtenberg
5.8 years ago by
United States
trakhtenberg160 wrote:


In accounting for gene length bias using GOseq, the PWF is used to adjust the p-value for determining the over-represented GO category.

Has anyone made a script which would use this PWF for adjusting the q or p value in the Cuffdiff output files (I use cufflinks) of all DEGs?

Then it would be straightforward to just select all DEGs below the updated significance value threshold for all sorts of downstream analysis (including GO enrichment, KEGG pathways, or just for having a more accurate list of the DEGs to select candidates for screening)?

Thank you.

rna-seq goseq gene length bias • 1.6k views
ADD COMMENTlink modified 5.8 years ago by Devon Ryan98k • written 5.8 years ago by trakhtenberg160
gravatar for Devon Ryan
5.8 years ago by
Devon Ryan98k
Freiburg, Germany
Devon Ryan98k wrote:

GOseq doesn't use the PWF to adjust the p-values. Rather, it uses the PWF to construct the background distribution, which is then used in the statistical test. To use a similar method with cuffdiff you'd have to completely rewrite it, this isn't a simple correction. I won't go into the fact that what GOseq is trying to compensate for isn't an issue at the level of analysis for which one uses cuffdiff.

What I think you're concerned about is GC/length bias in library construction and variance in that between samples/experimental groups. If that's something you're concerned about, then absolutely do not use cuffdiff. There's already a package to help compensate for these sorts of things, called CQN. You can't use that with cuffdiff, but you can with edgeR/DESeq2/limma.

ADD COMMENTlink written 5.8 years ago by Devon Ryan98k
Please log in to add an answer.


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