Grey seal pup dive behavior. Pup01; PTT 238272

Michelle Shero, November 18, 2023

Pup01

PTT 238272

Look at the ‘Histos’ File

setwd("~/Desktop/grey seal pup dive behavior/238272")
histos<- read.csv("238272-Histos.csv")
attach(histos)

The following objects are masked from histos (pos = 3):

    BadTherm, Bin1, Bin10, Bin11, Bin12, Bin13, Bin14, Bin15, Bin16, Bin17, Bin18, Bin19, Bin2, Bin20, Bin21, Bin22, Bin23, Bin24, Bin25, Bin26, Bin27, Bin28, Bin29,
    Bin3, Bin30, Bin31, Bin32, Bin33, Bin34, Bin35, Bin36, Bin37, Bin38, Bin39, Bin4, Bin40, Bin41, Bin42, Bin43, Bin44, Bin45, Bin46, Bin47, Bin48, Bin49, Bin5,
    Bin50, Bin51, Bin52, Bin53, Bin54, Bin55, Bin56, Bin57, Bin58, Bin59, Bin6, Bin60, Bin61, Bin62, Bin63, Bin64, Bin65, Bin66, Bin67, Bin68, Bin69, Bin7, Bin70,
    Bin71, Bin72, Bin8, Bin9, Count, Date, DeployID, DepthSensor, DOY, HistType, Hour, Instr, Latitude, LocationQuality, Longitude, NumBins, Ptt, Source, Sum,
    Time.Offset

The following objects are masked from histos (pos = 4):

    BadTherm, Bin1, Bin10, Bin11, Bin12, Bin13, Bin14, Bin15, Bin16, Bin17, Bin18, Bin19, Bin2, Bin20, Bin21, Bin22, Bin23, Bin24, Bin25, Bin26, Bin27, Bin28, Bin29,
    Bin3, Bin30, Bin31, Bin32, Bin33, Bin34, Bin35, Bin36, Bin37, Bin38, Bin39, Bin4, Bin40, Bin41, Bin42, Bin43, Bin44, Bin45, Bin46, Bin47, Bin48, Bin49, Bin5,
    Bin50, Bin51, Bin52, Bin53, Bin54, Bin55, Bin56, Bin57, Bin58, Bin59, Bin6, Bin60, Bin61, Bin62, Bin63, Bin64, Bin65, Bin66, Bin67, Bin68, Bin69, Bin7, Bin70,
    Bin71, Bin72, Bin8, Bin9, Count, Date, DeployID, DepthSensor, DOY, HistType, Hour, Instr, Latitude, LocationQuality, Longitude, NumBins, Ptt, Source, Sum,
    Time.Offset

The following objects are masked from histos (pos = 32):

    BadTherm, Bin1, Bin10, Bin11, Bin12, Bin13, Bin14, Bin15, Bin16, Bin17, Bin18, Bin19, Bin2, Bin20, Bin21, Bin22, Bin23, Bin24, Bin25, Bin26, Bin27, Bin28, Bin29,
    Bin3, Bin30, Bin31, Bin32, Bin33, Bin34, Bin35, Bin36, Bin37, Bin38, Bin39, Bin4, Bin40, Bin41, Bin42, Bin43, Bin44, Bin45, Bin46, Bin47, Bin48, Bin49, Bin5,
    Bin50, Bin51, Bin52, Bin53, Bin54, Bin55, Bin56, Bin57, Bin58, Bin59, Bin6, Bin60, Bin61, Bin62, Bin63, Bin64, Bin65, Bin66, Bin67, Bin68, Bin69, Bin7, Bin70,
    Bin71, Bin72, Bin8, Bin9, Count, Date, DeployID, DepthSensor, DOY, HistType, Hour, Instr, Latitude, LocationQuality, Longitude, NumBins, Ptt, Source, Sum,
    Time.Offset

The following objects are masked from histos (pos = 33):

    BadTherm, Bin1, Bin10, Bin11, Bin12, Bin13, Bin14, Bin15, Bin16, Bin17, Bin18, Bin19, Bin2, Bin20, Bin21, Bin22, Bin23, Bin24, Bin25, Bin26, Bin27, Bin28, Bin29,
    Bin3, Bin30, Bin31, Bin32, Bin33, Bin34, Bin35, Bin36, Bin37, Bin38, Bin39, Bin4, Bin40, Bin41, Bin42, Bin43, Bin44, Bin45, Bin46, Bin47, Bin48, Bin49, Bin5,
    Bin50, Bin51, Bin52, Bin53, Bin54, Bin55, Bin56, Bin57, Bin58, Bin59, Bin6, Bin60, Bin61, Bin62, Bin63, Bin64, Bin65, Bin66, Bin67, Bin68, Bin69, Bin7, Bin70,
    Bin71, Bin72, Bin8, Bin9, Count, Date, DeployID, DepthSensor, HistType, Instr, Latitude, LocationQuality, Longitude, NumBins, Ptt, Source, Sum, Time.Offset

head(histos)

# Load libraries
library(mgcv)
library(gamm4)
library(pspline)
library(suncalc)
library(lmerTest)
library(MuMIn)
library(parallel)
library(mgcViz)
library(rgl)
library(feather)
library(data.table)
library(car)
library(ggplot2)
library(grid)
library(animation)
library(itsadug)
library(visreg)
library(dplyr)
library(viridis)
library(hrbrthemes)

TAD

Reorganize data

TAD.1<- c(Bin1[which(HistType == "TAD")], Bin2[which(HistType == "TAD")], Bin3[which(HistType == "TAD")], Bin4[which(HistType == "TAD")], Bin5[which(HistType == "TAD")],
                 Bin6[which(HistType == "TAD")], Bin7[which(HistType == "TAD")], Bin8[which(HistType == "TAD")], Bin9[which(HistType == "TAD")], Bin10[which(HistType == "TAD")],
                 Bin11[which(HistType == "TAD")], Bin12[which(HistType == "TAD")], Bin13[which(HistType == "TAD")])
TAD.2<-TAD.1/100
n<-length(Bin1[which(HistType == "TAD")])
Bin<- c(replicate(n, Bin1[which(HistType == "TADLIMITS")]), replicate(n, Bin2[which(HistType == "TADLIMITS")]), replicate(n, Bin3[which(HistType == "TADLIMITS")]),
                  replicate(n, Bin4[which(HistType == "TADLIMITS")]), replicate(n, Bin5[which(HistType == "TADLIMITS")]), replicate(n, Bin6[which(HistType == "TADLIMITS")]),
                  replicate(n, Bin7[which(HistType == "TADLIMITS")]), replicate(n, Bin8[which(HistType == "TADLIMITS")]), replicate(n, Bin9[which(HistType == "TADLIMITS")]),
                  replicate(n, Bin10[which(HistType == "TADLIMITS")]), replicate(n, Bin11[which(HistType == "TADLIMITS")]), replicate(n, Bin12[which(HistType == "TADLIMITS")]),
                  replicate(n, Bin13[which(HistType == "TADLIMITS")]))
Day<- c(replicate(13, DOY[which(HistType == "TAD")]))
Hour.1<- c(replicate(13, Hour[which(HistType == "TAD")]))

df.TAD<- data.frame(TAD.2, Bin, Day, Hour.1)
head(df.TAD)

Plotting TAD across time by bin

x<- df.TAD$Day[which(Hour.1 == 0)]
y<- df.TAD$TAD.2[which(Hour.1 == 0)]
f<- df.TAD$Bin[which(Hour.1 == 0)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_area(alpha=1) +
  scale_fill_manual(values = c("purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Proportion") + ggtitle("TAD Bins, Hour 0") + labs(fill = "Bin")

x<- df.TAD$Day[which(Hour.1 == 6)]
y<- df.TAD$TAD.2[which(Hour.1 == 6)]
f<- df.TAD$Bin[which(Hour.1 == 6)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_area(alpha=1) +
  scale_fill_manual(values = c("purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Proportion") + ggtitle("TAD Bins, Hour 6") + labs(fill = "Bin")

x<- df.TAD$Day[which(Hour.1 == 12)]
y<- df.TAD$TAD.2[which(Hour.1 == 12)]
f<- df.TAD$Bin[which(Hour.1 == 12)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_area(alpha=1) +
  scale_fill_manual(values = c("purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Proportion") + ggtitle("TAD Bins, Hour 12") + labs(fill = "Bin")

x<- df.TAD$Day[which(Hour.1 == 18)]
y<- df.TAD$TAD.2[which(Hour.1 == 18)]
f<- df.TAD$Bin[which(Hour.1 == 18)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_area(alpha=1) +
  scale_fill_manual(values = c("purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Proportion") + ggtitle("TAD Bins, Hour 18") + labs(fill = "Bin")

Number of Dives by Hour Block

plot(Sum[which(HistType == "DiveDepth" & Hour == "0")] ~ DOY[which(HistType == "DiveDepth" & Hour == "0")], type = "b", pch = 16, col = "navy", ylim = c(0,300),
     xlab = "Day of Year", ylab = "Number of Dives")
lines(Sum[which(HistType == "DiveDepth" & Hour == "6")] ~ DOY[which(HistType == "DiveDepth" & Hour == "6")], type = "b", pch = 16, col = "blue")

lines(Sum[which(HistType == "DiveDepth" & Hour == "12")] ~ DOY[which(HistType == "DiveDepth" & Hour == "12")], type = "b", pch = 16, col = "turquoise3")
lines(Sum[which(HistType == "DiveDepth" & Hour == "18")] ~ DOY[which(HistType == "DiveDepth" & Hour == "18")], type = "b", pch = 16, col = "green")

legend(50, 300, legend=c("Hour 0","Hour 6", "Hour 12","Hour 18"),
       col=c("navy","blue","turquoise3","green"), lty=1, cex=0.8)

Number of Dives - Summed By Day

agg1<- aggregate(Sum[which(HistType == "DiveDepth")], list(DOY[which(HistType == "DiveDepth")]), FUN = sum)
plot(agg1$x ~ agg1$Group.1, type = "b", pch = 16, col = "blue", ylim = c(0,900),
     xlab = "Day of Year", ylab = "Number of Dives")

Depth

Reorganize data

Dep<- c(Bin1[which(HistType == "DiveDepth")], Bin2[which(HistType == "DiveDepth")], Bin3[which(HistType == "DiveDepth")], Bin4[which(HistType == "DiveDepth")], Bin5[which(HistType == "DiveDepth")], Bin6[which(HistType == "DiveDepth")], Bin7[which(HistType == "DiveDepth")], Bin8[which(HistType == "DiveDepth")], Bin9[which(HistType == "DiveDepth")], Bin10[which(HistType == "DiveDepth")],  Bin11[which(HistType == "DiveDepth")], Bin12[which(HistType == "DiveDepth")], Bin13[which(HistType == "DiveDepth")], Bin14[which(HistType == "DiveDepth")])
n<-length(Bin1[which(HistType == "DiveDepth")])
Bin<- c(replicate(n, Bin1[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin2[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin3[which(HistType == "DiveDepthLIMITS")]),
                  replicate(n, Bin4[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin5[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin6[which(HistType == "DiveDepthLIMITS")]),
                  replicate(n, Bin7[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin8[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin9[which(HistType == "DiveDepthLIMITS")]),
                  replicate(n, Bin10[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin11[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin12[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin13[which(HistType == "DiveDepthLIMITS")]), replicate(n, Bin14[which(HistType == "DiveDepthLIMITS")]))
Day<- c(replicate(14, DOY[which(HistType == "DiveDepth")]))
Hour.1<- c(replicate(14, Hour[which(HistType == "DiveDepth")]))

df.Dep<- data.frame(Dep, Bin, Day, Hour.1)
head(df.Dep)

Plotting Depth across time by bin

x<- df.Dep$Day[which(Hour.1 == 0)]
y<- df.Dep$Dep[which(Hour.1 == 0)]
f<- df.Dep$Bin[which(Hour.1 == 0)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_bar(position = "stack", stat = "identity") +
  scale_fill_manual(values = c("purple4","purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Count") + ggtitle("Depth Bins, Hour 0") + labs(fill = "Bin")

x<- df.Dep$Day[which(Hour.1 == 6)]
y<- df.Dep$Dep[which(Hour.1 == 6)]
f<- df.Dep$Bin[which(Hour.1 == 6)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_bar(position = "stack", stat = "identity") +
  scale_fill_manual(values = c("purple4","purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Count") + ggtitle("Depth Bins, Hour 6") + labs(fill = "Bin")

x<- df.Dep$Day[which(Hour.1 == 12)]
y<- df.Dep$Dep[which(Hour.1 == 12)]
f<- df.Dep$Bin[which(Hour.1 == 12)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_bar(position = "stack", stat = "identity") +
  scale_fill_manual(values = c("purple4","purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Count") + ggtitle("Depth Bins, Hour 12") + labs(fill = "Bin")

x<- df.Dep$Day[which(Hour.1 == 18)]
y<- df.Dep$Dep[which(Hour.1 == 18)]
f<- df.Dep$Bin[which(Hour.1 == 18)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_bar(position = "stack", stat = "identity") +
  scale_fill_manual(values = c("purple4","purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Count") + ggtitle("Depth Bins, Hour 18") + labs(fill = "Bin")

Duration

Reorganize data

Dur<- c(Bin1[which(HistType == "DiveDuration")], Bin2[which(HistType == "DiveDuration")], Bin3[which(HistType == "DiveDuration")], Bin4[which(HistType == "DiveDuration")], Bin5[which(HistType == "DiveDuration")], Bin6[which(HistType == "DiveDuration")], Bin7[which(HistType == "DiveDuration")], Bin8[which(HistType == "DiveDuration")], Bin9[which(HistType == "DiveDuration")], Bin10[which(HistType == "DiveDuration")],  Bin11[which(HistType == "DiveDuration")], Bin12[which(HistType == "DiveDuration")], Bin13[which(HistType == "DiveDuration")], Bin14[which(HistType == "DiveDuration")])
n<-length(Bin1[which(HistType == "DiveDuration")])
Bin<- c(replicate(n, Bin1[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin2[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin3[which(HistType == "DiveDurationLIMITS")]),
                  replicate(n, Bin4[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin5[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin6[which(HistType == "DiveDurationLIMITS")]),replicate(n, Bin7[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin8[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin9[which(HistType == "DiveDurationLIMITS")]),  replicate(n, Bin10[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin11[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin12[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin13[which(HistType == "DiveDurationLIMITS")]), replicate(n, Bin14[which(HistType == "DiveDurationLIMITS")]))
Day<- c(replicate(14, DOY[which(HistType == "DiveDuration")]))
Hour.1<- c(replicate(14, Hour[which(HistType == "DiveDuration")]))

df.Dur<- data.frame(Dur, Bin, Day, Hour.1)
head(df.Dur)

Plotting Duration across time by bin

x<- df.Dur$Day[which(Hour.1 == 0)]
y<- df.Dur$Dur[which(Hour.1 == 0)]
f<- df.Dur$Bin[which(Hour.1 == 0)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_bar(position = "stack", stat = "identity") +
  scale_fill_manual(values = c("purple4","purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Count") + ggtitle("Duration Bins, Hour 0") + labs(fill = "Bin")

x<- df.Dur$Day[which(Hour.1 == 6)]
y<- df.Dur$Dur[which(Hour.1 == 6)]
f<- df.Dur$Bin[which(Hour.1 == 6)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_bar(position = "stack", stat = "identity") +
  scale_fill_manual(values = c("purple4","purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Count") + ggtitle("Duration Bins, Hour 6") + labs(fill = "Bin")

x<- df.Dur$Day[which(Hour.1 == 12)]
y<- df.Dur$Dur[which(Hour.1 == 12)]
f<- df.Dur$Bin[which(Hour.1 == 12)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_bar(position = "stack", stat = "identity") +
  scale_fill_manual(values = c("purple4","purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Count") + ggtitle("Duration Bins, Hour 12") + labs(fill = "Bin")

x<- df.Dur$Day[which(Hour.1 == 18)]
y<- df.Dur$Dur[which(Hour.1 == 18)]
f<- df.Dur$Bin[which(Hour.1 == 18)]
df1<- data.frame(x,y,f)
ggplot(df1, aes(x=x, y=y, fill=as.factor(f))) + 
  geom_bar(position = "stack", stat = "identity") +
  scale_fill_manual(values = c("purple4","purple","navy","blue","turquoise3","cyan","darkgreen","green","darkgoldenrod","yellow","sienna1","red","red4","black")) + theme_linedraw()+
 theme(axis.text=element_text(size=12)) + xlab("Day of Year") + ylab("Count") + ggtitle("Duration Bins, Hour 18") + labs(fill = "Bin")

Surface Time

Reorganize data

Surf<- c(Bin1[which(HistType == "Percent")], Bin2[which(HistType == "Percent")], Bin3[which(HistType == "Percent")], Bin4[which(HistType == "Percent")], Bin5[which(HistType == "Percent")], Bin6[which(HistType == "Percent")], Bin7[which(HistType == "Percent")], Bin8[which(HistType == "Percent")], Bin9[which(HistType == "Percent")], Bin10[which(HistType == "Percent")],  Bin11[which(HistType == "Percent")], Bin12[which(HistType == "Percent")], Bin13[which(HistType == "Percent")], Bin14[which(HistType == "Percent")], Bin15[which(HistType == "Percent")], Bin16[which(HistType == "Percent")], Bin17[which(HistType == "Percent")], Bin18[which(HistType == "Percent")], Bin19[which(HistType == "Percent")], Bin20[which(HistType == "Percent")], Bin21[which(HistType == "Percent")], Bin22[which(HistType == "Percent")], Bin23[which(HistType == "Percent")], Bin24[which(HistType == "Percent")])

n<-length(Bin1[which(HistType == "Percent")])
Hour.1<- c(replicate(n, 0),replicate(n, 1),replicate(n, 2),replicate(n, 3),replicate(n, 4),replicate(n, 5),replicate(n, 6),replicate(n, 7),replicate(n, 8),replicate(n, 9),replicate(n, 10),replicate(n, 11),replicate(n, 12),replicate(n, 13),replicate(n, 14),replicate(n, 15),replicate(n, 16),replicate(n, 17),replicate(n, 18),replicate(n, 19),replicate(n, 20),replicate(n, 21),replicate(n, 22),replicate(n, 23))
Day<- c(replicate(24, DOY[which(HistType == "Percent")]))

df.Surf<- data.frame(Surf, Hour.1, Day)
head(df.Surf)

Plotting %Surface by Hour and Day

p <-ggplot(df.Surf,aes(Day,Hour.1,fill=Surf))+
  geom_tile(color= "grey",size=0.1) + 
  scale_fill_viridis(name="Surface (%)",option ="turbo", limits = c(0, 100), oob = scales::squish)
p

Look at the ‘Behavior’ File

behav<- read.csv("238272-Behavior.csv")
attach(behav)

The following objects are masked from behav (pos = 3):

    Count, Deep, DeployID, DepthMax, DepthMin, DepthSensor, DurationMax, DurationMin, End, Instr, Number, Ptt, Shallow, Shape, Source, Start, What

The following objects are masked from histos (pos = 4):

    Count, DeployID, DepthSensor, DOY, Instr, Ptt, Source

The following objects are masked from histos (pos = 5):

    Count, DeployID, DepthSensor, DOY, Instr, Ptt, Source

The following objects are masked from histos (pos = 6):

    Count, DeployID, DepthSensor, DOY, Instr, Ptt, Source

The following objects are masked from histos (pos = 34):

    Count, DeployID, DepthSensor, DOY, Instr, Ptt, Source

The following objects are masked from histos (pos = 35):

    Count, DeployID, DepthSensor, Instr, Ptt, Source

head(behav)

Depth over deployment

boxplot(DepthMin[which(What == "Dive")] ~ DOY[which(What == "Dive")], main = "Min Depth", ylab = "Depth (m)", xlab = "Day of Year")

boxplot(DepthMax[which(What == "Dive")] ~ DOY[which(What == "Dive")], main = "Max Depth", ylab = "Depth (m)", xlab = "Day of Year")

Duration over deployment

boxplot(DurationMin[which(What == "Dive")] ~ DOY[which(What == "Dive")], main = "Min Duration", ylab = "Duration (min)", xlab = "Day of Year")

boxplot(DurationMax[which(What == "Dive")] ~ DOY[which(What == "Dive")], main = "Max Duration", ylab = "Duration (min)", xlab = "Day of Year")