expected_obs.stan

Go to the documentation of this file.

 
vector expected_obs(real tar_obs, vector lh, int l, int ref_as_p, int agg_probs, matrix agg_indicator) {
  vector[l] p;
  if (ref_as_p == 1) {
    p = lh;
  }else{
    profile("model_likelihood_expected_obs_inv_logit") {
    p = inv_logit(lh);
    }
    profile("model_likelihood_expected_obs_hazard_to_prob") {
    p = hazard_to_log_prob(p, l);
    }
  }
  if (agg_probs == 1) {
    p = log(agg_indicator * exp(p));
  }
  return(tar_obs + p);
}