Hello,
I'm trying to run the goseq on a not supported organism but unfortunately I already stuck for a while with the nullp function...
I have my DE genes in a vector:
head(degenes)
AT1G21660 AT1G11310 AT1G06220 AT1G01010 AT1G01020 AT1G01030
1 0 0 0 0 0
and also my gene_length:
head(gene_length)
AT1G21660 AT1G11310 AT1G06220 AT1G01010 AT1G01020 AT1G01030
3366 3447 4909 2268 2809 2065
if I run now the nullp function:
pwf<-nullp(degenes,bias.data=gene_length)
Error in if (hi <= low) { : missing value where TRUE/FALSE needed
Its not too good explained in the goseq manual so I hope anybody can help me.
Cheers,
Julian
I'm trying to run the goseq on a not supported organism but unfortunately I already stuck for a while with the nullp function...
I have my DE genes in a vector:
head(degenes)
AT1G21660 AT1G11310 AT1G06220 AT1G01010 AT1G01020 AT1G01030
1 0 0 0 0 0
and also my gene_length:
head(gene_length)
AT1G21660 AT1G11310 AT1G06220 AT1G01010 AT1G01020 AT1G01030
3366 3447 4909 2268 2809 2065
if I run now the nullp function:
pwf<-nullp(degenes,bias.data=gene_length)
Error in if (hi <= low) { : missing value where TRUE/FALSE needed
Its not too good explained in the goseq manual so I hope anybody can help me.
Cheers,
Julian