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

#this is for all dives
summary(pen_sp_nov)

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

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


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

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

h_pen_sp_novDec2006<-hist(pen_sp_novDec2006, breaks=12,xlim=c(0,20),ylim=c(0,500),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="December 2006 Peniagone sp. nov. (cm)")

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

#This was done to help with the axis
#pen_sp_novfix<-c(14)
# 
#pen_sp_novFeb2007<-c(pen_sp_nov1077,pen_sp_nov1080,pen_sp_novfix)
#pen_sp_novFeb2007

#h_pen_sp_novFeb2007<-hist(pen_sp_novFeb2007, breaks=12,xlim=c(0,20),ylim=c(0,500),
      main= "  ", border= "dark blue", col = "light blue",  
      xlab="December 2006 Peniagone sp. nov. (cm)")

summary(pen_sp_novDec2006)
#boxplot(pen_sp_novDec2006)

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

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

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

pen_sp_novFeb2007<-c(pen_sp_nov1077,pen_sp_nov1080)
pen_sp_novFeb2007

h_pen_sp_novFeb2007<-hist(pen_sp_novFeb2007, breaks=12,xlim=c(0,20),ylim=c(0,500),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2007 Peniagone sp. nov. (cm)")

#This was done to help with the axis
pen_sp_novfix<-c(14)
# 
#pen_sp_novFeb2007<-c(pen_sp_nov1077,pen_sp_nov1080,pen_sp_novfix)
#pen_sp_novFeb2007

#h_pen_sp_novFeb2007<-hist(pen_sp_novFeb2007, breaks=12,xlim=c(0,20),ylim=c(0,500),
main= "  ", border= "dark blue", col = "light blue",  
xlab="February 2007 Peniagone sp. nov. (cm)")

summary(pen_sp_novFeb2007)
#boxplot(pen_sp_novFeb2007)

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

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

h_pen_sp_novJune2007<-hist(pen_sp_novJune2007, breaks=12,xlim=c(0,20),ylim=c(0,500),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2007 Peniagone sp. nov. (cm)")
#This was done to help with the axis
#pen_sp_novfix2<-c(14)
# 
#pen_sp_novJune2007<-c(pen_sp_novJune2007,pen_sp_novfix2)
#pen_sp_novJune2007


summary(pen_sp_novJune2007)
#boxplot(pen_sp_novJune2007)

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

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

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

pen_sp_novSept2007<-c(pen_sp_nov1141,pen_sp_nov1143)
pen_sp_novSept2007

h_pen_sp_novSept2007<-hist(pen_sp_novSept2007, breaks=12,xlim=c(0,20),ylim=c(0,500),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="September 2007 Peniagone sp. nov. (cm)")

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

h_pen_sp_novFeb2007<-hist(pen_sp_novFeb2007, breaks=12,xlim=c(0,20),ylim=c(0,500),
main= "  ", border= "dark blue", col = "light blue",  
xlab="September 2007 Peniagone sp. nov. (cm)")


summary(pen_sp_novSept2007)
#boxplot(pen_sp_novSept2007)

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

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

h_pen_sp_novFeb2009<-hist(pen_sp_novFeb2009, breaks=12,xlim=c(0,20),ylim=c(0,500),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2009 Peniagone sp. nov. (cm)")

# This was done to help with the axis
# pen_sp_novfixf<-c(14)
# 
# pen_sp_novFeb2009<-c(pen_sp_novFeb2009,pen_sp_novfixf)
# pen_sp_novFeb2009
#h_pen_sp_novFeb2009<-hist(pen_sp_novFeb2009, breaks=12,xlim=c(0,20),ylim=c(0,500),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2009 Peniagone sp. nov. (cm)")
summary(pen_sp_novFeb2009)
boxplot(pen_sp_novFeb2009)

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

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

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

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

pen_sp_novMay2011<-c(pen_sp_nov230,pen_sp_nov231,pen_sp_nov232)
pen_sp_novMay2011

h_pen_sp_novMay2011<-hist(pen_sp_novMay2011, breaks=12,xlim=c(0,20),ylim=c(0,500),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="May 2011 Peniagone sp. nov. (cm)")

summary(pen_sp_novMay2011)
#boxplot(pen_sp_novMay2011)

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


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

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

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

pen_sp_novNov2011<-c(pen_sp_nov321,pen_sp_nov323,pen_sp_nov324)
pen_sp_novNov2011



h_pen_sp_novNov2011<-hist(pen_sp_novNov2011, breaks=12,xlim=c(0,20),ylim=c(0,500),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2011 Peniagone sp. nov. (cm)")

summary(pen_sp_novNov2011)
#boxplot(pen_sp_novNov2011)

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

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

h_pen_sp_novJune2012<-hist(pen_sp_novJune2012, breaks=12,xlim=c(0,20),ylim=c(0,500),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2012 Peniagone sp. nov. (cm)")


summary(pen_sp_novJune2012)
boxplot(pen_sp_novJune2012)

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

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

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

pen_sp_novNov2012<-c(pen_sp_nov442,pen_sp_nov443)
pen_sp_novNov2012

h_pen_sp_novNov2012<-hist(pen_sp_novNov2012, breaks=12,xlim=c(0,20),ylim=c(0,500),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2012 Peniagone sp. nov. (cm)")



summary(pen_sp_novNov2012)
#boxplot(pen_sp_novNov2012)