20array[]
int obs_rng(vector log_exp_obs, array[] real phi,
int model_obs) {
21 int n = num_elements(log_exp_obs);
24 pp = neg_binomial_2_log_rng(log_exp_obs, phi[1]);
26 pp = poisson_log_rng(log_exp_obs);
array[] int obs_rng(vector log_exp_obs, array[] real phi, int model_obs)