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

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

h_phanta_length<-hist(Length_cm, breaks=10,xlim=c(0,20),ylim=c(0,5),
                      main= "  ", border= "dark blue", col = "light blue",  
                      xlab="Oneirophanta mutabilis complex Length (cm)")

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

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

h_phantaDec2006<-hist(phantaDec2006, breaks=10,xlim=c(0,15),ylim=c(0,5),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="December 2006 Oneirophanta mutabilis complex (cm)")

#Because there are no observations in this time period for Oneirophanta, use
#June 2012 graph below with June 2007 labels so I can erase the line in Powerpoint
phantaJune2012<-Length_cm[DiveNumber=="403"]
phantaJune2012
h_phantaJune2012<-hist(phantaJune2012, breaks=10,xlim=c(0,20),ylim=c(0,5),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="December 2006 Oneirophanta mutabilis complex (cm)")

summary(phantaDec2006)
#boxplot(phantaDec2006)

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

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

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

phantaFeb2007<-c(phanta1077,phanta1080)
phantaFeb2007

h_phantaFeb2007<-hist(phantaFeb2007, breaks=10,xlim=c(0,15),ylim=c(0,5),
                      main= "  ", border= "dark blue", col = "light blue",  
                      xlab="February 2007 Oneirophanta mutabilis complex (cm)")

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

phantaJune2012<-Length_cm[DiveNumber=="403"]
phantaJune2012
h_phantaJune2012<-hist(phantaJune2012, breaks=10,xlim=c(0,20),ylim=c(0,5),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="February 2007 Oneirophanta mutabilis complex (cm)")


summary(phantaFeb2007)
#boxplot(phantaFeb2007)

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

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

h_phantaJune2007<-hist(phantaJune2007, breaks=10,xlim=c(0,20),ylim=c(0,5),
                    main= "  ", border= "dark blue", col = "light blue",  
                    xlab="June 2007 Oneirophanta mutabilis complex (cm)")

#Because there are no observations in this time period for Oneirophanta, use
#June 2012 graph below with June 2007 labels so I can erase the line in Powerpoint
phantaJune2012<-Length_cm[DiveNumber=="403"]
phantaJune2012
h_phantaJune2012<-hist(phantaJune2012, breaks=10,xlim=c(0,20),ylim=c(0,5),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="June 2007 Oneirophanta mutabilis complex (cm)")

summary(phantaJune2007)
#boxplot(phantaJune2007)

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

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

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

phantaSept2007<-c(phanta1141,phanta1143)
phantaSept2007

h_phantaSept2007<-hist(phantaSept2007, breaks=10,xlim=c(0,20),ylim=c(0,5),
                    main= "  ", border= "dark blue", col = "light blue",  
                    xlab="September 2007 Oneirophanta mutabilis complex (cm)")

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

phantaJune2012<-Length_cm[DiveNumber=="403"]
phantaJune2012
h_phantaJune2012<-hist(phantaJune2012, breaks=10,xlim=c(0,20),ylim=c(0,5),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="September 2007 Oneirophanta mutabilis complex (cm)")

summary(phantaSept2007)
#boxplot(phantaSept2007)

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

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

h_phantaFeb2009<-hist(phantaFeb2009, breaks=10,xlim=c(0,20),ylim=c(0,5),
               main= "  ", border= "dark blue", col = "light blue",  
               xlab="February 2009 Oneirophanta mutabilis complex (cm)")


summary(phantaFeb2009)
boxplot(phantaFeb2009)

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

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

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

# This was done to make a better axis then taken out in Powerpoint
#phanta232<-c(0,20)
#phanta232

phantaMay2011<-c(phanta230,phanta231,phanta232)
phantaMay2011

h_phantaMay2011<-hist(phantaMay2011, breaks=10,xlim=c(0,20),ylim=c(0,5),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="May 2011 Oneirophanta mutabilis complex (cm)")

summary(phantaMay2011)
#boxplot(phantaMay2011)

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


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

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

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

phantaNov2011<-c(phanta321,phanta323,phanta324)
phantaNov2011

h_phantaNov2011<-hist(phantaNov2011, breaks=10,xlim=c(0,20),ylim=c(0,5),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="November 2011 Oneirophanta mutabilis complex (cm)")

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

phantaJune2012<-Length_cm[DiveNumber=="403"]
phantaJune2012
h_phantaJune2012<-hist(phantaJune2012, breaks=10,xlim=c(0,20),ylim=c(0,5),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="November 2011 Oneirophanta mutabilis complex (cm)")
summary(phantaNov2011)
#boxplot(phantaNov2011)

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

#June 2012 length data for Oneirophanta
#must make a vector that is only the length measurements from that time period
phantaJune2012<-Length_cm[DiveNumber=="403"]
phantaJune2012
h_phantaJune2012<-hist(phantaJune2012, breaks=10,xlim=c(0,20),ylim=c(0,5),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="June 2012 Oneirophanta mutabilis complex (cm)")

#Below was done to have a better axis then fixed in powerpoint on the pic
#phantaadd<-c(20)
#phantaadd

#phantafix<-c(phantaJune2012,phantaadd)


#h_phantaJune2012<-hist(phantafix, breaks=10,xlim=c(0,20),ylim=c(0,5),
                    main= "  ", border= "dark blue", col = "light blue",  
                    xlab="June 2012 Oneirophanta mutabilis complex (cm)")


summary(phantaJune2012)
boxplot(phantaJune2012)

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

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

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

phantaNov2012<-c(phanta442,phanta443)
phantaNov2012

h_phantaNov2012<-hist(phantaNov2012, breaks=10,xlim=c(0,20),ylim=c(0,5),
                   main= "  ", border= "dark blue", col = "light blue",  
                   xlab="November 2012 Oneirophanta mutabilis complex (cm)")

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

phantaJune2012<-Length_cm[DiveNumber=="403"]
phantaJune2012
h_phantaJune2012<-hist(phantaJune2012, breaks=10,xlim=c(0,20),ylim=c(0,5),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="November 2012 Oneirophanta mutabilis complex (cm)")

summary(phantaNov2012)
#boxplot(phantaNov2012)