---
title: "Calculation of growing degree days (GDD) in the pollen package"
author: Jakub Nowosad
date: "`r Sys.Date()`"
output: rmarkdown::html_vignette
vignette: >
  %\VignetteIndexEntry{Calculation of growing degree days (GDD) in the pollen package}
  %\VignetteEngine{knitr::rmarkdown}
  %\VignetteEncoding{UTF-8}
bibliography: "`r system.file('refs.bib', package = 'pollen')`"
nocite: '@Baskerville1969'
---

```{r, include = FALSE}
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
```

Growing degree days (GDD) are a measure of heat accumulation used to predict plant and animal development rates.
It can be, however, calculated using several approaches.
The **pollen** package allows for three different versions of the GDD calculations.
The goal of this vignette is to describe and show examples of calculating GDD using different methods.

Let's start with attaching a set of package used in this vignette:

```{r setup, message=FALSE}
library(pollen)
library(ggplot2)
library(tidyr)
```

Next, we will use the `gdd_data` dataset - a data frame with three columns:

- `day` - this could be a date in real data
- `tmax` - this value would represent a maximum temperature in real data
- `tmin` - this value would represent a minimum temperature in real data

```{r}
data("gdd_data", package = "pollen")
head(gdd_data)
```

We can vizualize the `gdd_data` dataset using the **ggplot2** package:

```{r}
df_plot1 <- pivot_longer(gdd_data, tmax:tmin)
p1 <- ggplot(df_plot1) +
  geom_line(aes(day, value, color = name))
p1
```

The **pollen** package allows for calculations of growing degree days (GDD) using the `gdd()` function. 
This function accepts up to five arguments:

- `tmax` - daily maximum temperature
- `tmin` - daily minimum temperature
- `tbase` - base temperature
- `tbase_max` - maximum base temperature
- `type` - type of the GDD calculations. 
Either "B", "C", or "D". 
The default is "C".

The last argument is inspired by the article by Baskerville and Emin (1969) (see Figure 1 in the mentioned paper).

- Type `"B"` - The heat units are calculated based on the difference between the mean daily temperature and the threshold (`tbase`). 
In the case when the value of `tmin` is lower than `tbase`, then it is replaced by `tbase`
- Type `"C"` - same as type `"B"` and when the value of `tmax` is larger than `tbase_max`, then it is replaced by `tbase_max`
- Type `"D"`- same as type `"B"` and when the value of `tmax` is larger than `tbase_max`, then no heat units are added

```{r}
gdd_data$type_b <- gdd(tmax = gdd_data$tmax, tmin = gdd_data$tmin, 
                       tbase = 5, type = "B")
gdd_data$type_c <- gdd(tmax = gdd_data$tmax, tmin = gdd_data$tmin, 
                       tbase = 5, tbase_max = 20, type = "C")
gdd_data$type_d <- gdd(tmax = gdd_data$tmax, tmin = gdd_data$tmin, 
                       tbase = 5, tbase_max = 20, type = "D")
head(gdd_data)
```

Finally, let's compare these three types:

```{r}
df_plot2 <- pivot_longer(gdd_data, type_b:type_d)
p2 <- ggplot(df_plot2) +
  geom_line(aes(day, value, color = name))
p2
```

The above figure clearly shows that the selection of the calculation method (`type`) can make a large difference in the output values.
Therefore, it is crucial to select the `type` and other parameters (`tbase` and `tbase_max`) that are appropriate for the studied phase in the plant or animal development. 

# References