Goseq pwf nullp error
Hi all,
I got some trouble while use goseq to do enrichment analysis recently, the trouble is:
> pwf=nullp(de.data,bias.data=gene.length)
Error in if (min(fv) < lower_bound) fv = fv - min(fv) + lower_bound :
missing value where TRUE/FALSE needed
while using nullp function to calculate probability weight, I get the errors, could somebody help me?
ps: no NA existed in dataframe de.data, length(de.data) == length(gene.length)
Thanks a lot!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Blog: The Complete Guide to Metatranscriptome Sequencing
written by Novogene •
-
Goseq - Pwf: Error In If (Hi <= Low) { : Missing Value Where True/False Needed
written by ryme •
-
GOseq gene2cat format errors
written by et395 •
-
Too many enriched GO terms using Goseq
written by tianshenbio
-
GOseq analysis with non-native organism
written by Bin •
-
GOseq Probability Weight Function Interpretation and Use
written by matt.a.bennett25890 •
-
How to do gene length bias corrected KEGG pathway analysis for non-native species in goseq?
written by CandiceChuDVM
-
GOSeq (very few enriched genes with p-value less than 0.05)
written by yfyian •
-
Goseq: Pwf: Error In If (Hi <= Low) { : Missing Value Where True/False Needed
written by liupfskygre
What is your organism? How are you producing
gene.length?Hi Parham, sorry for my delay, the organism is sheep, and the gene.length is as below:
Did I make any errors?
Thanks a lot!
If I can help you I have to see the codes you have before
pwf=nullp(de.data,bias.data=gene.length).Hi Parham,
The codes before
pwf-nullp()is as below:Its not very clear to me what you are doing. One thing I can say is, do the
bias.dataandde.datahave the samelength?Yes, no
NAexisted in dataframede.data, andlength(de.data) == length(gene.length)Could you please provide an example of :
genes_lisths.fa.len