paelo<-read.table(file.choose(),sep=",",header=TRUE)
paelo
attach(paelo)

#this is for all dives
summary(paelo)

h_paelo_length<-hist(Length_cm, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="Paelopatides confundens (cm)")

###############################################################################

#December 2006 length data 
#must make a vector that is only the length measurements from that time period
paeloDec2006<-Length_cm[DiveNumber=="1067"]
paeloDec2006

#Because there are no observations in this time period for Paelopatides confundens, use
#Nov 2012 graph below with June 2007 labels so I can erase the line in Powerpoint

h_paelo_length<-hist(Length_cm, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="December 2006 Paelopatides confundens (cm)")


summary(paeloDec2006)
#boxplot(paeloDec2006)

#############################################################################

#February 2007 length data 
#must make a vector that is only the length measurements from that time period
paelo1077<-Length_cm[DiveNumber=="1077"]
paelo1077

paelo1080<-Length_cm[DiveNumber=="1080"]
paelo1080

paeloFeb2007<-c(paelo1077,paelo1080)
paeloFeb2007

h_paeloFeb2007<-hist(paeloFeb2007, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="February 2007 Paelopatides confundens (cm)")

#Because there are no observations in this time period for Paelopatides confundens, use
#Nov 2012 graph below with June 2007 labels so I can erase the line in Powerpoint

h_paelo_length<-hist(Length_cm, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="February 2007 Paelopatides confundens (cm)")

summary(paeloFeb2007)
#boxplot(paeloFeb2007)

###########################################################################

#June 2007 length data 
#must make a vector that is only the length measurements from that time period
paeloJune2007<-Length_cm[DiveNumber=="1094"]
paeloJune2007

h_paeloJune2007<-hist(paeloJune2007, breaks=12,xlim=c(0,30),ylim=c(0,10),
                      main= "  ", border= "dark blue", col = "light blue",  
                      xlab="June 2007 Paelopatides confundens (cm)")

#Because there are no observations in this time period for Paelopatides confundens, use
#Nov 2012 graph below with June 2007 labels so I can erase the line in Powerpoint

h_paelo_length<-hist(Length_cm, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="June 2007 Paelopatides confundens (cm)")

summary(paeloJune2007)
#boxplot(paeloJune2007)

############################################################################

#Sept 2007 length data 
#must make a vector that is only the length measurements from that time period
paelo1141<-Length_cm[DiveNumber=="1141"]
paelo1141

paelo1143<-Length_cm[DiveNumber=="1143"]
paelo1143

paeloSept2007<-c(paelo1141,paelo1143)
paeloSept2007

h_paeloSept2007<-hist(paeloSept2007, breaks=12,xlim=c(0,30),ylim=c(0,10),
                      main= "  ", border= "dark blue", col = "light blue",  
                      xlab="September 2007 Paelopatides confundens (cm)")

#Because there are no observations in this time period for Paelopatides confundens, use
#Nov 2012 graph below with June 2007 labels so I can erase the line in Powerpoint

h_paelo_length<-hist(Length_cm, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="September 2007 Paelopatides confundens (cm)")


summary(paeloSept2007)
#boxplot(paeloSept2007)

##############################################################################

#February 2009 length data for paeloidia
#must make a vector that is only the length measurements from that time period
paeloFeb2009<-Length_cm[DiveNumber=="8"]
paeloFeb2009

h_paeloFeb2009<-hist(paeloFeb2009, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="February 2009 Paelopatides confundens (cm)")

#Because there are no observations in this time period for Paelopatides confundens, use
#Nov 2012 graph below with June 2007 labels so I can erase the line in Powerpoint

h_paelo_length<-hist(Length_cm, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="February 2009 Paelopatides confundens (cm)")

summary(paeloFeb2009)
boxplot(paeloFeb2009)

###############################################################################

#May 2011 length data 
#must make a vector that is only the length measurements from that time period
paelo230<-Length_cm[DiveNumber=="230"]
paelo230

paelo231<-Length_cm[DiveNumber=="231"]
paelo231

paelo232<-Length_cm[DiveNumber=="232"]
paelo232

paeloMay2011<-c(paelo230,paelo231,paelo232)
paeloMay2011

h_paeloMay2011<-hist(paeloMay2011, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="May 2011 Paelopatides confundens (cm)")

#Because there are no observations in this time period for Paelopatides confundens, use
#Nov 2012 graph below with June 2007 labels so I can erase the line in Powerpoint

h_paelo_length<-hist(Length_cm, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="May 2011 Paelopatides confundens (cm)")

summary(paeloMay2011)
#boxplot(paeloMay2011)

################################################################################


#Nov 2011 length data 
#must make a vector that is only the length measurements from that time period
paelo321<-Length_cm[DiveNumber=="321"]
paelo321

paelo323<-Length_cm[DiveNumber=="323"]
paelo323

paelo324<-Length_cm[DiveNumber=="324"]
paelo324

paeloNov2011<-c(paelo321,paelo323,paelo324)
paeloNov2011

h_paeloNov2011<-hist(paeloNov2011, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="November 2011 Paelopatides confundens (cm)")

#this is done to fix the axis
#paelofix<-c(30)
#paeloNov2011a<-c(paeloNov2011,paelofix)

#h_paeloNov2011<-hist(paeloNov2011a, breaks=12,xlim=c(0,30),ylim=c(0,10),
main= "  ", border= "dark blue", col = "light blue",  
xlab="November 2011 Paelopatides confundens (cm)")

summary(paeloNov2011)
#boxplot(paeloNov2011)

#################################################################################

#June 2012 length data for paeloidia
#must make a vector that is only the length measurements from that time period
paeloJune2012<-Length_cm[DiveNumber=="403"]
paeloJune2012

h_paeloJune2012<-hist(paeloJune2012, breaks=12,xlim=c(0,30),ylim=c(0,10),
                      main= "  ", border= "dark blue", col = "light blue",  
                      xlab="June 2012 Paelopatides confundens (cm)")

#this is done to fix the axis
#paelofix<-c(30)
#paeloJune2012a<-c(paeloJune2012,paelofix)

#h_paeloJune2012<-hist(paeloJune2012a, breaks=12,xlim=c(0,30),ylim=c(0,10),
main= "  ", border= "dark blue", col = "light blue",  
xlab="June 2012 Paelopatides confundens (cm)")

summary(paeloJune2012)
boxplot(paeloJune2012)

#############################################################################

#November 2012 length data 
#must make a vector that is only the length measurements from that time period
paelo442<-Length_cm[DiveNumber=="442"]
paelo442

paelo443<-Length_cm[DiveNumber=="443"]
paelo443

paeloNov2012<-c(paelo442,paelo443)
paeloNov2012


h_paeloNov2012<-hist(paeloNov2012, breaks=12,xlim=c(0,30),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="November 2012 Paelopatides confundens (cm)")

#this is done to fix the axis
#paelofix1<-c(0)
#paeloNov2012a<-c(paeloNov2012,paelofix1)
#paeloNov2012a

#h_paeloNov2012<-hist(paeloNov2012a, breaks=12,xlim=c(0,30),ylim=c(0,10),
main= "  ", border= "dark blue", col = "light blue",  
xlab="November 2012 Paelopatides confundens (cm)")

summary(paeloNov2012)
#boxplot(paeloNov2012)