renamed stan files

Author: PeterLi2016

bugs_examples/vol1/litter/litter.stan

@@ -6,12 +6,17 @@ data {
 }
 parameters {
   matrix<lower=0,upper=1>[G,N] p;
-  vector<lower=0.1>[G] a;
-  vector<lower=0.1>[G] b;
+  vector<lower=0,upper=1>[G] mu;
+  vector<lower=0.1>[G] a_plus_b;
+}
+transformed parameters {
+  vector[G] a;
+  vector[G] b;
+  a <- mu .* a_plus_b;
+  b <- (1 - mu) .* a_plus_b;
 }
 model {
-  a ~ gamma(1,0.001);
-  b ~ gamma(1,0.001);
+  a_plus_b ~ pareto(0.1,1.5);
   for (g in 1:G) {
     for (i in 1:N) { 
       p[g,i] ~ beta(a[g],b[g]);
@@ -20,10 +25,7 @@ model {
   }      
 }
 generated quantities {
-  vector<lower=0,upper=1>[G] mu;
   vector<lower=0>[G] theta;
-  for (g in 1:G)
-    mu[g] <- a[g] / (a[g] + b[g]);
   for (g in 1:G)
     theta[g] <- 1 / (a[g] + b[g]);
-}
+}

…examples/vol1/litter/litter_reparam.stan …amples/vol1/litter/litter_old_param.stanbugs_examples/vol1/litter/litter_reparam.stan renamed to bugs_examples/vol1/litter/litter_old_param.stan bugs_examples/vol1/litter/litter_reparam.stan renamed to bugs_examples/vol1/litter/litter_old_param.stan

@@ -6,17 +6,12 @@ data {
 }
 parameters {
   matrix<lower=0,upper=1>[G,N] p;
-  vector<lower=0,upper=1>[G] mu;
-  vector<lower=0.1>[G] a_plus_b;
-}
-transformed parameters {
-  vector[G] a;
-  vector[G] b;
-  a <- mu .* a_plus_b;
-  b <- (1 - mu) .* a_plus_b;
+  vector<lower=0.1>[G] a;
+  vector<lower=0.1>[G] b;
 }
 model {
-  a_plus_b ~ pareto(0.1,1.5);
+  a ~ gamma(1,0.001);
+  b ~ gamma(1,0.001);
   for (g in 1:G) {
     for (i in 1:N) { 
       p[g,i] ~ beta(a[g],b[g]);
@@ -25,7 +20,10 @@ model {
   }      
 }
 generated quantities {
+  vector<lower=0,upper=1>[G] mu;
   vector<lower=0>[G] theta;
+  for (g in 1:G)
+    mu[g] <- a[g] / (a[g] + b[g]);
   for (g in 1:G)
     theta[g] <- 1 / (a[g] + b[g]);
 }