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

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

h_psychro_length<-hist(Length_cm, breaks=12,xlim=c(0,15),ylim=c(0,5),
                       main= "  ", border= "dark blue", col = "light blue",  
                       xlab="Psychropotes longicauda Length (cm)")
########################################################################

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

h_psychroDec2006<-hist(psychroDec2006, breaks=12,xlim=c(0,15),ylim=c(0,5),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="December 2006 Psychropotes longicauda (cm)")

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


h_psychroJune2012<-hist(psychroJune2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="December 2006 Psychropotes longicauda (cm)")

summary(psychroDec2006)
#boxplot(psychroDec2006)

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

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

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

psychroFeb2007<-c(psychro1077,psychro1080)
psychroFeb2007

h_psychroFeb2007<-hist(psychroFeb2007, breaks=12,xlim=c(0,15),ylim=c(0,5),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="February 2007 Psychropotes longicauda (cm)")


#Because there are no observations in this time period for Psychropotes longicauda, use
#June2012 graph below with June 2007 labels so I can erase the line in Powerpoint
                        
                        
h_psychroJune2012<-hist(psychroJune2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="February 2007 Psychropotes longicauda (cm)")

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

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

h_psychroJune2007<-hist(psychroJune2007, breaks=12,xlim=c(0,15),ylim=c(0,5),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2007 Psychropotes longicauda (cm)")

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

h_psychroFeb2007<-hist(psychroFeb2007, breaks=12,xlim=c(0,15),ylim=c(0,5),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="June 2007 Psychropotes longicauda (cm)")

 #Because there are no observations in this time period for Psychropotes longicauda, use
 #June2012 graph below with June 2007 labels so I can erase the line in Powerpoint
                          
                          
h_psychroJune2012<-hist(psychroJune2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="June 2007 Psychropotes longicauda (cm)")
                         
summary(psychroJune2007)
#boxplot(psychroJune2007)

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

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

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

psychroSept2007<-c(psychro1141,psychro1143)
psychroSept2007

h_psychroSept2007<-hist(psychroSept2007, breaks=12,xlim=c(0,15),ylim=c(0,5),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="September 2007 Psychropotes longicauda (cm)")

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


h_psychroJune2012<-hist(psychroJune2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="September 2007 Psychropotes longicauda (cm)")


summary(psychroSept2007)
#boxplot(psychroSept2007)

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

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

h_psychroFeb2009<-hist(psychroFeb2009, breaks=12,xlim=c(0,15),ylim=c(0,5),
                  main= "  ", border= "dark blue", col = "light blue",  
                  xlab="February 2009 Psychropotes longicauda (cm)")

#Because there are no observations in this time period for Psychropotes longicauda, use
#June2012 graph below with June 2007 labels so I can erase the line in Powerpoint
                         
                         
h_psychroJune2012<-hist(psychroJune2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
               main= "  ", border= "dark blue", col = "light blue",  
               xlab="February 2009 Psychropotes longicauda (cm)")
                         
summary(psychroFeb2009)
boxplot(psychroFeb2009)

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

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

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

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

psychroMay2011<-c(psychro230,psychro231,psychro232)
psychroMay2011

h_psychroMay2011<-hist(psychroMay2011, breaks=12,xlim=c(0,15),ylim=c(0,5),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="May 2011 Psychropotes longicauda (cm)")

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


h_psychroJune2012<-hist(psychroJune2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="May 2011 Psychropotes longicauda (cm)")

summary(psychroMay2011)
#boxplot(psychroMay2011)

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


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

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

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

psychroNov2011<-c(psychro321,psychro323,psychro324)
psychroNov2011

h_psychroNov2011<-hist(psychroNov2011, breaks=12,xlim=c(0,15),ylim=c(0,5),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2011 Psychropotes longicauda (cm)")

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


h_psychroJune2012<-hist(psychroJune2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
                        main= "  ", border= "dark blue", col = "light blue",  
                        xlab="November 2011 Psychropotes longicauda (cm)")

summary(psychroNov2011)
#boxplot(psychroNov2011)

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

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

h_psychroJune2012<-hist(psychroJune2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
                           main= "  ", border= "dark blue", col = "light blue",  
                           xlab="June 2012 Psychropotes longicauda (cm)")


# This was done to help with the axis
# psychrofix2<-c(0)
# 
# psychroJune2012a<-c(psychroJune2012,psychrofix2)
# psychroJune2012a
h_psychroJune2012<-hist(psychroJune2012a, breaks=12,xlim=c(0,15),ylim=c(0,5),
                         main= "  ", border= "dark blue", col = "light blue",  
                         xlab="June 2012 Psychropotes longicauda (cm)")


summary(psychroJune2012)
boxplot(psychroJune2012)

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

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

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

psychroNov2012<-c(psychro442,psychro443)
psychroNov2012
# 
# This was done to help with the axis
# psychrofix2<-c(14)
# 
# psychroNov2012a<-c(psychroNov2012,psychrofix2)
# psychroNov2012a
# h_psychroNov2012<-hist(psychroNov2012a, breaks=12,xlim=c(0,15),ylim=c(0,5),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2012 Psychropotes longicauda (cm)")

h_psychroNov2012<-hist(psychroNov2012, breaks=12,xlim=c(0,15),ylim=c(0,5),
                          main= "  ", border= "dark blue", col = "light blue",  
                          xlab="November 2012 Psychropotes longicauda (cm)")



summary(psychroNov2012)
#boxplot(psychroNov2012)