Question: GOseq Probability Weight Function Interpretation and Use
gravatar for matt.a.bennett25890
15 days ago by
matt.a.bennett258900 wrote:


Looking to get some feedback on a probability weight function (PWF) curve I've generated using GOseq. As I understand it, as long as the curves fit the data points reasonably then the downstream analysis will be ok however mine look a bit strange compared to examples I've seen elsewhere. This seems to be mainly due to a few extreme outliers and I'm wondering whether these will be affecting the results unduly and whether there is a way to account for these within the nullp() function. Also whether I've made some mistake in the gene length determination... Here is my process

A) Get length data for Hg38 from Ensembl and filter for only my expressed genes (filtered by all genes with count >1 in treatment of interest - too low?):

txdb <- makeTxDbFromEnsembl("Homo sapiens", release = 89)
lengthfilt <- lengthData[names(lengthData) %in% names(gene_list)]

B) Use nullp() to generate a PWF for the expressed genes based on their respective lengths

pwf = nullp(clustVsAll, "hg38", "ensGene", = lengthfilt)

pwf curve

Any help/feedback much appreciated!

rna-seq goseq gene ontology • 60 views
ADD COMMENTlink modified 15 days ago • written 15 days ago by matt.a.bennett258900
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: 991 users visited in the last hour