Ошибка в коде JAGS при использовании пакета R2jags
Я пытаюсь сделать байесовскую полиномиальную логистическую регрессию с помощью пакета R2jags на R. Так я объявляю свои данные...
library(R2jags)
library(mlbench)
data(Glass)
attach(Glass)
n <- as.double(nrow(Glass)) # Sample size
p <- as.double(ncol(Glass)) # Regressors + 1
M <- as.double(length(levels(Type))) # Number of classes
y <- cbind(as.numeric(Type==1), as.numeric(Type==2), as.numeric(Type==3),
as.numeric(Type==5), as.numeric(Type==6), as.numeric(Type==7)) # Response
X <- as.matrix(Glass[,-10]) # Regressors + 1 Matrix
DataList <- list('n', 'p', 'M', 'y', 'X')
Par <- 'B'
IniVal <- function(){list('B' = rbind( rep(0,p) , matrix(rnorm((M-1)*p),
ncol = 10) ))}
output <- jags2(data = DataList, inits = IniVal, parameters.to.save = Par,
n.chains = 1, n.iter = 200, working.directory = path,
clearWD=TRUE, model.file="mod.txt")
Когда я запускаю этот код (путь указывает на каталог, где находится файл mod.txt), появляется следующее сообщение об ошибке:
Welcome to JAGS 4.2.0 on Sat Nov 25 16:51:30 2017
JAGS is free software and comes with ABSOLUTELY NO WARRANTY
Loading module: basemod: ok
Loading module: bugs: ok
Deleting model
Clearing data table
Loading module: dic: ok
Reading data file jagsdata.txt
syntax error, unexpected LIST, expecting DOUBLE or NA or ASINTEGER or 'c'
Error in file(file, "rt") : cannot open the connection
In addition: Warning message:
In file(file, "rt") :
cannot open file 'CODAindex.txt': No such file or directory
Пожалуйста, помогите мне решить это. Я понятия не имею, что означает эта ошибка. Кроме того, я подозреваю, что есть более короткий и эффективный способ написания модели. Спасибо!