
```{r}
#| label = "{{DataName}}.lm"
{{DataName}}.lm = lm({{ResponseName}} ~ {{RowFactor}}+{{ColFactor}}+{{FactorName}}, data = {{DataName}})
{{DataName}}.lm |> anova()
```


```{r}
#| label = "plot{{DataName}}.lm",
#| fig.cap = "Residual analysis for {{DataName}}.lm"
library(ggfortify)
{{DataName}}.lm |> autoplot()
```


```{GetMeans}
library(emmeans)
emmeans({{DataName}}.lm, pairwise ~{{FactorName}})
```
