search for: codapkg

...nits.beta0,taueps=1.0E-3)} parameters <- list("sigma","beta","beta0","y.star") fitm <- bugs(data,inits,parameters,model.file="model.bug", n.chains=3, n.iter=n.iter, n.burnin=n.burnin, n.thin = n.thin, debug=FALSE,DIC=FALSE,digit=5,codaPkg=FALSE, bugs.directory="C:/Program Files/WinBUGS14/") but I always get the following error: Error in FUN(X[[1L]], ...) : .C(..): 'type' must be "real" for this format I tried the web, but failed. Could anyone give me a clue? Best! -- View this message in con...

...Bugs and run the same model and get the basic summary back. However, I cannot build the coda object. From what I have read "read.openbugs" should do the trick. However, I have discovered the current version of OpenBUGS (v 3.03) doesn't write any output files in spite of setting codaPkg=TRUE in the "bugs" call. "Str(openbugs.object)" reveals a large collection of elements, one of which is an array that contains the MCMC chains. Is there a function to build the coda object? Thanks in advance for any suggestions, Bill Halteman

...at, digits = 0) parm <- c("lbda") bugs(dat, inits=list(NULL), parm, "d2.bug", n.chains = 1, n.iter = 5000, n.burnin = floor(n.iter/2), n.thin = max(1, floor(n.chains * (n.iter - n.burnin)/1000)), bin = (n.iter - n.burnin) / n.thin, debug = TRUE, DIC = TRUE, digits = 5, codaPkg = FALSE, bugs.directory = "C:/WinBUGS14/", working.directory = NULL, clearWD = FALSE) The objective of the program is to compare means of two independent samples that results in Beherens-Fisher posterior and in the model.file of WinBUGS "d2.bug" there is the following c...

..."Error: cannot allocate vector of size 216.0 Mb". This is an extract of my R syntax: > memory.limit(4000) NULL > simu.BYM<-bugs(data, inits, parameters, model.file="Modelo.txt", n.chains=3, n.iter=200000, n.burnin=20000, n.thin=180, debug=FALSE, DIC=FALSE, digits=5, codaPkg=FALSE, bugs.directory="c:/Archivos de programa/WinBUGS14", working.directory=NULL) Error: cannot allocate vector of size 216.0 Mb I have tried: -Change size of virtual memory inside Windows XP -change the R memory settings -run from DOS (R CMD BATCH...) -reviewed the R Windows FAQ and R-...

...get the basic >> summary back. However, I cannot build the coda object. From what I have >> read "read.openbugs" should do the trick. However, I have discovered the >> current version of OpenBUGS (v 3.03) doesn't write any output files in spite >> of setting codaPkg=TRUE in the "bugs" call. "Str(openbugs.object)" reveals >> a large collection of elements, one of which is an array that contains the >> MCMC chains. Is there a function to build the coda object? >> >> Thanks in advance for any suggestions, >> &gt...

...and diff2) are returned. The code used to obtain these estimates is as follows: sim <- bugs(data,init,model.file="drift.bugs",parameters=c("theta1","theta2","diff1","diff2"),n.chains=n.chains,n.iter=n.reps,n.burnin=1000,bin=1,n.thin=1,digits=3, codaPkg=T) sim.list <- read.bugs(sim) However, when I read the coda files using read.bugs the order of the columns in the parameter matrix are as follows: [1] "deviance" "diff1[1]" "diff1[10]" "diff1[100]" "diff1[11]" "diff1[12]&quot...

...eta2=0, beta3=0)} parameters <- c("beta0", "beta1", "beta2", "beta3") #call winbugs using bugs() function model <- bugs(data, inits, parameters, model.file="C:/Jinmodel/test_codes/model_poisson.txt", n.chains=1, n.iter=30, n.burnin=0,n.thin=1, codaPkg=TRUE, DIC=TRUE, debug=TRUE, bugs.directory="C:/Users/Yiyi/Desktop/winbugs14/WinBUGS14") ########################################### I can get pretty good traceplots & parameter estimate outputs in WinBUGS. However, when I switched back to the R Console (by pressing ESC key), a warn...

...(){list(a0=rbeta(1,1,1), qx=matrix(rgamma(K^2,alpha,1),nrow=K), sigma0.r=rbeta(K,1,1))} model.sim=bugs(data,inits,model.file="model.txt",parameters=c("mu","sigma"), n.chains=3,n.iter=3500,n.burnin=500,n.thin=1,bugs.directory="C:/Users/t-wec/Desktop/WinBUGS14",codaPkg=T,debug=T) mcmcout=read.bugs(model.sim) summary(mcmcout) ------------------------------------------------------------------------------------------------------------------------ model { # cluster parameters mu and tau tau1 <- (1-a*a)*taue taue <- 1/sigmae for (k in 1:K)...

...<- c("gamma0", "delta0", "beta1", "alpha1","alpha2", "index") model <- bugs(data, inits, parameters, model.file = "C:/Documents and Settings/del125/My Documents/PhD data/ETBF - BUGS/CodeR2BUGS/MinusMonthSeaWiFS.bug", codaPkg = FALSE, bugs.directory = "C:/Program Files/WinBUGS14/", debug =TRUE) ##And the Error Message from WinBUGS model is syntactically correct data(C:/Program Files/R/R-2.6.0/data.txt) data loaded compile(3) made use of undefined node nSeaW inits(1,C:/Program Files/R/R-2.6.0/ini...

...quot;,"m") parameters<-c("b","tau") a<-bugs(data=data, inits=inits, parameters, model.file="c:/try/sur.txt", n.chains = 3, n.iter = 1000, bugs.directory = "c:/Program Files/WinBUGS14/", working.directory = "c:/try", clearWD = FALSE,codaPkg = FALSE,debug=T) model { for (i in 1:n) { Y[i,1:2] ~ dmnorm(mu[i,],P[1:2,1:2]) # means in separate time series mu[i,1] <- inprod(X1[i,],b[,1]) mu[i,2] <- inprod(X2[i,],b[,2]) P[1:2,1:2]<-tau*I[1:2,1:2] } # priors on regr...

...inits=function(){list(alpha=rgamma(1,1.5,1),theta=mvrnorm(1,theta0,S2),T=rwish(3,S3),K=6)} dp_norm.sim=bugs(data,inits,model.file="dp_norm.txt",parameters=c("alpha","K"), n.chains=1,n.iter=500,n.burnin=50,n.thin=1,bugs.directory="C:/Program Files/WinBUGS14/",codaPkg=T,debug=T) -- View this message in context: http://n4.nabble.com/a-small-question-about-R-with-Winbugs-tp1788668p1788668.html Sent from the R help mailing list archive at Nabble.com.

...im = bugs(data,inits, model.file="model.bug", parameters=c("T","best","d","lor","mu","or","p"), n.chains=1,n.iter=20000,n.burnin=5000,n.thin=1, bugs.directory="c:/Program Files/WinBUGS14/", codaPkg=FALSE, debug=FALSE) Thank you in advance for any help you may be able to provide! Sincerely, Isabella Isabella R. Ghement, Ph.D. Ghement Statistical Consulting Company 301-7031 Blundell Road, Richmond, B.C., Canada, V6Y 1J5 Tel: 604-767-1250 Fax: 604-270-3922...