Agregar sims de coeficientes

felipegonzalez · Apr 23, 2024 · 0c0f1f5 · 0c0f1f5
1 parent bd43e30
commit 0c0f1f5
Show file tree

Hide file tree

Showing 3 changed files with 12 additions and 1 deletion.
diff --git a/notas/src/bangladesh-3.stan b/notas/src/bangladesh-3.stan
@@ -43,5 +43,8 @@ generated quantities {
     prob_distrito_urbano[i] = inv_logit(alpha[i] + beta[i]);
     prob_distrito_rural[i] = inv_logit(alpha[i]);
   }
-
+  // Simular de a priori poblacional
+  vector[2] beta_sim;
+  beta_sim[1] = normal_rng(alpha_bar, sigma_alpha);
+  beta_sim[2] = normal_rng(beta_bar, sigma_beta);
 }
diff --git a/notas/src/bangladesh-4.stan b/notas/src/bangladesh-4.stan
@@ -46,4 +46,7 @@ generated quantities {
     prob_distrito_urbano[i] = inv_logit(beta[i][1] + beta[i][2]);
     prob_distrito_rural[i] = inv_logit(beta[i][1]);
   }
+    // Simular de a priori poblacional
+  vector[2] beta_sim;
+  beta_sim = multi_normal_rng(beta_bar, Sigma);
 }
diff --git a/notas/src/bangladesh-5.stan b/notas/src/bangladesh-5.stan
@@ -56,4 +56,9 @@ generated quantities {
     prob_distrito_urbano[i] = inv_logit(beta[1,i] + beta[2,i]);
     prob_distrito_rural[i] = inv_logit(beta[1,i]);
   }
+
+  // Simular de a priori poblacional
+  vector[2] beta_sim;
+  beta_sim = multi_normal_rng(beta_bar, Sigma);
+
 }