search for: samplesstat

...9;d like to run a bunch of models via BRugs, saving the state after each one, like so: for (m in models) { fit <- BRugsFit(paste(m$modelname, ".bug"), ...) save(fit, paste(m$modelname, ".RData)) } After I load the saved fitted objects in a new R process, BRugs functions like samplesStats() no longer work. Is it possible to somehow load the saved objects and use samplesStats and other BRugs functions?

Hello R users, I have used the R package BRugs in some applications but one thing that annoy me is the functions to retrieve samples ("samplesSample") or summary statistics ("samplesStats"). They take a really long time to return the values if the size of the chains in question is moderate, for example 10.000 samples. Does this happens with every one or there is a better way to obtain such values? Best -- Thiago Guerrera. [[alternative HTML version deleted]]

...ated, model initialized > modelUpdate(25000) 25000 updates took 1 s > samplesSet(c("sigma2.with","sigma2.btw")) monitor set for variable 'sigma2.with' monitor set for variable 'sigma2.btw' > #dicSet() > modelUpdate(100000) 100000 updates took 5 s > samplesStats("*") mean sd MC_error val2.5pc median val97.5pc start sample sigma2.btw 681.9 1161 10.89 0.7016 253.8 4232 25001 100000 sigma2.with 4266.0 1246 4.92 2480.0000 4057.0 7262 25001 100000 > ############################################ #############...

Dear all, I just try to run the example in R brugs packages, but not only once. Loop is added in this example. After several times (7, 11, or other random number), there is an error message "r for windows gui front-end has stopped working". This happened in two laptops with windows 7 and vista. I have tried different versions of R (2.14.1 and 2.13) and Brugs packages (0.7.1 and 0.7-5).

...round(nBurnin/(nThin * 100)) for (i in 1:100) { modelUpdate(burninIncSize, thin = nThin) } samplesSet(parametersToSave) iterIncSize <- round(nIter/(nThin * 100)) for (i in 1:100) { modelUpdate(iterIncSize, thin = nThin) } sims <- samplesStats("*") return(list(Stats = sims)) } require(BRugs) numKnots <- 25 range.x <- c(0,202) ox <- order(x) x <- x[ox] delta <- delta[ox] n <- length(x) xTil <- as.numeric(names(table(x))) cnts <- as.num...

Hello, I want to know if there are some functions or packages to solve differential and integral equation using R. Thanks. Shao chunxuan. [[alternative HTML version deleted]]