[DATALAD] Recorded changes

code/models/dev_siblings.stan

@@ -333,7 +333,7 @@ model {
 
         for (j in 1:n_classes) {
             mus[i,j] ~ exponential(1);
-            alphas[i,j] ~ lognormal(0, 1);
+            alphas[i,j] ~ exponential(0.1);
             // mus[i,j] ~ exponential(1);
             // alphas[i,j] ~ exponential(1);
             for (c in 1:n_groups) {
@@ -344,11 +344,11 @@ model {
 
     // speech rates
     mu_pop_level ~ exponential(4); // 250 vocs/hour
-    alpha_pop_level ~ gamma(8, 4); // sd = 0.35 x \mu
-    alpha_pop ~ gamma(10, 10);
+    alpha_pop_level ~ gamma(16, 4); // dispersion of corpora within population sd = 0.5 x \mu
+    alpha_pop ~ gamma(16, 4); // dispersion of dispersion of individuals within corpora
     for (i in 1:n_classes-1) {
-        alpha_corpus_level[1,i,:] ~ gamma(4, 4/alpha_pop[i]);
-        alpha_corpus_level[2,i,:] ~ gamma(4, 4/alpha_pop[i]);
+        alpha_corpus_level[1,i,:] ~ gamma(8, 8/alpha_pop[i]);
+        alpha_corpus_level[2,i,:] ~ gamma(8, 8/alpha_pop[i]);
         mu_corpus_level[i,:] ~ gamma(alpha_pop_level[i],alpha_pop_level[i]/mu_pop_level[i+1]);
     }
 

code/models/dev_siblings_effect.stan

@@ -343,8 +343,8 @@ model {
 
     // speech rates
     mu_pop_level ~ exponential(4); // 250 vocs/hour
-    alpha_pop_level ~ gamma(16, 4); // sd = 0.5 x \mu
-    alpha_pop ~ gamma(10, 10);
+    alpha_pop_level ~ gamma(16, 4); // dispersion of corpora within population sd = 0.5 x \mu
+    alpha_pop ~ gamma(16, 4); // dispersion of dispersion of individuals within corpora
     for (i in 1:n_classes-1) {
         alpha_corpus_level[1,i,:] ~ gamma(8, 8/alpha_pop[i]);
         alpha_corpus_level[2,i,:] ~ gamma(8, 8/alpha_pop[i]);