For the section:
# Median incubation time: med_incubation <- median(linelist$incubation, na.rm = TRUE)
Add documentation on ?median. Explain that my default, the function median does not remove missing values na.rm = FALSE. This would result that the median would return a NA value. Thus we need to specify na.rm = TRUE
For the section:
# Median incubation time: med_incubation <- median(linelist$incubation, na.rm = TRUE)Add documentation on ?median. Explain that my default, the function median does not remove missing values
na.rm = FALSE. This would result that the median would return a NA value. Thus we need to specifyna.rm = TRUE