R/mean.wt.catch.R

In ices-tools-dev/NephAssess: SG Nephrops assessment functions

mean.wt.catch<- function (wk.dir, stock.list) 
{
    miny <- max(sapply(stock.list, function(x) x@range["minyear"]))
    maxy <- max(sapply(stock.list, function(x) x@range["maxyear"]))
    years <- miny:maxy
    new.list <- lapply(stock.list, trim, year = years)
    wt.in.catch <- lapply(new.list, function(x) {
        tot.wt <- seasonSums(quantSums(unitSums(x@catch.n * 
            x@catch.wt)))
        tot.num <- seasonSums(quantSums(unitSums(x@catch.n)))
        return(round(1000 * tot.wt/tot.num, 2))
    })
    tmp <- cbind(years, data.frame(matrix(unlist(wt.in.catch), 
        length(years), length(new.list))))
    colnames(tmp) <- c("Year", lapply(new.list, name))
    write.table(tmp, paste(wk.dir, "Mean weights in catch.csv", 
        sep = ""), row.names = FALSE, sep = ",")
    return(tmp)
    print("weights saved in file:")
    cat("\n")
    print(paste(Wkdir, "output/Mean weights in catch.csv", 
        sep = ""))
    cat("\n")
}

mean.wt.catch2<- function (stock.list) 
{
  miny <- max(sapply(stock.list, function(x) x@range["minyear"]))
  maxy <- max(sapply(stock.list, function(x) x@range["maxyear"]))
  years <- miny:maxy
  new.list <- lapply(stock.list, trim, year = years)
  wt.in.catch <- lapply(new.list, function(x) {
    tot.wt <- seasonSums(quantSums(unitSums(x@catch.n * 
                                              x@catch.wt)))
    tot.num <- seasonSums(quantSums(unitSums(x@catch.n)))
    return(round(1000 * tot.wt/tot.num, 2))
  })
  tmp <- cbind(years, data.frame(matrix(unlist(wt.in.catch), 
                                        length(years), length(new.list))))
  colnames(tmp) <- c("Year", lapply(new.list, name))

  return(tmp)
}