在R中运行R2WinBUGS时出错:对于此格式，“type”必须为"real“

Question

我在R中运行R2WinBUGS时遇到错误。这里是可重复性的代码。

#Set working directory 
setwd("H://AChaudhuri/Testing/CSVS")

matrix=NULL 
csvs <- paste("MVN", 1:2, ".csv", sep="") 
for(i in 1:length(csvs)){ 
matrix[[i]] <- read.csv(file=csvs[i], header=T) 
print(matrix[[i]]) 
} 

现在我已经读入了2个模拟数据集，它们看起来像

  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 

我的下一步是在这些数据集上运行多变量logit正态模型，并为这些模拟数据集自动执行此过程。下面是模型语句：

model 
 { 
for (j in 1 : Nf) 

  { 
  p1[j, 1:2 ] ~ dmnorm(gamma[1:2 ], T[1:2 ,1:2 ]) 

  for (i in 1:2) 
  { 
 logit(p[j,i])<-p1[j,i] 

Y[j,i] ~ dbin(p[j,i],n) 
}} 

我正在尝试使用以下代码在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 
inits=NULL, 
   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

> traceback()
6: .C("str_signif", x = x, n = n, mode = as.character(mode), width =       as.integer(width), 
   digits = as.integer(digits), format = as.character(format), 
   flag = as.character(flag), result = blank.chars(i.strlen), 
   PACKAGE = "base")
5: FUN(X[[1L]], ...)
4: lapply(data, formatC, digits = digits, format = "E")
3: write.datafile(lapply(data, formatC, digits = digits, format = "E"), 
   file.path(dir, data.file))
2: bugs.data(data, dir = getwd(), digits)
1: bugs(data = list(Y = Y, Nf = 5), 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 = "C://Tina/USB_Backup_042213 /winbugs14/WinBUGS14", 
   working.directory = NULL)

任何建议都会很有帮助。谢谢!

Answer 1

我也犯了同样的错误。在我的例子中，问题是Y是True/False格式，而不是0/1格式。