pen_papillata<-read.table(file.choose(),sep=",",header=TRUE)
pen_papillata
attach(pen_papillata)

#this is for all dives
summary(pen_papillata)

#only histogram for length done because this is grouped with the
#Peniagone complex for biomass purposes

h_pen_papillata_length<-hist(Length_cm, breaks=12,xlim=c(0,25),ylim=c(0,70),,
                             main= "  ", border= "dark blue", col = "light blue",  
                             xlab="Peniagone papillata Length (cm)")

########################################################################

#December 2006 length data 
#must make a vector that is only the length measurements from that time period
pen_papillataDec2006<-Length_cm[DiveNumber=="1067"]
pen_papillataDec2006

h_pen_papillataDec2006<-hist(pen_papillataDec2006, breaks=12,xlim=c(0,25),ylim=c(0,70),,
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="December 2006 Peniagone papillata (cm)")
#This was done to help with the axis
#pen_papillatafix<-c(24)
# 
#pen_papillataDec2006<-c(pen_papillataDec2006,pen_papillatafix)
#pen_papillataDec2006

#h_pen_papillataDec2006<-hist(pen_papillataDec2006, breaks=12,xlim=c(0,25),ylim=c(0,70),,
                             main= "  ", border= "dark blue", col = "light blue",  
                             xlab="December 2006 Peniagone papillata (cm)")
summary(pen_papillataDec2006)
#boxplot(pen_papillataDec2006)

#############################################################################

#February 2007 length data 
#must make a vector that is only the length measurements from that time period
pen_papillata1077<-Length_cm[DiveNumber=="1077"]
pen_papillata1077

pen_papillata1080<-Length_cm[DiveNumber=="1080"]
pen_papillata1080

pen_papillataFeb2007<-c(pen_papillata1077,pen_papillata1080)
pen_papillataFeb2007

h_pen_papillataFeb2007<-hist(pen_papillataFeb2007, breaks=12,xlim=c(0,25),ylim=c(0,70),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2007 Peniagone papillata (cm)")

#This was done to help with the axis
#pen_papillatafix<-c(24)
# 
#pen_papillataFeb2007<-c(pen_papillata1077,pen_papillata1080,pen_papillatafix)
#pen_papillataFeb2007

#h_pen_papillataFeb2007<-hist(pen_papillataFeb2007, breaks=12,xlim=c(0,25),ylim=c(0,70),
#main= "  ", border= "dark blue", col = "light blue",  
#xlab="February 2007 Peniagone papillata (cm)")

summary(pen_papillataFeb2007)
#boxplot(pen_papillataFeb2007)

###########################################################################

#June 2007 length data 
#must make a vector that is only the length measurements from that time period
pen_papillataJune2007<-Length_cm[DiveNumber=="1094"]
pen_papillataJune2007

h_pen_papillataJune2007<-hist(pen_papillataJune2007, breaks=12,xlim=c(0,25),ylim=c(0,70),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2007 Peniagone papillata (cm)")

#This was done to help with the axis
pen_papillatafix<-c(24)
# 
#pen_papillataJune2007a<-c(pen_papillataJune2007,pen_papillatafix)
#pen_papillataJune2007a

#h_pen_papillataFeb2007<-hist(pen_papillataJune2007a, breaks=12,xlim=c(0,25),ylim=c(0,70),
main= "  ", border= "dark blue", col = "light blue",  
xlab="June 2007 Peniagone papillata (cm)")

summary(pen_papillataJune2007)
#boxplot(pen_papillataJune2007)

############################################################################

#Sept 2007 length data 
#must make a vector that is only the length measurements from that time period
pen_papillata1141<-Length_cm[DiveNumber=="1141"]
pen_papillata1141

pen_papillata1143<-Length_cm[DiveNumber=="1143"]
pen_papillata1143

pen_papillataSept2007<-c(pen_papillata1141,pen_papillata1143)
pen_papillataSept2007

h_pen_papillataSept2007<-hist(pen_papillataSept2007, breaks=12,xlim=c(0,25),ylim=c(0,70),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="September 2007 Peniagone papillata (cm)")

#This was done to help with the axis
#pen_papillatafix<-c(24)
# 
#pen_papillataSept2007a<-c(pen_papillataSept2007,pen_papillatafix)
#pen_papillataSept2007a

#h_pen_papillataSept2007<-hist(pen_papillataSept2007a, breaks=12,xlim=c(0,25),ylim=c(0,70),
main= "  ", border= "dark blue", col = "light blue",  
xlab="September 2007 Peniagone papillata (cm)")


summary(pen_papillataSept2007)
#boxplot(pen_papillataSept2007)

##############################################################################

#February 2009 length data for pen_papillataidia
#must make a vector that is only the length measurements from that time period
pen_papillataFeb2009<-Length_cm[DiveNumber=="8"]
pen_papillataFeb2009

h_pen_papillataFeb2009<-hist(pen_papillataFeb2009, breaks=12,xlim=c(0,25),ylim=c(0,70),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2009 Peniagone papillata (cm)")

#This was done to help with the axis
#pen_papillatafix<-c(24)
# 
#pen_papillataFeb2009a<-c(pen_papillataFeb2009,pen_papillatafix)
#pen_papillataFeb2009a

#h_pen_papillataFeb2009<-hist(pen_papillataFeb2009a, breaks=12,xlim=c(0,25),ylim=c(0,70),
main= "  ", border= "dark blue", col = "light blue",  
xlab="February 2009 Peniagone papillata (cm)")

summary(pen_papillataFeb2009)
boxplot(pen_papillataFeb2009)

###############################################################################

#May 2011 length data 
#must make a vector that is only the length measurements from that time period
pen_papillata230<-Length_cm[DiveNumber=="230"]
pen_papillata230

pen_papillata231<-Length_cm[DiveNumber=="231"]
pen_papillata231

pen_papillata232<-Length_cm[DiveNumber=="232"]
pen_papillata232

pen_papillataMay2011<-c(pen_papillata230,pen_papillata231,pen_papillata232)
pen_papillataMay2011

h_pen_papillataMay2011<-hist(pen_papillataMay2011, breaks=12,xlim=c(0,25),ylim=c(0,70),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="May 2011 Peniagone papillata (cm)")

summary(pen_papillataMay2011)
#boxplot(pen_papillataMay2011)

################################################################################


#Nov 2011 length data 
#must make a vector that is only the length measurements from that time period
pen_papillata321<-Length_cm[DiveNumber=="321"]
pen_papillata321

pen_papillata323<-Length_cm[DiveNumber=="323"]
pen_papillata323

pen_papillata324<-Length_cm[DiveNumber=="324"]
pen_papillata324

pen_papillataNov2011<-c(pen_papillata321,pen_papillata323,pen_papillata324)
pen_papillataNov2011

# This was done to help with the axis
# pen_papillatafix<-c(24)
# 
# pen_papillataNov2011<-c(pen_papillata321,pen_papillata323,pen_papillata324,pen_papillatafix)
# pen_papillataNov2011

h_pen_papillataNov2011<-hist(pen_papillataNov2011, breaks=12,xlim=c(0,25),ylim=c(0,70),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2011 Peniagone papillata (cm)")

summary(pen_papillataNov2011)
#boxplot(pen_papillataNov2011)

#################################################################################

#June 2012 length data for pen_papillataidia
#must make a vector that is only the length measurements from that time period
pen_papillataJune2012<-Length_cm[DiveNumber=="403"]
pen_papillataJune2012

h_pen_papillataJune2012<-hist(pen_papillataJune2012, breaks=12,xlim=c(0,25),ylim=c(0,70),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2012 Peniagone papillata (cm)")


# This was done to help with the axis
# pen_papillatafix2<-c(24)
# 
# pen_papillataJune2012a<-c(pen_papillataJune2012,pen_papillatafix2)
# pen_papillataJune2012a
# h_pen_papillataJune2012<-hist(pen_papillataJune2012a, breaks=12,xlim=c(0,25),ylim=c(0,70),
                         main= "  ", border= "dark blue", col = "light blue",  
                          xlab="June 2012 Peniagone papillata (cm)")


summary(pen_papillataJune2012)
boxplot(pen_papillataJune2012)

#############################################################################

#November 2012 length data 
#must make a vector that is only the length measurements from that time period
pen_papillata442<-Length_cm[DiveNumber=="442"]
pen_papillata442

pen_papillata443<-Length_cm[DiveNumber=="443"]
pen_papillata443

pen_papillataNov2012<-c(pen_papillata442,pen_papillata443)
pen_papillataNov2012
# 
# This was done to help with the axis
# pen_papillatafix<-c(24)
# 
# pen_papillataNov2012a<-c(pen_papillataNov2012,pen_papillatafix)
# pen_papillataNov2012a
# h_pen_papillataNov2012<-hist(pen_papillataNov2012a, breaks=12,xlim=c(0,25),ylim=c(0,70),
                          main= "  ", border= "dark blue", col = "light blue",  
                        xlab="November 2012 Peniagone papillata (cm)")

h_pen_papillataNov2012<-hist(pen_papillataNov2012, breaks=12,xlim=c(0,25),ylim=c(0,70),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2012 Peniagone papillata (cm)")



summary(pen_papillataNov2012)
#boxplot(pen_papillataNov2012)
