polish up index

Author: probably-jaden

ch7ModelComparison.qmd

@@ -173,17 +173,22 @@ b7.0 <- brm(data = speciesBrain,
             prior = c(
               prior(normal(0.5, 1), class = Intercept),
               prior(normal(0, 1), class = b),
-              prior(exponential(1), class = sigma)
+              prior(lognormal(0,1), class = sigma)
             ),
             iter = 2000, warmup = 500, seed = 4, cores = 4,
-            backend = "cmdstanr", silent = 2, file = "fits/b07.0.4"
+            backend = "cmdstanr", silent = 2, file = "fits/b07.0.5"
             )
 
 b7.0_sim <- as_tibble(b7.0) %>% 
   mutate(simMass = seq(from = -3, to = 3, length.out = n()),
          simBrainEst = Intercept + (b_mass_z * simMass),
          simBrain = rnorm(n(), simBrainEst, sd = sigma))
 
+```
+
+```{r}
+#| fig-width: 6
+#| fig-height: 4
 
 ggplot() +
   stat_density_2d(data = b7.0_sim, 
@@ -192,9 +197,10 @@ ggplot() +
   scale_fill_viridis_c(option = "magma") +
   geom_point(data = speciesBrain,
              aes(x = mass_z, y = brain_z), 
-             shape = 21, color = "white", fill = "black", lwd = 3, alpha = .8)+
+             shape = 21, color = "white", fill = "black", lwd = 3, alpha = 1, size =2.5)+
   labs(title = "brain ~ mass", subtitle = "Mu estimate") +
   theme_minimal()+
+  ggrepel::geom_text_repel(data = speciesBrain, aes(x = mass_z, y = brain_z, label = species), color = "white", size = 3)+
   # xlim(c(-3,3))+
   # ylim(c(-3,3))+
   scale_x_continuous("body mass (kg)",
@@ -212,3 +218,21 @@ ggplot() +
 
 ```
 
+```{r}
+mass_brain <- speciesBrain %>% 
+  rename("mass_obs" = mass_z,
+         "brain_obs" = brain_z) %>%
+  dplyr::select(c(mass_obs, brain_obs))
+
+
+as_tibble(b7.0) %>%
+  # Add row numbers to track individual draws
+  mutate(draw_id = row_number()) %>%
+  # Create all combinations of draws and new X values
+  crossing(mass_brain) %>% 
+  mutate(brain_mu_pred = b_Intercept + (b_mass_z * brain_obs),
+         brain_mu_resid= brain_obs - brain_pred_mu,
+         brain_mu_
+         )
+```
+

index.qmd

@@ -1,17 +1,18 @@
 ---
-title: "Statistical Rethinking Homepage"
+title: "Statistical Rethinking"
 author: "Jaden Earl"
 format: html
 editor: visual
 ---
 
-# Welcome
+[Statistical Rethinking](https://oceanrep.geomar.de/id/eprint/55819/1/Statistical%20Rethinking%202nd%20Edition.pdf) has been a fun unlearning of the frequentist techniques that I practised in my undergrad. I appreciated it so much that I wanted to document my current statistical skillset using his textbook as the structure and basebone. Richard McElreath also has a great series on [youtube](https://www.youtube.com/watch?v=FdnMWdICdRs&list=PLDcUM9US4XdPz-KxHM4XHt7uUVGWWVSus) going through each chapter of his book that I highly reccommend as well.
 
-This site showcases projects from my **Statistical Rethinking** work in R and Quarto.
+## 📈 [Linear Regression](ch4_linear_regression_polished.qmd)
 
--   📈 [Linear Regression](ch4_linear_regression_polished.qmd)\
--   🔀 [Multi Regression](ch5_multilinear.qmd)\
--    ➡️ [Causal Paths](ch6_DAG.qmd)
--    🎯️ [Model Scoring](ch7ModelComparison.qmd)
+## 📊 [Multi Regression](ch5_multilinear.qmd) 
+
+## 🔀️ [Causal Paths](ch6_DAG.qmd)
+
+## 🎯️ [Model Scoring](ch7ModelComparison.qmd)
 
 Source code is on [GitHub](https://github.com/probably-jaden/Statistical-Rethinking).