robertness
6/25/2017 - 11:42 PM

binomial_stan

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]);
}
"