---
title: "227.212 Biostatistics Assignment"
author: "CHANGE ME!!!"
---

```{r setup, message=FALSE, echo=FALSE}
# Don't change this chunk unless you know what you're doing :)
# Make sure you've installed each of these packages (you'll be prompted for this on first load!)
library(ggplot2)
library(broom)
library(dplyr)
all <- read.csv("https://www.massey.ac.nz/~jcmarsha/227212/data/dairy.csv")
all$Month <- factor(all$Month)
all$MonthLabel <- factor(all$Month, labels=month.abb[c(6:12,1:5)])
dairy <- na.omit(all)
```

## Exercise 1: Milk through time

### 1. Milk per cow per day by Month

```{r}
ggplot(data=dairy) +
  geom_boxplot(mapping=aes(x=Month, y=Milk))
```

### 2. Milk per cow per day by Year

```{r}

```

### 3. Outliers in May

```{r}

```

### 4. Milk per cow per day by region

```{r}

```

## Exercise 2: Somatic cell counts
  
### 1. Distribution of somatic cell counts

```{r}

```

### 2. Distribution of log somatic cell counts

```{r}

```

### 3. Confidence intervals for average somatic cell count by island

```{r}

```

### 4. Somatic cell counts by month for North and South islands

```{r}

```

## Exercise 3: Relationships with milk volume

### 1. Milk fat versus milk volume

```{r}

```

### 2. Average milk volume by island

```{r}

```
