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

#this is for all dives
summary(pen_sp1)

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="Peniagone sp. 1 Length (cm)")

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

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

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="December 2006 Peniagone sp. 1 (cm)")


summary(pen_sp1Dec2006)
#boxplot(pen_sp1Dec2006)

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

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

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

pen_sp1Feb2007<-c(pen_sp11077,pen_sp11080)
pen_sp1Feb2007

h_pen_sp1Feb2007<-hist(pen_sp1Feb2007, breaks=12,xlim=c(0,8),ylim=c(0,60),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2007 Peniagone sp. 1 (cm)")

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2007 Peniagone sp. 1 (cm)")

summary(pen_sp1Feb2007)
#boxplot(pen_sp1Feb2007)

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

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

h_pen_sp1June2007<-hist(pen_sp1June2007, breaks=12,xlim=c(0,8),ylim=c(0,60),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2007 Peniagone sp. 1 (cm)")

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="June 2007 Peniagone sp. 1 (cm)")

summary(pen_sp1June2007)
#boxplot(pen_sp1June2007)

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

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

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

pen_sp1Sept2007<-c(pen_sp11141,pen_sp11143)
pen_sp1Sept2007

h_pen_sp1Sept2007<-hist(pen_sp1Sept2007, breaks=12,xlim=c(0,8),ylim=c(0,60),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="September 2007 Peniagone sp. 1 (cm)")

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="September 2007 Peniagone sp. 1 (cm)")


summary(pen_sp1Sept2007)
#boxplot(pen_sp1Sept2007)

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

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

h_pen_sp1Feb2009<-hist(pen_sp1Feb2009, breaks=12,xlim=c(0,8),ylim=c(0,60),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2009 Peniagone sp. 1 (cm)")

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2009 Peniagone sp. 1 (cm)")

summary(pen_sp1Feb2009)
boxplot(pen_sp1Feb2009)

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

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

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

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

pen_sp1May2011<-c(pen_sp1230,pen_sp1231,pen_sp1232)
pen_sp1May2011

h_pen_sp1May2011<-hist(pen_sp1May2011, breaks=12,xlim=c(0,8),ylim=c(0,60),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="May 2011 Peniagone sp. 1 (cm)")

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="May 2011 Peniagone sp. 1 (cm)")

summary(pen_sp1May2011)
#boxplot(pen_sp1May2011)

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


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

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

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

pen_sp1Nov2011<-c(pen_sp1321,pen_sp1323,pen_sp1324)
pen_sp1Nov2011

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="November 2011 Peniagone sp. 1 (cm)")

summary(pen_sp1Nov2011)
#boxplot(pen_sp1Nov2011)

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

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

h_pen_sp1June2012<-hist(pen_sp1June2012, breaks=12,xlim=c(0,8),ylim=c(0,60),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2012 Peniagone sp. 1 (cm)")

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

h_pen_sp1_length<-hist(Length_cm, breaks=12,xlim=c(0,8),ylim=c(0,60),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="June 2012 Peniagone sp. 1 (cm)")

summary(pen_sp1June2012)
boxplot(pen_sp1June2012)

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

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

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

pen_sp1Nov2012<-c(pen_sp1442,pen_sp1443)
pen_sp1Nov2012


h_pen_sp1Nov2012<-hist(pen_sp1Nov2012, breaks=12,xlim=c(0,8),ylim=c(0,60),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2012 Peniagone sp. 1 (cm)")



summary(pen_sp1Nov2012)
#boxplot(pen_sp1Nov2012)