behavior_model_priors_nonlinear.stan 702 B

12345678910111213141516171819
  1. alpha_child_level ~ gamma(4,1);
  2. mu_pop_level ~ exponential(4); // 250 vocs/hour
  3. alpha_pop_level ~ gamma(8, 4); // sd = 0.35 x \mu
  4. //alpha_pop ~ gamma(10, 10);
  5. for (i in 1:n_classes-1) {
  6. //alpha_corpus_level[1,i] ~ gamma(4, 4/alpha_pop[i]);
  7. //alpha_corpus_level[2,i] ~ gamma(4, 4/alpha_pop[i]);
  8. alpha_corpus_level[1,i] ~ gamma(4, 2);
  9. alpha_corpus_level[2,i] ~ gamma(4, 2);
  10. mu_corpus_level[i,:] ~ gamma(alpha_pop_level[i],alpha_pop_level[i]/mu_pop_level[i+1]);
  11. }
  12. has_siblings ~ binomial(has_siblings+no_siblings, p_sib);
  13. p_sib ~ uniform(0, 1);
  14. beta_sib_och ~ normal(0, 1);
  15. beta_sib_adu ~ normal(0, 1);
  16. d_beta_age ~ normal(0, 0.05); // 5% variation per month
  17. beta_dev ~ normal(0, 1);