[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