scoto<-read.table(file.choose(),sep=",",header=TRUE)
scoto
attach(scoto)

#this is for all dives
summary(scoto)
#plot(Biomass~Length_cm)

h_scoto_length<-hist(Length_cm, breaks=12,xlim=c(0,20),ylim=c(0,150),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="Scotoplanes globosa Length (cm)")

########################################################################

#December 2006 length data 
#must make a vector that is only the length measurements from that time period
scotoDec2006<-Length_cm[DiveNumber=="1067"]
scotoDec2006

h_scotoDec2006<-hist(scotoDec2006, breaks=12,xlim=c(0,20),ylim=c(0,150),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="December 2006 Scotoplanes globosa (cm)")

#Because there are no observations in this time period for Scotoplanes globosa, use
#xxxx graph below with June 2007 labels so I can erase the line in Powerpoint

h_scotoNov2012<-hist(scotoNov2012, breaks=12,xlim=c(0,20),ylim=c(0,150),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="December 2006 Scotoplanes globosa (cm)")


summary(scotoDec2006)
#boxplot(scotoDec2006)

#############################################################################

#February 2007 length data 
#must make a vector that is only the length measurements from that time period
scoto1077<-Length_cm[DiveNumber=="1077"]
scoto1077

scoto1080<-Length_cm[DiveNumber=="1080"]
scoto1080

scotoFeb2007<-c(scoto1077,scoto1080)
scotoFeb2007

h_scotoFeb2007<-hist(scotoFeb2007, breaks=12,xlim=c(0,20),ylim=c(0,150),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2007 Scotoplanes globosa (cm)")

#Because there are no observations in this time period for Scotoplanes globosa, use
#xxxx graph below with February 2007 labels so I can erase the line in Powerpoint

h_scotoNov2012<-hist(scotoNov2012, breaks=12,xlim=c(0,20),ylim=c(0,150),
                  main= "  ", border= "dark blue", col = "light blue",  
                  xlab="February 2007 Scotoplanes globosa (cm)")



summary(scotoFeb2007)
#boxplot(scotoFeb2007)

###########################################################################

#June 2007 length data 
#must make a vector that is only the length measurements from that time period
scotoJune2007<-Length_cm[DiveNumber=="1094"]
scotoJune2007

h_scotoJune2007<-hist(scotoJune2007, breaks=12,xlim=c(0,20),ylim=c(0,150),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2007 Scotoplanes globosa (cm)")

#Because there are no observations in this time period for Scotoplanes globosa, use
#xxxx graph below with June 2007 labels so I can erase the line in Powerpoint

h_scotoNov2012<-hist(scotoNov2012, breaks=12,xlim=c(0,20),ylim=c(0,150),
                  main= "  ", border= "dark blue", col = "light blue",  
                  xlab="June 2007 Scotoplanes globosa (cm)")

summary(scotoJune2007)
#boxplot(scotoJune2007)

############################################################################

#Sept 2007 length data 
#must make a vector that is only the length measurements from that time period
scoto1141<-Length_cm[DiveNumber=="1141"]
scoto1141

scoto1143<-Length_cm[DiveNumber=="1143"]
scoto1143

scotoSept2007<-c(scoto1141,scoto1143)
scotoSept2007

h_scotoSept2007<-hist(scotoSept2007, breaks=12,xlim=c(0,20),ylim=c(0,150),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="September 2007 Scotoplanes globosa (cm)")

#Because there are no observations in this time period for Scotoplanes globosa, use
#xxxx graph below with June 2007 labels so I can erase the line in Powerpoint

h_scotoNov2012<-hist(scotoNov2012, breaks=12,xlim=c(0,20),ylim=c(0,150),
                  main= "  ", border= "dark blue", col = "light blue",  
                  xlab="September 2007 Scotoplanes globosa (cm)")


summary(scotoSept2007)
#boxplot(scotoSept2007)

##############################################################################

#February 2009 length data
#must make a vector that is only the length measurements from that time period
scotoFeb2009<-Length_cm[DiveNumber=="8"]
scotoFeb2009

h_scotoFeb2009<-hist(scotoFeb2009, breaks=12,xlim=c(0,20),ylim=c(0,150),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2009 Scotoplanes globosa (cm)")

#Because there are no observations in this time period for Scotoplanes globosa, use
#xxxx graph below with June 2007 labels so I can erase the line in Powerpoint

h_scotoNov2012<-hist(scotoNov2012, breaks=12,xlim=c(0,20),ylim=c(0,150),
                  main= "  ", border= "dark blue", col = "light blue",  
                  xlab="February 2009 Scotoplanes globosa (cm)")


summary(scotoFeb2009)
boxplot(scotoFeb2009)

###############################################################################

#May 2011 length data 
#must make a vector that is only the length measurements from that time period
scoto230<-Length_cm[DiveNumber=="230"]
scoto230

scoto231<-Length_cm[DiveNumber=="231"]
scoto231

scoto232<-Length_cm[DiveNumber=="232"]
scoto232

scotoMay2011<-c(scoto230,scoto231,scoto232)
scotoMay2011

h_scotoMay2011<-hist(scotoMay2011, breaks=12,xlim=c(0,20),ylim=c(0,150),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="May 2011 Scotoplanes globosa (cm)")

summary(scotoMay2011)
#boxplot(scotoMay2011)

################################################################################


#Nov 2011 length data 
#must make a vector that is only the length measurements from that time period
scoto321<-Length_cm[DiveNumber=="321"]
scoto321

scoto323<-Length_cm[DiveNumber=="323"]
scoto323

scoto324<-Length_cm[DiveNumber=="324"]
scoto324

scotoNov2011<-c(scoto321,scoto323,scoto324)
scotoNov2011


h_scotoNov2011<-hist(scotoNov2011, breaks=12,xlim=c(0,20),ylim=c(0,150),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2011 Scotoplanes globosa (cm)")

summary(scotoNov2011)
#boxplot(scotoNov2011)

#################################################################################

#June 2012 length data for scotoidia
#must make a vector that is only the length measurements from that time period
scotoJune2012<-Length_cm[DiveNumber=="403"]
scotoJune2012

h_scotoJune2012<-hist(scotoJune2012, breaks=12,xlim=c(0,20),ylim=c(0,150),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2012 Scotoplanes globosa (cm)")




summary(scotoJune2012)
boxplot(scotoJune2012)

#############################################################################

#November 2012 length data 
#must make a vector that is only the length measurements from that time period
scoto442<-Length_cm[DiveNumber=="442"]
scoto442

scoto443<-Length_cm[DiveNumber=="443"]
scoto443

scotoNov2012<-c(scoto442,scoto443)
scotoNov2012

h_scotoNov2012<-hist(scotoNov2012, breaks=12,xlim=c(0,20),ylim=c(0,150),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2012 Scotoplanes globosa (cm)")



summary(scotoNov2012)
#boxplot(scotoNov2012)