Hi All:

I had a basic question to ask. I am running R2WinBUGS so that I could
automate the running of my model using 1000 simulated datasets. Below is the
code I am using. The only problem I am having is the bugs output that comes
out shows my parameters as nos with 1 decimal place after. I would want to
have the parameters with 5 places after decimal. How would I specify that in
my code for R2WinBUGS.

trial.data <- read.table("trial_data.txt", header=T)
bugs.output <- list()
for(i in 1:5){
       nausea <- as.integer(trial.data[i,])
       bugs.output[[i]] <- bugs(
       data=list(nausea=nausea, N=63),
       inits=list(
               list(alpha0=0,tau=1),
               list(alpha0=0.9, tau=1),
               list(alpha0=0.8,tau=1)
               ),
       model.file="B-LN_model_trial.txt",
       parameters.to.save = c("sigma", "tau","alpha0","p"),
n.chains=3, n.iter=12000, n.burnin=5000,
bugs.directory="E:/AChaudhuri/winbugs14/WinBUGS14")}
print(bugs.output)

The output comes out as
Inference for Bugs model at "B-LN_model_trial.txt", fit using WinBUGS,
 2 chains, each with 12000 iterations (first 5000 discarded), n.thin = 14
 n.sims = 1000 iterations saved
          mean   sd  2.5%   25%   50%   75% 97.5% Rhat n.eff
sigma      1.2  0.2   1.0   1.1   1.2   1.3   1.6    1   630
tau        0.7  0.2   0.4   0.6   0.7   0.8   1.0    1   630
alpha0    -2.0  0.2  -2.3  -2.1  -2.0  -1.9  -1.7    1   330
p[1]       0.2  0.0   0.1   0.1   0.2   0.2   0.3    1  1000
p[2]       0.1  0.0   0.0   0.1   0.1   0.1   0.2    1  1000
p[3]       0.0  0.0   0.0   0.0   0.0   0.1   0.1    1  1000
etc

> Thanks
>
Anamika

	[[alternative HTML version deleted]]

