binomial_stan
stan_str <- " data { int N; int D; int Y[N]; int M[N]; vector[D] X[N]; } parameters { matrix[1, D] beta; } model { for (n in 1:N) dead[n] ~ binomial_logit(M[n], beta * X[n]); } "