@@ -39,7 +39,7 @@ input_data <- list(
3939pars_poisson <- c(' alpha_national' , ' pop_post_pred' )
4040
4141fit_poisson <- stan(
42- file = here(' tutorials' , ' refresher' ,' poisson.stan' ),
42+ file = here(' tutorials' , ' refresher' , ' section1 ' , ' poisson.stan' ),
4343 data = input_data ,
4444 iter = iter + warmup ,
4545 warmup = warmup ,
@@ -103,7 +103,7 @@ ggplot(comp_pop, aes(x=value, fill=model))+
103103pars_lognormal <- c(' alpha_national' , ' sigma' , ' pop_post_pred' )
104104
105105fit_poisson_lognormal <- stan(
106- file = here(' tutorials' , ' refresher' ,' poisson_lognormal.stan' ),
106+ file = here(' tutorials' , ' refresher' ,' section2 ' , ' poisson_lognormal.stan' ),
107107 data = input_data ,
108108 iter = iter + warmup ,
109109 warmup = warmup ,
@@ -180,7 +180,7 @@ pars_lognormal_independent <- c('pop_post_pred', 'alpha_settlement', 'sigma')
180180
181181
182182fit_poisson_lognormal_independent <- stan(
183- file = here(' tutorials' , ' refresher' ,' poisson_independent.stan' ),
183+ file = here(' tutorials' , ' refresher' ,' section2 ' , ' poisson_independent.stan' ),
184184 data = input_data ,
185185 iter = iter + warmup ,
186186 warmup = warmup ,
@@ -237,7 +237,7 @@ ggplotly(ggplot(pop_posterior_lognormal_independent, aes(x=value, fill=source, a
237237pars_lognormal_hierarchical_cp <- c(' alpha_national' ,' alpha_settlement' , ' u_alpha_settlement' , ' sigma' , ' pop_post_pred' )
238238
239239fit_poisson_lognormal_hierarchical_cp <- stan(
240- file = here(' tutorials' , ' refresher' ,' poisson_hierarchical_cp.stan' ),
240+ file = here(' tutorials' , ' refresher' ,' section2 ' , ' poisson_hierarchical_cp.stan' ),
241241 data = input_data ,
242242 iter = iter + warmup ,
243243 warmup = warmup ,
@@ -350,7 +350,7 @@ pars_lognormal_hierarchical_var <- c('alpha_national','alpha_settlement',
350350 ' sigma_settlement' , ' pop_post_pred' )
351351
352352fit_poisson_lognormal_hierarchical_var <- stan(
353- file = here(' tutorials' , ' refresher' ,' poisson_hierarchical_variance.stan' ),
353+ file = here(' tutorials' , ' refresher' ,' section3 ' , ' poisson_hierarchical_variance.stan' ),
354354 data = input_data ,
355355 iter = iter + warmup ,
356356 warmup = warmup ,
@@ -456,7 +456,7 @@ pars_lognormal_hierarchical_nclp <- c('alpha_national','delta_settlement', 'u_de
456456 ' sigma' , ' pop_post_pred' )
457457
458458fit_poisson_lognormal_hierarchical_nclp <- stan(
459- file = here(' tutorials' , ' refresher' ,' poisson_hierarchical_nclp.stan' ),
459+ file = here(' tutorials' , ' refresher' ,' section3 ' , ' poisson_hierarchical_nclp.stan' ),
460460 data = input_data ,
461461 iter = iter + warmup ,
462462 warmup = warmup ,
@@ -540,7 +540,7 @@ pars_lognormal_hierarchical_nclsp <- c('alpha_national','u_delta_settlement',
540540 ' sigma' , ' pop_post_pred' )
541541
542542fit_poisson_lognormal_hierarchical_nclsp <- stan(
543- file = here(' tutorials' , ' refresher' ,' poisson_hierarchical_nclsp.stan' ),
543+ file = here(' tutorials' , ' refresher' ,' section3 ' , ' poisson_hierarchical_nclsp.stan' ),
544544 data = input_data ,
545545 iter = iter + warmup ,
546546 warmup = warmup ,
0 commit comments