[BioC] GRENITS won't run with example data

Tyler Funnell tyler.funnell at gmail.com
Mon Mar 26 04:44:38 CEST 2012


Hello,

I've been having trouble getting the GRENITS package to run so I tried running just their example code on their included sample data.
Here is the outputP

> library(GRENITS)
Loading required package: Rcpp
Loading required package: RcppArmadillo
Loading required package: ggplot2
> data(Athaliana_ODE)
> output.folder <- paste(tempdir(), "/Example_LinearNet", sep="")
> LinearNet(output.folder, Athaliana_ODE)
Created folder /var/folders/ty/9l9tb3g54614cx23_clq3tnr0000gn/T//RtmpTPnsMW/Example_LinearNet
Started MCMC chain 1  ============= 
Estimated runtime = 0.0531906 min
25% completed
50% completed
75% completed
MCMC chain 1 finished!
Started MCMC chain 2  ============= 
Estimated runtime = 0.0521771 min
25% completed
50% completed
75% completed
MCMC chain 2 finished!
> analyse.output(output.folder)
Error in UseMethod("melt", data) : 
  no applicable method for 'melt' applied to an object of class "c('matrix', 'double', 'numeric')"
> 

I noticed on the project page that it says it imports the reshape library. When I load the reshape library and try the example it works.
However, when I set the output.folder to a directory in my desktop and run the analysis I get a segfault.

> output.folder <- paste("~/Desktop", "/Example_LinearNet", sep="")
> LinearNet(output.folder, Athaliana_ODE)
Created folder ~/Desktop/Example_LinearNet
Started MCMC chain 1  ============= 

 *** caught segfault ***
address 0x10, cause 'memory not mapped'

Traceback:
 1: .Call("callAR1", timeSeries, resultsFolder.i, ParamVec, fixMe,     PACKAGE = "GRENITS")
 2: LinearNet(output.folder, Athaliana_ODE)


Cheers,
Tyler



More information about the Bioconductor mailing list