Displaying 3 results from an estimated 3 matches for "ln_model_tri".
Did you mean:
ln_model_trial
2010 Aug 16
1
Specify decimal places for parameters in BUGS output
...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...
2010 Aug 10
0
[BUGS] [R-BUGS] error while plotting
...<- 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)
> > plot(bugs.out...
2013 Jun 24
0
Running MCMC using R2WinBUGS
...;-structure(.Data = c(.001,0,0,.001),.Dim = c(2, 2))
}
# R2WinBUGS code to run MCMC
library("R2WinBUGS")
bugs.output <- list()
for(i in 1:2){
Y<-(matrix[[i]])
bugs.output[[i]] <- bugs(
data=list(Y=Y, Nf=5), # change for no of sites
model.file="M-LN_model_trial.txt",
parameters.to.save = c("p","rho","sigma2"),
n.chains=1, n.iter=12000, n.burnin=5000,
bugs.directory="H://AChaudhuri/winbugs14/WinBUGS14",
working.directory=NULL)
}
Error in FUN(X[[1L]], ...) :
.C(..): 'type' must be "real&qu...