Stan 2.10 syntax

stan-dev · Nov 3, 2016 · 2015e2a · 2015e2a
1 parent 8b7e50b
commit 2015e2a
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 17 deletions.
diff --git a/BPA/Ch.03/GLM_Binomial.stan b/BPA/Ch.03/GLM_Binomial.stan
@@ -8,7 +8,7 @@ data {
 transformed data {
   vector[nyears] year_squared;
 
-  year_squared <- year .* year;
+  year_squared = year .* year;
 }
 
 parameters {
@@ -21,9 +21,9 @@ transformed parameters {
   vector[nyears] logit_p;
 
   // Linear predictor
-  logit_p <- alpha +
-             beta1 * year +
-             beta2 * year_squared;
+  logit_p = alpha
+          + beta1 * year
+          + beta2 * year_squared;
 }
 
 model {
@@ -41,5 +41,5 @@ generated quantities {
   real<lower=0,upper=1> p[nyears];
 
   for (i in 1:nyears)
-    p[i] <- inv_logit(logit_p[i]);
+    p[i] = inv_logit(logit_p[i]);
 }
diff --git a/BPA/Ch.03/GLM_Poisson.stan b/BPA/Ch.03/GLM_Poisson.stan
@@ -8,8 +8,8 @@ transformed data {
   vector[n] year_squared;
   vector[n] year_cubed;
 
-  year_squared <- year .* year;
-  year_cubed <- year .* year .* year;
+  year_squared = year .* year;
+  year_cubed = year .* year .* year;
 }
 
 parameters {
@@ -22,18 +22,14 @@ parameters {
 transformed parameters {
   vector[n] log_lambda;
 
-  log_lambda <- alpha +
-                beta1 * year +
-                beta2 * year_squared +
-                beta3 * year_cubed;
+  log_lambda = alpha
+             + beta1 * year +
+             + beta2 * year_squared +
+             + beta3 * year_cubed;
 }
 
 model {
-  // Priors
-  alpha ~ uniform(-20, 20);
-  beta1 ~ uniform(-10, 10);
-  beta2 ~ uniform(-10, 10);
-  beta3 ~ uniform(-10, 10);
+  // Implicit uniform priors are used.
 
   // Likelihood
   C ~ poisson_log(log_lambda);
@@ -42,5 +38,5 @@ model {
 generated quantities {
   vector[n] lambda;
 
-  lambda <- exp(log_lambda);
+  lambda = exp(log_lambda);
 }