changes to data and model

R/data.R

@@ -1419,9 +1419,9 @@
 "df_lagoE1"
 
 
-#' Working memory capacity scores and speed of the rapid automatized naming task from the exploratory data of \insertCite{NicenboimEtAlCogSci2018;textual}{bcogsci}
+#' Working memory capacity scores and speed of the rapid automatized naming task from the confirmatory data of \insertCite{NicenboimEtAlCogSci2018;textual}{bcogsci}
 #'
-#'Data from from the exploratory experiment  of \insertCite{NicenboimEtAlCogSci2018;textual}{bcogsci}.
+#'Data from from the confirmatory experiment  of \insertCite{NicenboimEtAlCogSci2018;textual}{bcogsci}.
 #'The data shows (1) partial-credit unit scores, which indicate the mean proportion of correctly recalled items within the sets of an operation span task, as an index of working memory capacity \insertCite{conway2005working}{bcogsci} together with their SE, and (2) the mean speeds in the rapid automatized naming task, RAN, \insertCite{DenklaRudel1976}{bcogsci} as a measure of reading fluency (in characters/second) together with their SE.
 #'
 #'

data-raw/data_repos.R

@@ -271,21 +271,21 @@ df_blp <- df_blp_complete %>% filter(subj <= 21) %>%
 
 #df_blp %>% filter(subj ==1, lex =="word")
 
-GET("https://osf.io/vgudp//?action=download",
-    write_disk("data-raw/data_repos/exp_subj.Rda", overwrite = TRUE),
+GET("https://osf.io/49q2t//?action=download",
+    write_disk("data-raw/data_repos/conf_subj.Rda", overwrite = TRUE),
     progress()
 )
 
-load("data-raw/data_repos/exp_subj.Rda")
-rs_expl <- summarise(group_by(subj_expl$ran,subj),
+load("data-raw/data_repos/conf_subj.Rda")
+rs <- summarise(group_by(subj_conf$ran,subj),
                                mean_rspeed = mean(1/rt),
                                se_rspeed = sd(1/rt)/sqrt(n()))
 
-pcu_expl <- summarise(group_by(subj_expl$ospan,subj),
+pcu <- summarise(group_by(subj_conf$ospan,subj),
                       mean_pcu= mean(memory / level),
                       se_pcu = sd(memory / level)/sqrt(n()))
 
-df_indiv <-  left_join(rs_expl, pcu_expl)
+df_indiv <-  left_join(rs, pcu)
 
 usethis::use_data(df_pupil,
                   df_pupil_complete,

inst/stan_models/me.stan

@@ -21,11 +21,11 @@ model {
   target += normal_lpdf(y | alpha + beta * x_true,
                         sqrt(square(SE_y) + square(sigma)));
   target += std_normal_lpdf(z);
-  target += normal_lpdf(alpha | 50, 50);
-  target += normal_lpdf(beta | 0, 10);
-  target += normal_lpdf(chi | 0, 10);
-  target += normal_lpdf(sigma| 0, 30)
-    - normal_lccdf(0 | 0, 30);
-  target += normal_lpdf(tau | 0, 30)
-    - normal_lccdf(0 | 0, 30);
+  target += normal_lpdf(alpha | 0, 0.5);
+  target += normal_lpdf(beta | 0, 0.5);
+  target += normal_lpdf(chi | 0, 0.5);
+  target += normal_lpdf(sigma| 0, 0.5)
+    - normal_lccdf(0 | 0, 0.5);
+  target += normal_lpdf(tau | 0, 0.5)
+    - normal_lccdf(0 | 0, 0.5);
 }