create linear models

Author: wlangera

source/dataset_bias_cv.Rmd

@@ -928,6 +928,20 @@ grouped_lm(
 )
 ```
 
+```{r}
+m <- lm(
+  rmse ~ evenness * rarity,
+  data = trend_dataset_rmse
+)
+
+# Diagnostics
+plot(m)
+```
+
+```{r}
+summary(m)
+```
+
 ## Trends in error: MRE
 
 We look at trends in CV error measures related to:
@@ -1088,6 +1102,20 @@ grouped_lm(
 )
 ```
 
+```{r}
+m2 <- lm(
+  mre ~ evenness * rarity,
+  data = trend_dataset_mre
+)
+
+# Diagnostics
+plot(m2)
+```
+
+```{r}
+summary(m2)
+```
+
 ```{r, echo=FALSE}
 col_rmse <- plot_grid(
   p_rmse_rarity +