M Quillen Fur Seal Data Analysis

Read in Hg concentrations file from local directory

# read in file from local directory
dir <- "/Users/bmeyer/Google Drive/UAF Semesters/Summer 2020/M_Quillen_Draft/THg_BodyLocation_NFS.csv"
x <- read.csv(dir) 

Create Hg concentrations boxplot

# create boxplot of THg value distribution by body part
x %>%
  gather(body_part,value,RHIP,SAT,VHF) %>%
  filter(!is.na(value)) %>%
  ggplot(aes(body_part,value)) +
    geom_boxplot() +
  xlab("Body Location") +
  ylab("Total Mercury Concentratin (µg/g)") +
  ggtitle("Mercury Concentrations\nMeasured From Different Parts of The Body Of Northern Fur Seals")

Other exploratory data analyses (in progress – Based on the available data, these might be interesting questions):

• How do THg concentrations vary by year?
• Do C and N stable isotope values predict THg concentrations?
• Others?