我只有一点使用 Stan 的经验,并且正在尝试构建一个分层模型。我遇到了一个我以前从未见过的错误,这对我来说完全没有信息。有人知道这里发生了什么吗?我的代码和错误消息如下。
model_default_prior = "
data {
int<lower=0> N;
int<lower=0> n[N];
int<lower=0> y[N];
}
parameters {
real<lower=0> alpha;
real<lower=0> beta;
real<lower=0,upper=1> theta[N];
}
model {
// default prior
target += -5*log(alpha+beta)/2;
// implicit joint distributions
theta ~ beta(alpha,beta);
y ~ binomial(n,theta);
}
"
rats.dat = list(y = rats$y, n = rats$n, N = nrow(rats))
m = stan_model(model_code=model_default_prior)
r.rats = sampling(m, rats.dat, c("alpha","beta","theta"), iter=10000,
control = list(adapt_delta = 0.9))
"Error in sampler$call_sampler(args_list[[i]]) : " " c++ exception (unknown reason)"
error occurred during calling the sampler; sampling not done