Changed priors for generic weakly informative priors

BenMerSci · Apr 1, 2024 · e12d98d · e12d98d
1 parent b9e322d
commit e12d98d
Show file tree

Hide file tree

Showing 5 changed files with 16 additions and 16 deletions.
diff --git a/stan/08_stan_model0.stan b/stan/08_stan_model0.stan
@@ -15,8 +15,8 @@ parameters {
 model {
 
   // Prios:
-  Intercept ~ normal(0,1.5);
-  sigma ~ exponential(5);
+  Intercept ~ normal(0, 1);
+  sigma ~ exponential(3);
 
 
   log_biomass_flow ~ normal(Intercept, sigma);

diff --git a/stan/08_stan_model1.stan b/stan/08_stan_model1.stan
@@ -23,8 +23,8 @@ model {
 vector[n] mu_flow;
 
   // Priors:
-  alpha ~ normal(-4, 2);
-  sigma ~ exponential(5);
+  alpha ~ normal(0, 1);
+  sigma ~ exponential(3);
 
   // Likelihood:
   // Computing target's mean

diff --git a/stan/08_stan_model2.stan b/stan/08_stan_model2.stan
@@ -26,10 +26,10 @@ model {
   vector[n] mu_flow;
 
   // Priors:
-  mu_alpha ~ normal(-4, 2);
+  mu_alpha ~ normal(0, 1);
   alpha ~ normal(mu_alpha, sd_alpha);
   sd_alpha ~ exponential(3);
-  sigma ~ exponential(5);
+  sigma ~ exponential(3);
 
   // Computing target's mean
    mu_flow = alpha[pred_id] + log_biomass_prey + log_abundance_predator;

diff --git a/stan/08_stan_model3.stan b/stan/08_stan_model3.stan
@@ -30,13 +30,13 @@ model {
   vector[n] pred_factor;
 
   // Priors:
-  mu_alpha ~ normal(-5,2);
+  mu_alpha ~ normal(0, 1);
   alpha ~ normal(mu_alpha, sd_alpha);
-  sd_alpha ~ exponential(1);
-  mu_ht ~ normal(-3,2);
+  sd_alpha ~ exponential(3);
+  mu_ht ~ normal(0, 1);
   ht ~ normal(mu_ht, sd_ht);
-  sd_ht ~ exponential(1);
-  sigma ~ exponential(5);
+  sd_ht ~ exponential(3);
+  sigma ~ exponential(3);
 
 
   // Computing predators part for the numerator and denominator

diff --git a/stan/08_stan_model4.stan b/stan/08_stan_model4.stan
@@ -32,13 +32,13 @@ model {
   vector[n] pred_factor;
 
   // Priors:
-  mu_alpha ~ normal(-5,2);
+  mu_alpha ~ normal(0, 1);
   alpha ~ normal(mu_alpha, sd_alpha);
-  sd_alpha ~ exponential(1);
-  mu_ht ~ normal(-3,2);
+  sd_alpha ~ exponential(3);
+  mu_ht ~ normal(0, 1);
   ht ~ normal(mu_ht, sd_ht);
-  sd_ht ~ exponential(1);
-  sigma ~ exponential(5);
+  sd_ht ~ exponential(3);
+  sigma ~ exponential(3);
 
   // Computing predators part for the numerator and denominator