This code calculates the average and range of jellyfish density inside aggregations in Sinclair Inlet and Quartermaster Harbor.

Markdown file is here: https://rpubs.com/HailaSchultz/field_jelly-density

load packages

library(knitr)
library(dplyr)

Set up data

Read database into R. If current database changes, just change the path

Database <- read.csv("/Users/hailaschultz/Dropbox/Other studies/Aurelia project/Data Analysis/data/current_data/Final_Aurelia_Database_Jan11_2023.csv")

Subset data

Subset to field data

Field_stations<- subset(Database, Trial.Type=='Field')
Field_stations_IN<-subset(Field_stations,Location=="IN")

Rename variables

#convert sample year to a factor
Field_stations_IN$Sample.Year<- as.factor(Field_stations_IN$Sample.Year)
#remove commas from jelly density
Field_stations_IN$Jelly.Density <- as.numeric(gsub(",","",Field_stations_IN$Jelly.Density..g.m3.))

Summarize

#summarize by station
Jelly_density_sum<-Field_stations_IN %>%
  group_by(Site,Sample.Code) %>%
  summarise(
    Jelly.Density = mean(Jelly.Density))
#get mean, min,max
Jelly_density_sum<-Jelly_density_sum %>%
  group_by(Site) %>%
  summarise(
    mean = mean(Jelly.Density),
    min=min(Jelly.Density),
    max=max(Jelly.Density))

make table

knitr::kable(Jelly_density_sum)
Site mean min max
QM 497.8780 35.65 2305.88
SC 488.3207 48.03 1382.86