inst/Requests/BZ.LengthWeight.R

In LobsterScience/bio.lobster: Lobster Stock Assessment Tools for DFO Maritimes

##Length Weight Request
require(tidyverse)

p = bio.lobster::load.environment()
#la()
#Plot LW Regressions for Lobster by Sex
 
savdir=file.path(project.datadirectory("bio.lobster","requests","length.weight"))

dir.create( savdir, recursive = TRUE, showWarnings = FALSE ) #Change as Appropriate

#Plot Seperately

sexes=c(1:3)
max.cw=140 #change as appropriate
for (s in sexes){
  
  L=50:max.cw
  B = lobLW(L, sex=s)
  MLW<-cbind(L,B)
  MLW<-as.data.frame(MLW)
  MLW$Weight_lbs<-MLW$B/453.6
  
  names(MLW)[names(MLW) == 'L'] <- 'CarapaceLength_mm'
  names(MLW)[names(MLW) == 'B'] <- 'Weight_g'

  titl=NA
  if (s=="1") {titl="Male Lobster Weight"}
  if (s=="2") {titl="Female Lobster Weight"}
  if (s=="3") {titl="Female Lobster (Berried) Weight"}


png(filename=paste0(savdir,titl,".png"),width=8, height=6.5, units = "in", res = 800)
plot(MLW$CarapaceLength_mm, MLW$Weight_lbs,xlab="Carapace Length (mm)", ylab= "Weight (lbs)", main= titl, type="l")
dev.off()
}


#Plot on one


  L=82.5:110
  B = lobLW(L, sex=1)
  m.MLW<-cbind(L,B)
  m.MLW<-as.data.frame(m.MLW)
  m.MLW$Weight_lbs<-m.MLW$B/453.6
 
 
  F = lobLW(L, sex=2)
  f.MLW<-cbind(L,F)
  f.MLW<-as.data.frame(f.MLW)
  f.MLW$Weight_lbs<-f.MLW$F/453.6
  
 
  b = lobLW(L, sex=3)
  b.MLW<-cbind(L,b)
  b.MLW<-as.data.frame(b.MLW)
  b.MLW$Weight_lbs<-b.MLW$b/453.6
  
  
  titl="Lobster Weights at Length mls change"
  png(filename=paste0(savdir,titl,".png"),width=8, height=6.5, units = "in", res = 800)
  plot(m.MLW$L, m.MLW$Weight_lbs,xlab="Carapace Length (mm)", ylab= "Weight (lbs)", main= titl, type="n", col="black")
  lines(f.MLW$L, f.MLW$Weight_lbs, col="red")
  lines(b.MLW$L, b.MLW$Weight_lbs, col="darkgreen")
  lines(m.MLW$L, m.MLW$Weight_lbs, col="blue")
  text(x=130, y=1.85, col="red", "Female (No Eggs)", pos=4)
  text(x=130, y=1.5, col="darkgreen", "Female (Berried)", pos=4)
  text(x=130, y=1.15, col="blue", "Male", pos=4)
  #abline(v = 82.5, col = 'black', lwd = 2, lty = 'dashed')
  #abline(v = 86, col = 'black', lwd = 2, lty = 'dashed')
  #rect(xleft = 82.5, xright = 86, ybottom = 0, ytop = 1.35, col="darkred", density=30, angle=-30, lwd=1)
  dev.off()

  
 #bring in (limited) 2024 port sampling data
 l24=read.csv(file.path(savdir, "master.nov.24.csv"))
 l24=l24[,colSums(is.na(l24))<nrow(l24)]
 colnames(l24) <- tolower(colnames(l24))
 l24= l24 %>% rename("L" = "carapace", "B" = "wt.g.final")
 
 l24=subset(l24, B< 3000 & B>100 & L>80 & L<200 )
 
#males
 m24=subset(l24, sex=="1")
 lo.m=loess(B~L, m24)
 l=order(m24$L)
 plot(m.MLW$L, m.MLW$B,xlab="Carapace Length (mm)", ylab= "Weight (grams)", main= "", type="l", col="black")
 lines(m24$L[l],lo.m$fitted[l],col="red",lwd=1)
 
#females
 f24=subset(l24, sex=="2")
 lo.f=loess(B~L, f24)
 lf=order(f24$L)
 plot(f.MLW$L, f.MLW$F,xlab="Carapace Length (mm)", ylab= "Weight (grams)", main= "Female", type="l", col="black")
 lines(f24$L[lf],lo.f$fitted[lf],col="red",lwd=1)