Question: SPIA plotP giving error
gravatar for docpriya21
2.4 years ago by
docpriya210 wrote:


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 


However, when I use


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)

[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       




spia plotp threshold • 334 views
ADD COMMENTlink written 2.4 years ago by docpriya210

I have the same problem/question. Looking at the code:

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:

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


ADD REPLYlink written 12 weeks ago by maxininfa0
Please log in to add an answer.


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