Hi All:
Not sure why my previous question never got posted. Here I am seeking some
help on my code. I am using the following code to run MCMC simulation on
the following data using the model below:
# Data
matrix<-NULL> csvs<-paste("MVN", 1:2,".csv",sep="")
> for (i in 1:length(csvs)){
+ matrix[[i]]<- read.csv(file=csvs[i],header=TRUE)
+ print(matrix[[i]])
+ }
Y1 Y2
1 11 6
2 8 5
3 25 13
4 1 13
5 8 22
Y1 Y2
1 9 1
2 7 9
3 25 13
4 1 18
5 9 12
#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))
expit[1]<-exp(gamma[1])/(1+exp(gamma[1]))
expit[2]<-exp(gamma[2])/(1+exp(gamma[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))
}
# 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" for this format
Thanks for your help.
Anamika
[[alternative HTML version deleted]]
