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

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

h_pen_vitrea_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,250),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="Peniagone vitrea Length (cm)")

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

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

h_pen_vitreaDec2006<-hist(pen_vitreaDec2006, breaks=12,xlim=c(0,15),ylim=c(0,250),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="December 2006 Peniagone vitrea (cm)")


summary(pen_vitreaDec2006)
#boxplot(pen_vitreaDec2006)

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

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

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

pen_vitreaFeb2007<-c(pen_vitrea1077,pen_vitrea1080)
pen_vitreaFeb2007

h_pen_vitreaFeb2007<-hist(pen_vitreaFeb2007, breaks=12,xlim=c(0,15),ylim=c(0,250),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="February 2007 Peniagone vitrea (cm)")

#This was done to help with the axis
pen_vitreafix<-c(14)
# 
#pen_vitreaFeb2007<-c(pen_vitrea1077,pen_vitrea1080,pen_vitreafix)
#pen_vitreaFeb2007

#h_pen_vitreaFeb2007<-hist(pen_vitreaFeb2007, breaks=12,xlim=c(0,15),ylim=c(0,250),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2007 Peniagone vitrea (cm)")
summary(pen_vitreaFeb2007)
#boxplot(pen_vitreaFeb2007)

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

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

h_pen_vitreaJune2007<-hist(pen_vitreaJune2007, breaks=12,xlim=c(0,15),ylim=c(0,250),
                         main= "  ", border= "dark blue", col = "light blue",  
                         xlab="June 2007 Peniagone vitrea (cm)")

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

h_pen_vitreaFeb2007<-hist(pen_vitreaFeb2007, breaks=12,xlim=c(0,15),ylim=c(0,250),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="June 2007 Peniagone vitrea (cm)")

summary(pen_vitreaJune2007)
#boxplot(pen_vitreaJune2007)

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

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

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

pen_vitreaSept2007<-c(pen_vitrea1141,pen_vitrea1143)
pen_vitreaSept2007

h_pen_vitreaSept2007<-hist(pen_vitreaSept2007, breaks=12,xlim=c(0,15),ylim=c(0,250),
                         main= "  ", border= "dark blue", col = "light blue",  
                         xlab="September 2007 Peniagone vitrea (cm)")

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

h_pen_vitreaFeb2007<-hist(pen_vitreaFeb2007, breaks=12,xlim=c(0,15),ylim=c(0,250),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="September 2007 Peniagone vitrea (cm)")


summary(pen_vitreaSept2007)
#boxplot(pen_vitreaSept2007)

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

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

h_pen_vitreaFeb2009<-hist(pen_vitreaFeb2009, breaks=12,xlim=c(0,15),ylim=c(0,250),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="February 2009 Peniagone vitrea (cm)")


summary(pen_vitreaFeb2009)
boxplot(pen_vitreaFeb2009)

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

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

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

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

pen_vitreaMay2011<-c(pen_vitrea230,pen_vitrea231,pen_vitrea232)
pen_vitreaMay2011

h_pen_vitreaMay2011<-hist(pen_vitreaMay2011, breaks=12,xlim=c(0,15),ylim=c(0,250),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="May 2011 Peniagone vitrea (cm)")

summary(pen_vitreaMay2011)
#boxplot(pen_vitreaMay2011)

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


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

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

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

pen_vitreaNov2011<-c(pen_vitrea321,pen_vitrea323,pen_vitrea324)
pen_vitreaNov2011

# This was done to help with the axis
# pen_vitreafix<-c(14)
# 
# pen_vitreaNov2011<-c(pen_vitrea321,pen_vitrea323,pen_vitrea324,pen_vitreafix)
# pen_vitreaNov2011

h_pen_vitreaNov2011<-hist(pen_vitreaNov2011, breaks=12,xlim=c(0,15),ylim=c(0,250),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="November 2011 Peniagone vitrea (cm)")

summary(pen_vitreaNov2011)
#boxplot(pen_vitreaNov2011)

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

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

h_pen_vitreaJune2012<-hist(pen_vitreaJune2012, breaks=12,xlim=c(0,15),ylim=c(0,250),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2012 Peniagone vitrea (cm)")


# This was done to help with the axis
# pen_vitreafix2<-c(14)
# 
# pen_vitreaJune2012a<-c(pen_vitreaJune2012,pen_vitreafix2<-c(14))
# pen_vitreaJune2012a
# h_pen_vitreaJune2012<-hist(pen_vitreaJune2012a, breaks=12,xlim=c(0,15),ylim=c(0,250),
#                          main= "  ", border= "dark blue", col = "light blue",  
#                          xlab="June 2012 Peniagone vitrea (cm)")


summary(pen_vitreaJune2012)
boxplot(pen_vitreaJune2012)

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

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

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

# pen_vitreaNov2012<-c(pen_vitrea442,pen_vitrea443)
# pen_vitreaNov2012
# 
# This was done to help with the axis
# pen_vitreafix<-c(14)
# 
# pen_vitreaNov2012a<-c(pen_vitreaNov2012,pen_vitreafix)
# pen_vitreaNov2012a
# h_pen_vitreaJune2012<-hist(pen_vitreaJune2012a, breaks=12,xlim=c(0,15),ylim=c(0,250),
#                          main= "  ", border= "dark blue", col = "light blue",  
#                          xlab="November 2012 Peniagone vitrea (cm)")

h_pen_vitreaNov2012<-hist(pen_vitreaNov2012, breaks=12,xlim=c(0,15),ylim=c(0,250),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="November 2012 Peniagone vitrea (cm)")



summary(pen_vitreaNov2012)
#boxplot(pen_vitreaNov2012)