I’ve succesfully complied two simpler models, but the following and other models with similar complexity fail to compile
writeLines(readLines("binomb.stan"))
data {
int<lower=0> N;
int<lower=0> y;
}
parameters {
real alpha;
}
transformed parameters {
real theta;
theta = inv_logit(alpha);
}
model {
alpha ~ normal(0,1.5);
y ~ binomial_logit(N,alpha);
}
d_bin <- list(N = 100, y = 70)
fit_bin <- stan(file = 'binomb.stan', data = d_bin)