Hi All:
I am running R2WinBUGS so that I could automate the running of my
Multivariate Logit Normal model on 300 simulated datasets. The problem I
have is I am not able to read in the data correctly and then code it in the
automation part.
*Heres my model:*
model
{
for (j in 1 : Nf)
{
p1[j, 1:2 ] ~ dmnorm(gamma[1:2], T[1:2 ,1:2])
# T is the precision matrix or inverse sigma for MVN
for (i in 1:2)
{
logit(p[j,i])<-p1[j,i]
Y[j,i] ~ dbin(p[j,i],n)
}
}
# Hyper-priors:
gamma[1:2] ~ dmnorm(mn[1:2],prec[1:2 ,1:2])
mn<-c(-1.59,-2.44)
prec<-structure(.Data = c(.001,0,0,.001),.Dim = c(2, 2))
T[1:2 ,1:2] ~ dwish(R[1:2 ,1:2], 2)
sigma2[1:2, 1:2] <-inverse(T[,])
rho <- sigma2[1,2]/sqrt(sigma2[1,1]*sigma2[2,2])
R<-structure(.Data = c(.001,0,0,.001),.Dim = c(2, 2))
}
As of now my *Y is in 2 seperate datasets*,generated using the 2
different p's.I am using the first 5 simulated datsets for now.My question
is should I put the two Y's together in one dataset and then read that. If
so, how do I read it?
*Heres my automation code:*
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, Nf=20),
model.file="model.txt", n.iter=12000, n.burnin=5000,
bugs.directory="E:/AChaudhuri/winbugs14/WinBUGS14")}
Appreciate any help.
Thanks
Anamika
[[alternative HTML version deleted]]