fixed file names so they work with wiki links

Author: PeterLi2016

ARM/Ch.10/10.6_IVinaRegressionFramework.R

@@ -44,7 +44,7 @@ watched.hat <- beta.mean3a[1] + beta.mean3a[2] * encouraged + beta.mean3a[3] * p
 ## (sesame_multi_preds_3b.stan)
 ## lm (y ~ watched.hat + pretest + as.factor(site) + setting)
 dataList.4 <- list(N=length(watched.hat), watched=y,encouraged=watched.hat,pretest=pretest, site=site,setting=setting)
-sesame_multi_pred_3b.sf1 <- stan(file='sesame_multi_preds_3a.stan',
+sesame_multi_pred_3b.sf1 <- stan(file='sesame_multi_preds_3b.stan',
                                  data=dataList.4,
                                  iter=1000, chains=4)
 print(sesame_multi_pred_3b.sf1)

ARM/Ch.10/sesame_mult_preds_3b.stan ARM/Ch.10/sesame_multi_preds_3b.stanARM/Ch.10/sesame_mult_preds_3b.stan renamed to ARM/Ch.10/sesame_multi_preds_3b.stan

@@ -1,6 +1,17 @@
 library(rstan)
 library(ggplot2)
 
+### Data
+
+source("kidiq.data.R", echo = TRUE)
+
+### Model: kid_score ~ mom_hs + mom_iq
+data.list <- c("N", "kid_score", "mom_hs", "mom_iq", "mom_hs_new", "mom_iq_new")
+kidiq_prediction.sf <- stan(file='kidiq_prediction.stan', data=data.list,
+                            iter=1000, chains=4)
+print(kidiq_prediction.sf, pars = c("kid_score_pred"), prob = c(0.025, 0.975))
+
+
 ### Data
 
 source("kids_before1987.data.R", echo = TRUE)