SNWR_FishSAC_2022

Setup

knitr::opts_chunk$set(echo = TRUE)

library(ggplot2)
library(vegan)
library(ggrepel)
library(magrittr)

Create Species Accumulation Curves for each unit

setwd("C:/Users/16mfr/OneDrive/Desktop/Capstone/Data and Code - Copy/SAC")

mcsac = read.csv("2022_MCSAC.csv")
mcsac <- mcsac %>% replace(is.na(.), 0)
sp1mc = specaccum(mcsac, "random")
mcsacplot <- plot(sp1mc,ci.type="poly",ci.col="pink",col="blue",xlab="No. Sites",ylab="Avg. # Species", main="MC 2022", asp = 1)
boxplot(sp1mc,add=TRUE)

mcsacplot

## NULL

mnsac = read.csv("2022_MNSAC.csv")
mnsac <- mnsac %>% replace(is.na(.), 0)
sp1mn = specaccum(mnsac, "random")
mnsacplot <- plot(sp1mn,ci.type="poly",ci.col="yellow",col="blue",xlab="No. Sites",ylab="Avg. # Species", main="MN 2022", asp =1)
boxplot(sp1mn,add=TRUE)

mnsacplot

## NULL

mssac = read.csv("2022_MSSAC.csv")
mssac <- mssac %>% replace(is.na(.), 0)
sp1ms = specaccum(mssac, "random")
mssacplot <- plot(sp1ms,ci.type="poly",ci.col="green",col="blue",xlab="No. Sites",ylab="Avg. # Species", main="MS 2022", asp = 1)
boxplot(sp1ms,add=TRUE)

mssacplot

## NULL

p1asac = read.csv("2022_P1ASAC.csv")
p1asac <- p1asac %>% replace(is.na(.), 0)
sp1p1a = specaccum(p1asac, "random")
p1aplot <- plot(sp1p1a,ci.type="poly",ci.col="blue",col="black",xlab="No. Sites",ylab="Avg. # Species", main="P1A 2022", asp = 1)
boxplot(sp1p1a,add=TRUE)

p1aplot

## NULL

shrsac = read.csv("2022_SHRSAC.csv")
shrsac <- shrsac %>% replace(is.na(.), 0)
sp1shr = specaccum(shrsac, "random")
shrsacplot <- plot(sp1shr,ci.type="poly",ci.col="purple",col="blue",xlab="No. Sites",ylab="Avg. # Species", main="SHR 2022", asp = 1)
boxplot(sp1shr,add=TRUE)

shrsacplot

## NULL