search for: parameterstosav

...AL DECOMPOSITION") B <- spline.des(allKnots, x, derivs = rep(drv, length(x)), outer.ok = TRUE)$design Z <- B %*% LZ attr(Z, "range.x") <- range.x attr(Z, "intKnots") <- intKnots return(Z) } BRugsMCMC<- function (data, inits, parametersToSave, nBurnin, nIter, nThin, modelFile) { if ((nBurnin < 100) | (nIter < 100)) stop("currently only working for chains longer than 100") if ((100 * round(nBurnin/100) != nBurnin) | (100 * round(nIter/100) != nIter)) warning("chain length...

...eater than array upper bound for a1 [1] "C:\\DOCUME~1\\Owner\\LOCALS~1\\Temp\\RtmpNvSdyb/inits1.txt" Initializing chain 1: model must be compiled before initial values loaded model must be initialized before updating model must be initialized before DIC an be monitored Error in samplesSet(parametersToSave) : model must be initialized before monitors used I checked the main effects alone, and it works fine, so I don't really understand why it's saying, "array index is greater than array upper bound for a1." Anyone who could help me with this would be greatly appreciated. Thanks...

...ot defined in model or in data set [1] "C:\\DOCUME~1\\maomao\\LOCALS~1\\Temp\\RtmpqJk9R3/inits1.txt" Initializing chain 1: model must be compiled before initial values loaded model must be initialized before updating model must be initialized before DIC an be monitored Error in samplesSet(parametersToSave) : model must be initialized before monitors used ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ I did define variable "ww" in my data and model (they are listed below). I am not sure if this is due to some errors in my code (please see below)...

...tial values loaded Initializing chain 3: model must be compiled before initial values loaded model must be compiled before generating initial values Sampling has been started... model must be initialized before updating model must be initialized before DIC an be monitored Error in BRugs::samplesSet(parametersToSave) : model must be initialized before monitors used Thanks for your help! John R commands datalist<-list("r", "N","n","Lev", "Light", "Soil1", "Soil2", "Soil3", "Consp", "Plot", "DConsp&...