// bioassay.stan // // An example in BDA3, Ch. 3 // // 2015-06-07 Yuki Yanai data { int J; int y[J]; int n[J]; real x[J]; } parameters { real alpha; real beta; } transformed parameters { real theta[J]; for (j in 1:J) { theta[j] <- 1 / (1 + exp(-alpha - beta * x[j])); } } model { y ~ binomial(n, theta); }