obs_lpmf.stan

Go to the documentation of this file.

 
real obs_lpmf(array[] int obs, vector log_exp_obs, array[] real phi,
                   int model_obs) {
  real tar = 0;
  if (model_obs == 0) {
    tar = poisson_log_lpmf(obs | log_exp_obs);
  }else if (model_obs == 1) {
    tar = neg_binomial_2_log_lpmf(obs | log_exp_obs, phi[1]);
  }else{
    int m = num_elements(log_exp_obs);
    vector[m] log_phi_nb1 = log_exp_obs + log(phi[1]);
    tar = neg_binomial_2_log_lpmf(obs | log_exp_obs, exp(log_phi_nb1));
  }
  return(tar);
}
 
real obs_lpmf(int obs, real log_exp_obs, array[] real phi,
         int model_obs) {
  real tar = 0;
  if (model_obs == 0) {
    tar = poisson_log_lpmf(obs | log_exp_obs);
  }else if (model_obs == 1) {
    tar = neg_binomial_2_log_lpmf(obs | log_exp_obs, phi[1]);
  }else{
    real log_phi_nb1 = log_exp_obs + log(phi[1]);
    tar = neg_binomial_2_log_lpmf(obs | log_exp_obs, exp(log_phi_nb1));
  }
  return(tar);
}