Displaying 1 result from an estimated 1 matches for "mcmcout".
2010 Jun 23
1
A question about R2Winbugs
...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)
{
# cluster...