model_comparison_trees.Rmd

---
title: 'Model comparison'
author: 'Francisco Rodríguez-Sánchez'
institute: 'https://frodriguezsanchez.net'
aspectratio: 43  # use 169 for wide format
fontsize: 10pt
output: 
  binb::metropolis:
    keep_tex: no
    incremental: yes
    fig_caption: no
    pandoc_args: ['--lua-filter=hideslide.lua']
urlcolor: blue
linkcolor: blue
header-includes:
  - \definecolor{shadecolor}{RGB}{230,230,230}
  - \setbeamercolor{frametitle}{bg=gray}
---


```{r knitr_setup, include=FALSE, cache=FALSE}

library('knitr')

### Chunk options ###

## Text results
opts_chunk$set(echo = TRUE, warning = FALSE, message = FALSE, size = 'tiny')

## Code decoration
opts_chunk$set(tidy = FALSE, comment = NA, highlight = TRUE, prompt = FALSE, crop = TRUE)

# ## Cache
# opts_chunk$set(cache = TRUE, cache.path = 'knitr_output/cache/')

# ## Plots
# opts_chunk$set(fig.path = 'knitr_output/figures/')
opts_chunk$set(fig.align = 'center', out.width = '90%')

### Hooks ###
## Crop plot margins
knit_hooks$set(crop = hook_pdfcrop)

## Reduce font size
## use tinycode = TRUE as chunk option to reduce code font size
# see http://stackoverflow.com/a/39961605
knit_hooks$set(tinycode = function(before, options, envir) {
  if (before) return(paste0('\n \\', options$size, '\n\n'))
  else return('\n\n \\normalsize \n')
  })

```


## Trees dataset

```{r}
trees <- read.csv('data/trees.csv')
head(trees)
```

```{r echo=FALSE}
trees$site <- as.factor(trees$site)
```

## Four models

```{r}
m1 <- lm(height ~ dbh, data = trees)
```

```{r}
m2 <- lm(height ~ sex, data = trees)
```

```{r}
m3 <- lm(height ~ site, data = trees)
```

```{r}
m4 <- lm(height ~ site*dbh, data = trees)
```


## Compare model performance

\scriptsize

```{r}
library('performance')
compare_performance(m1, m2, m3, m4)
```


## Compare model performance

```{r out.width='70%'}
library('see')
plot(compare_performance(m1, m2, m3, m4))
```


## Compare parameters

\scriptsize

```{r}
library('parameters')
compare_parameters(m3, m4)
```


## Compare parameters

```{r}
library('parameters')
plot(compare_parameters(m3, m4))
```