# 線形回帰の基本プロット
ここではbaseパッケージを用いた線形回帰のための二次元散布図と回帰直線を描いてみます．

## データの読み込み
```{r}
library(vegan)
phyto_metadata <- readRDS("phyto_metadata.obj")
species_ryuko_data <- readRDS("phyto_ryuko_data.obj")
species_richness <- apply(species_ryuko_data > 0, 1, sum)
total_abundance <- apply(species_ryuko_data, 1, sum)
```
## 二次元散布図＋回帰直線
植物プランクトンのデータを使った場合
```{r}
#Classical plot (base)
model01 <- lm(species_richness ~ phyto_metadata$temp)
plot(
  species_richness ~ phyto_metadata$temp, 
  type = "p", 
  cex = 3,
  xlab = "temperature",
  ylab = "species richness"
)
abline(model01,col = 4)　#回帰直線
```

# 線形回帰のggplot
```{r}

```

