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

#this is for all dives
summary(abyss)

h_abyss_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,10),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="Abyssocucumis abyssorum (cm)")


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

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

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

h_abyss_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="December 2006 Abyssocucumis abyssorum (cm)")


summary(abyssDec2006)
#boxplot(abyssDec2006)

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

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

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

abyssFeb2007<-c(abyss1077,abyss1080)
abyssFeb2007

h_abyssFeb2007<-hist(abyssFeb2007, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2007 Abyssocucumis abyssorum (cm)")

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

h_abyss_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2007 Abyssocucumis abyssorum (cm)")

summary(abyssFeb2007)
#boxplot(abyssFeb2007)

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

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

h_abyssJune2007<-hist(abyssJune2007, breaks=12,xlim=c(0,15),ylim=c(0,10),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="June 2007 Abyssocucumis abyssorum (cm)")

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

h_abyss_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="June 2007 Abyssocucumis abyssorum (cm)")

summary(abyssJune2007)
#boxplot(abyssJune2007)

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

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

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

abyssSept2007<-c(abyss1141,abyss1143)
abyssSept2007

h_abyssSept2007<-hist(abyssSept2007, breaks=12,xlim=c(0,15),ylim=c(0,10),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="September 2007 Abyssocucumis abyssorum (cm)")

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

h_abyss_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="September 2007 Abyssocucumis abyssorum (cm)")


summary(abyssSept2007)
#boxplot(abyssSept2007)

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

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

h_abyssFeb2009<-hist(abyssFeb2009, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2009 Abyssocucumis abyssorum (cm)")

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

h_abyss_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2009 Abyssocucumis abyssorum (cm)")

summary(abyssFeb2009)
boxplot(abyssFeb2009)

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

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

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

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

abyssMay2011<-c(abyss230,abyss231,abyss232)
abyssMay2011

h_abyssMay2011<-hist(abyssMay2011, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="May 2011 Abyssocucumis abyssorum (cm)")

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

h_abyss_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="May 2011 Abyssocucumis abyssorum (cm)")

summary(abyssMay2011)
#boxplot(abyssMay2011)

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


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

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

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

abyssNov2011<-c(abyss321,abyss323,abyss324)
abyssNov2011

h_abyssNov2011<-hist(abyssNov2011, breaks=12,xlim=c(0,15),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="November 2011 Abyssocucumis abyssorum (cm)")

#this is done to fix the axis
#abyssfix<-c(0,14)
#abyssNov2011a<-c(abyssNov2011,abyssfix)

#h_abyssNov2011<-hist(abyssNov2011a, breaks=12,xlim=c(0,15),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="November 2011 Abyssocucumis abyssorum (cm)")

summary(abyssNov2011)
#boxplot(abyssNov2011)

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

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

h_abyssJune2012<-hist(abyssJune2012, breaks=12,xlim=c(0,15),ylim=c(0,10),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="June 2012 Abyssocucumis abyssorum (cm)")


summary(abyssJune2012)
boxplot(abyssJune2012)

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

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

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

abyssNov2012<-c(abyss442,abyss443)
abyssNov2012


h_abyssNov2012<-hist(abyssNov2012, breaks=12,xlim=c(0,15),ylim=c(0,10),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="November 2012 Abyssocucumis abyssorum (cm)")

#this is done to fix the axis
#abyssfix1<-c(0)
#abyssNov2012a<-c(abyssNov2012,abyssfix1)
#abyssNov2012a

#h_abyssNov2012<-hist(abyssNov2012a, breaks=12,xlim=c(0,15),ylim=c(0,10),
                     main= "  ", border= "dark blue", col = "light blue",  
                     xlab="November 2012 Abyssocucumis abyssorum (cm)")

summary(abyssNov2012)
#boxplot(abyssNov2012)