search for: numchains

...2[i] Y[i] ~ dnorm(c,d) It seems doesn't work in my model. The problem is between two ######. The error message is > modelCheck("BayesBioMarker.BUGS") model is syntactically correct > modelData(paste("BUGS_data.txt",sep="")) data loaded > modelCompile(numChains=1) multiple definitions of node bm[1] =================================== model { for(iter in 1:numSubj){ bmd[iter] ~ dnorm(u[trt[iter]+1],sdbmd.precision); ############### bm[iter] <- inprod(biomarker[iter,1:4], b...

...quot;mu","tau") inits<-function() {list(beta=beta0, rho_o=rho_o_0, rho_d=rho_d_0, rho_w=rho_w_0)} parameters<-c("beta", "rho_o", "rho_d", "rho_w") probit.sim<-BRugsFit(data,inits,parameters,modelFile="spatial.openbugs.txt",numChains=1,nIter=2000) # my model model { for (i in 1:676) { y[i] ~ dbern(p[i]) wwy[i]<- inprod(ww[i, 1:676] , y[]) woy[i]<- inprod(wo[i, 1:676] , y[]) wdy[i]<- inprod(wd[i, 1:676] , y[]) probit(p[i])<- rho_o * woy[i] + rho_d * wdy[i] + rho_w * wwy[i] + beta[1] + beta[2] * cap2[i] + be...

...rtation/R/Dyes/Model.txt' > model.data.name = 'C:/Users/ethan/Desktop/Dissertation/R/Dyes/data.txt' > n.chains = 1 > write.model(model, model.file.name) > modelCheck(model.file.name) model is syntactically correct > modelData(model.data.name) data loaded > modelCompile(numChains=n.chains) model compiled > model.init.name = bugsInits(inits) > modelInits(rep(model.init.name,n.chains)) Initializing chain 1: initial values loaded but chain contain uninitialized variables > modelGenInits() initial values generated, model initialized > modelUpdate(25000) 25000 update...

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).

Hi All, I am trying to figure out how to use R-Coda with the output from OpenBugs. I have installed and loaded the packages BRugs and R2WinBUGS. I have successfully run a simple Bayes model in WinBUGS using R2WinBUGS' "bugs" and have used "read.bugs" to build the coda object. I can successfully switch to OpenBugs and run the same model and get the basic summary

....Monitor [01C59510H] .sampleMonitor SamplesMonitors.Monitor [01C0A4E0H] UpdaterActions.UpdateMonitors [0000011CH] .chain INTEGER 0 .monitor UpdaterActions.Monitor [01C59510H] BugsCmds.Action.Step [00000BCAH] .a BugsCmds.Action [011D40C0H] .chain INTEGER 0 .i INTEGER 8 .j INTEGER 1000 .numChains INTEGER 1 BugsCmds.Action.Do [00000C58H] .a BugsCmds.Action [011D40C0H] .elapsedTime LONGINT 62EC21F162EB0384H .profileList ARRAY 2024 OF CHAR "" ... .res INTEGER 1635734111 .s ARRAY 120 OF CHAR 3X ... Services.Exec [00000136H] .a Services.Action [011D40C0H] .t POINTER [6...

...gam<-rep(1,Lstra), u<-matrix(rep(0, 30), nrow=2) ) } bi.bmi.model.parameters=c( 'a0', 'a1', 'a2', 'a3', 'a12') bi.bmi.model.bugs=BRugsFit(modelFile='bi.bmi.model2.txt', data=bi.bmi.model.data, inits=bi.bmi.model.init, numChains=1, para=bi.bmi.model.parameters, nBurnin=20, nIter=40) When I run this I get this message. model is syntactically correct data loaded array index is greater than array upper bound for a1 [1] "C:\\DOCUME~1\\Owner\\LOCALS~1\\Temp\\RtmpNvSdyb/inits1.txt" Initializing chain 1:...

...ins longer than 100") if ((100 * round(nBurnin/100) != nBurnin) | (100 * round(nIter/100) != nIter)) warning("chain lengths not multiples of 100; truncation may occur.") modelCheck(modelFile) modelData(bugsData(data)) modelCompile(numChains = 1) modelInits(bugsInits(inits)) burninIncSize <- 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) {...