1. Download data

Downloaded to bio folder

2. Create new RMD file

Created as “Salinity Data Project”

3. Read data into R

data <- read.csv("DataSalinity_SFbar_FortPoint.csv")

summary(data)
##      Date             WaterTemp        Salinity       Location        
##  Length:14967       Min.   :11.61   Min.   :27.65   Length:14967      
##  Class :character   1st Qu.:14.15   1st Qu.:30.56   Class :character  
##  Mode  :character   Median :14.82   Median :31.94   Mode  :character  
##                     Mean   :14.67   Mean   :31.34                     
##                     3rd Qu.:15.20   3rd Qu.:32.12                     
##                     Max.   :17.37   Max.   :32.41                     
##                                     NA's   :9175

4. Make sure dates are recognized as dates

library(chron)

data$Date <- as.chron(data$Date)

summary(data)
##       Date                       WaterTemp        Salinity    
##  Min.   :(10/10/21 01:00:00)   Min.   :11.61   Min.   :27.65  
##  1st Qu.:(10/17/21 06:05:00)   1st Qu.:14.15   1st Qu.:30.56  
##  Median :(10/24/21 07:15:00)   Median :14.82   Median :31.94  
##  Mean   :(10/25/21 11:28:33)   Mean   :14.67   Mean   :31.34  
##  3rd Qu.:(11/02/21 09:37:30)   3rd Qu.:15.20   3rd Qu.:32.12  
##  Max.   :(11/11/21 12:55:00)   Max.   :17.37   Max.   :32.41  
##                                                NA's   :9175   
##    Location        
##  Length:14967      
##  Class :character  
##  Mode  :character  
##                    
##                    
##                    
##

5. Plot date vs salinity

plot(x=data$Date, y=data$Salinity, col="aquamarine3")

5. Add line for Oct 24th 2021

plot(x=data$Date, y=data$Salinity, col="aquamarine3")
abline(v=as.chron("10/24/2021"))

ggplot facet_wrap

library(ggplot2)

myplot <- ggplot(data, aes(Date, WaterTemp))
myplot + geom_point(col="deepskyblue4")+
  scale_x_chron()+
  geom_vline(xintercept = as.chron("10/24/2021"),linetype = 2,color = "coral1", size = 1) +
  facet_wrap(~Location)

ggplot color for location

library(ggplot2)

myplot <- ggplot(data, aes(Date, WaterTemp))
myplot + geom_point(aes(col=Location), alpha=0.1) +
  scale_x_chron()+
  geom_vline(xintercept = as.chron("10/24/2021"),linetype = 2,color = "darkorchid4", size = 1)