Я пишу модель, используя функцию ulam в пакете переосмысления в rstan. Список данных, который я предоставляю, состоит из цифр c (is.numeri c возвращает true для каждой переменной) сам по себе, но не в целом. Как я могу сделать так, чтобы модель могла работать? Вот список данных и код модели, а также полное сообщение об ошибке:
Список данных:
datBIG <- list(
R = d$R,
P = d$P,
I = d$I,
N = d$N,
M = d$M,
B = d$B,
AID = d$AID)
Модель:
mBIG<-ulam(alist(
V ~ dnorm(mu, sigma),
mu <- a[AID] + bR*R + bP*P + bI*I + bN*N + bM*M + bB*B,
a[AID] ~ dnorm(abar, sigma_a),
abar ~ dnorm(0, 0.2),
bR ~ dnorm(0, 0.2),
bP ~ dnorm(0, 0.2),
bI ~ dnorm(0, 0.2),
bN ~ dnorm(0, 0.2),
bM ~ dnorm(0, 0.2),
bB ~ dnorm(0, 0.2),
sigma_a ~ dexp(1),
sigma ~ dexp(1)
), data=datBIG, chains=1, cores=1, log_lik=TRUE)
И сообщение об ошибке:
Ошибка в validObject (.Object): недопустимый объект класса «ulam»: недопустимый объект для слота «coef» в классе «ulam»: получил класс «NULL», должен быть или расширять класс «цифр» c "