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

#this is for all dives
summary(pen_sp2)

#only histogram for length done because this is grouped with the
#Peniagone complex for biomass purposes

h_pen_sp2_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="Peniagone sp. 2 Length (cm)")

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

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

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

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

h_pen_sp2_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="December 2006 Peniagone sp. 2 (cm)")


summary(pen_sp2Dec2006)
#boxplot(pen_sp2Dec2006)

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

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

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

pen_sp2Feb2007<-c(pen_sp21077,pen_sp21080)
pen_sp2Feb2007

h_pen_sp2Feb2007<-hist(pen_sp2Feb2007, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2007 Peniagone sp. 2 (cm)")

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

h_pen_sp2_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2007 Peniagone sp. 2 (cm)")

summary(pen_sp2Feb2007)
#boxplot(pen_sp2Feb2007)

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

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

h_pen_sp2June2007<-hist(pen_sp2June2007, breaks=12,xlim=c(0,15),ylim=c(0,20),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="June 2007 Peniagone sp. 2 (cm)")

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

h_pen_sp2_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="June 2007 Peniagone sp. 2 (cm)")

summary(pen_sp2June2007)
#boxplot(pen_sp2June2007)

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

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

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

pen_sp2Sept2007<-c(pen_sp21141,pen_sp21143)
pen_sp2Sept2007

h_pen_sp2Sept2007<-hist(pen_sp2Sept2007, breaks=12,xlim=c(0,15),ylim=c(0,20),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="September 2007 Peniagone sp. 2 (cm)")

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

h_pen_sp2_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="September 2007 Peniagone sp. 2 (cm)")


summary(pen_sp2Sept2007)
#boxplot(pen_sp2Sept2007)

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

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

h_pen_sp2Feb2009<-hist(pen_sp2Feb2009, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2009 Peniagone sp. 2 (cm)")

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

h_pen_sp2_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2009 Peniagone sp. 2 (cm)")

summary(pen_sp2Feb2009)
boxplot(pen_sp2Feb2009)

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

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

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

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

pen_sp2May2011<-c(pen_sp2230,pen_sp2231,pen_sp2232)
pen_sp2May2011

h_pen_sp2May2011<-hist(pen_sp2May2011, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="May 2011 Peniagone sp. 2 (cm)")


summary(pen_sp2May2011)
#boxplot(pen_sp2May2011)

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


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

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

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

pen_sp2Nov2011<-c(pen_sp2321,pen_sp2323,pen_sp2324)
pen_sp2Nov2011

h_pen_sp2Nov2011<-hist(pen_sp2Nov2011, breaks=12,xlim=c(0,15),ylim=c(0,20),
                   main= "  ", border= "dark blue", col = "light blue",
                   xlab="November 2011 Peniagone sp. 2 (cm)")

#this is to better fit the axis
#pen_sp2fix<-c(14)

#pen_sp2Nov2011a<-c(pen_sp2Nov2011,pen_sp2fix)

#h_pen_sp2Nov2011<-hist(pen_sp2Nov2011a, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",
                       xlab="November 2011 Peniagone sp. 2 (cm)")

summary(pen_sp2Nov2011)
#boxplot(pen_sp2Nov2011)

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

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

h_pen_sp2June2012<-hist(pen_sp2June2012, breaks=12,xlim=c(0,15),ylim=c(0,20),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="June 2012 Peniagone sp. 2 (cm)")

#this is to better fit the axis
#pen_sp2fix<-c(14)

#pen_sp2June2012a<-c(pen_sp2June2012,pen_sp2fix)
#pen_sp2June2012a

#h_pen_sp2June2012<-hist(pen_sp2June2012a, breaks=12,xlim=c(0,15),ylim=c(0,20),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="June 2012 Peniagone sp. 2 (cm)")

summary(pen_sp2June2012)
boxplot(pen_sp2June2012)

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

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

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

pen_sp2Nov2012<-c(pen_sp2442,pen_sp2443)
pen_sp2Nov2012


h_pen_sp2Nov2012<-hist(pen_sp2Nov2012, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="November 2012 Peniagone sp. 2 (cm)")

#this is to better fit the axis
#pen_sp2fix<-c(14)

#pen_sp2Nov2012a<-c(pen_sp2Nov2012,pen_sp2fix)
#pen_sp2Nov2012a

h_pen_sp2Nov2012<-hist(pen_sp2Nov2012a, breaks=12,xlim=c(0,15),ylim=c(0,20),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="November 2012 Peniagone sp. 2 (cm)")


summary(pen_sp2Nov2012)
#boxplot(pen_sp2Nov2012)