obs_rng.stan

Go to the documentation of this file.

 
array[] int obs_rng(vector log_exp_obs, array[] real phi, int model_obs) {
  int n = num_elements(log_exp_obs);
  array[n] int pp;
  if (model_obs == 0) {
    pp = poisson_log_rng(log_exp_obs);
  }else if (model_obs == 1) {
    pp = neg_binomial_2_log_rng(log_exp_obs, phi[1]);
  }else{
    int m = num_elements(log_exp_obs);
    vector[m] log_phi_nb1 = log_exp_obs + log(phi[1]);
    pp = neg_binomial_2_log_rng(log_exp_obs, exp(log_phi_nb1));
  }
  return(pp);
}