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

#this is for all dives
summary(pen_gracilis)
#plot(Biomass~Length_cm)

h_pen_gracilis_length<-hist(Length_cm, breaks=12,xlim=c(0,10),ylim=c(0,150),
                            main= "  ", border= "dark blue", col = "light blue",  
                            xlab="Peniagone gracilis (cm)")

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

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

h_pen_gracDec2006<-hist(pen_gracDec2006, breaks=12,xlim=c(0,10),ylim=c(0,150),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="December 2006 Peniagone gracilis (cm)")

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

h_pen_gracJune2012<-hist(pen_gracJune2012, breaks=12,xlim=c(0,10),ylim=c(0,150),
                         main= "  ", border= "dark blue", col = "light blue",  
                         xlab="December 2006 Peniagone gracilis (cm)")

summary(pen_gracDec2006)
#boxplot(pen_gracDec2006)

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

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

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

pen_gracFeb2007<-c(pen_grac1077,pen_grac1080)
pen_gracFeb2007

h_pen_gracFeb2007<-hist(pen_gracFeb2007, breaks=12,xlim=c(0,10),ylim=c(0,150),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="February 2007 Peniagone gracilis (cm)")

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

h_pen_gracJune2012<-hist(pen_gracJune2012, breaks=12,xlim=c(0,10),ylim=c(0,150),
                         main= "  ", border= "dark blue", col = "light blue",  
                         xlab="February 2007 Peniagone gracilis (cm)")


summary(pen_gracFeb2007)
#boxplot(pen_gracFeb2007)

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

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

h_pen_gracJune2007<-hist(pen_gracJune2007, breaks=12,xlim=c(0,10),ylim=c(0,150),
                    main= "  ", border= "dark blue", col = "light blue",  
                    xlab="June 2007 Peniagone gracilis (cm)")

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

h_pen_gracJune2012<-hist(pen_gracJune2012, breaks=12,xlim=c(0,10),ylim=c(0,150),
                         main= "  ", border= "dark blue", col = "light blue",  
                         xlab="June 2007 Peniagone gracilis (cm)")

summary(pen_gracJune2007)
#boxplot(pen_gracJune2007)

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

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

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

pen_gracSept2007<-c(pen_grac1141,pen_grac1143)
pen_gracSept2007

h_pen_gracSept2007<-hist(pen_gracSept2007, breaks=12,xlim=c(0,10),ylim=c(0,150),
                    main= "  ", border= "dark blue", col = "light blue",  
                    xlab="September 2007 Peniagone gracilis (cm)")

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

h_pen_gracJune2012<-hist(pen_gracJune2012, breaks=12,xlim=c(0,10),ylim=c(0,150),
                         main= "  ", border= "dark blue", col = "light blue",  
                         xlab="September 2007 Peniagone gracilis (cm)")


summary(pen_gracSept2007)
#boxplot(pen_gracSept2007)

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

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

h_pen_gracFeb2009<-hist(pen_gracFeb2009, breaks=12,xlim=c(0,10),ylim=c(0,150),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="February 2009 Peniagone gracilis (cm)")

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

h_pen_gracJune2012<-hist(pen_gracJune2012, breaks=12,xlim=c(0,10),ylim=c(0,150),
                         main= "  ", border= "dark blue", col = "light blue",  
                         xlab="February 2009 Peniagone gracilis (cm)")


summary(pen_gracFeb2009)
boxplot(pen_gracFeb2009)

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

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

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

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

pen_gracMay2011<-c(pen_grac230,pen_grac231,pen_grac232)
pen_gracMay2011

h_pen_gracMay2011<-hist(pen_gracMay2011, breaks=12,xlim=c(0,10),ylim=c(0,150),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="May 2011 Peniagone gracilis (cm)")

summary(pen_gracMay2011)
#boxplot(pen_gracMay2011)

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


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

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

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

pen_gracNov2011<-c(pen_grac321,pen_grac323,pen_grac324)
pen_gracNov2011

#This was done to help with the axis
#pen_gracfix<-c(9)

#pen_gracNov2011<-c(pen_grac321,pen_grac323,pen_grac324,pen_gracfix)
#pen_gracNov2011

h_pen_gracNov2011<-hist(pen_gracNov2011, breaks=12,xlim=c(0,10),ylim=c(0,150),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="November 2011 Peniagone gracilis (cm)")

summary(pen_gracNov2011)
#boxplot(pen_gracNov2011)

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

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


h_pen_gracJune2012<-hist(pen_gracJune2012, breaks=12,xlim=c(0,10),ylim=c(0,150),
                    main= "  ", border= "dark blue", col = "light blue",  
                    xlab="June 2012 Peniagone gracilis (cm)")


summary(pen_gracJune2012)
boxplot(pen_gracJune2012)

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

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

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

pen_gracNov2012<-c(pen_grac442,pen_grac443)
pen_gracNov2012



h_pen_gracNov2012<-hist(pen_gracNov2012, breaks=12,xlim=c(0,10),ylim=c(0,150),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="November 2012 Peniagone gracilis (cm)")

#This was done to help with the axis
# pen_gracfix2<-c(9)
# 
# pen_gracNov2012<-c(pen_grac321,pen_grac323,pen_grac324,pen_gracfix2)
# pen_gracNov2012
# 
# h_pen_gracNov2012<-hist(pen_gracNov2012, breaks=12,xlim=c(0,10),ylim=c(0,150),
#                         main= "  ", border= "dark blue", col = "light blue",  
#                         xlab="November 2012 Peniagone gracilis (cm)")


summary(pen_gracNov2012)
#boxplot(pen_gracNov2012)