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]]