SPIA plotP giving error
0
1
Entering edit mode
docpriya21 ▴ 10
@docpriya21-11876
Last seen 8.0 years ago

Hello

I'm using SPIA on my toptable generated by limma. I have R version 3.3.1, and updated biocLite to the latest version. It makes the results file using 

res=spia(de=DE,all=ALL,organism="hsa",nB=2000,plots=FALSE,beta=NULL,combine="fisher")

However, when I use

plotP(res,threshold=0.05)

I get the following error message, and it doesn't plot. 

Error in plotP(res, threshold = 0.05) : 
  The threshold value should be 0.173107798207408 or higher!!!

It plots when I use 0.1732 for threshold but no points fall to the right of the oblique line. I'm not sure what else to do. 

Here is the output of sessionInfo():=

R version 3.3.1 (2016-06-21)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.12.1 (Sierra)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] SPIA_2.26.0          KEGGgraph_1.32.0     BiocInstaller_1.24.0

loaded via a namespace (and not attached):
[1] parallel_3.3.1      tools_3.3.1         BiocGenerics_0.20.0
[4] stats4_3.3.1        XML_3.98-1.5        graph_1.52.0       

 

Thanks

Priya

SPIA plotP threshold • 1.3k views
ADD COMMENT
0
Entering edit mode

I have the same problem/question. Looking at the code: https://rdrr.io/bioc/SPIA/src/R/plotP.R

if(threshold<x[1,"pgfdr"]){ msg&lt;-paste("the="" threshold="" value="" should="" be",x[1,"pgfdr"],"or="" higher!!!");="" stop(msg);="" }<="" p="">

Then at the explanation: https://www.rdocumentation.org/packages/SPIA/versions/2.24.0/topics/spia

"pGFdr and pGFWER are the False Discovery Rate and respectively Bonferroni adjusted global p-values."

So...

ADD REPLY

Login before adding your answer.

Traffic: 768 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6