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

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

h_elp_length<-hist(Length_cm, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="Elpidia sp. nov. (cm)")

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

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

h_elpDec2006<-hist(elpDec2006, breaks=12,xlim=c(0,6),ylim=c(0,800),
               main= "  ", border= "dark blue", col = "light blue",  
               xlab="December 2006 Elpidia sp. nov. (cm)")

# if you want to zoom in on the values to know where they are (Close-up)

h_elpDec2006<-hist(elpDec2006, breaks=12,xlim=c(0,6),ylim=c(0,10),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="December 2006 Elpidia sp. nov. (cm)")

summary(elpDec2006)
#boxplot(elpDec2006)

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

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

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

elpFeb2007<-c(elp1077,elp1080)
elpFeb2007

h_elpFeb2007<-hist(elpFeb2007, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="February 2007 Elpidia sp. nov. (cm)")

# if you want to zoom in on the values to know where they are (Close-up)

h_elpFeb2007<-hist(elpFeb2007, breaks=12,xlim=c(0,6),ylim=c(0,10),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="February 2007 Elpidia sp. nov. (cm)")

summary(elpFeb2007)
#boxplot(elpFeb2007)

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

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

h_elpJune2007<-hist(elpJune2007, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="June 2007 Elpidia sp. nov. (cm)")

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

elpDec2006<-Length_cm[DiveNumber=="1067"]
elpDec2006

h_elpDec2006<-hist(elpDec2006, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="June 2007 Elpidia sp. nov. (cm)")

summary(elpJune2007)
#boxplot(elpJune2007)

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

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

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

elpSept2007<-c(elp1141,elp1143)
elpSept2007

h_elpSept2007<-hist(elpSept2007, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="September 2007 Elpidia sp. nov. (cm)")

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

elpDec2006<-Length_cm[DiveNumber=="1067"]
elpDec2006

h_elpDec2006<-hist(elpDec2006, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="September 2007 Elpidia sp. nov. (cm)")

summary(elpSept2007)
#boxplot(elpSept2007)

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

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

h_elpFeb2009<-hist(elpFeb2009, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="February 2009 Elpidia sp. nov. (cm)")

# if you want to zoom in on the values to know where they are (Close-up)

h_elpFeb2009<-hist(elpFeb2009, breaks=12,xlim=c(0,6),ylim=c(0,10),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="February 2009 Elpidia sp. nov. (cm)")


summary(elpFeb2009)
boxplot(elpFeb2009)

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

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

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

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

elpMay2011<-c(elp230,elp231,elp232)
elpMay2011

h_elpMay2011<-hist(elpMay2011, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="May 2011 Elpidia sp. nov. (cm)")

summary(elpMay2011)
#boxplot(elpMay2011)

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


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

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

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

elpNov2011<-c(elp321,elp323,elp324)
elpNov2011

h_elpNov2011<-hist(elpNov2011, breaks=12,xlim=c(0,6),ylim=c(0,800),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="November 2011 Elpidia sp. nov. (cm)")

summary(elpNov2011)
#boxplot(elpNov2011)

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

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

h_elpJune2012<-hist(elpJune2012, breaks=12,xlim=c(0,6),ylim=c(0,800),
               main= "  ", border= "dark blue", col = "light blue",  
               xlab="June 2012 Elpidia sp. nov. (cm)")

# if you want to zoom in on the values to know where they are (Close-up)

h_elpJune2012<-hist(elpJune2012, breaks=12,xlim=c(0,6),ylim=c(0,100),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="June2012 Elpidia sp. nov. (cm)")


summary(elpJune2012)
boxplot(elpJune2012)

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

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

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

elpNov2012<-c(elp442,elp443)
elpNov2012

h_elpNov2012<-hist(elpNov2012, breaks=12,xlim=c(0,6),ylim=c(0,800),
                    main= "  ", border= "dark blue", col = "light blue",  
                    xlab="November 2012 Elpidia sp. nov. (cm)")


summary(elpNov2012)
#boxplot(elpNov2012)