R/proportion.legal.r

In jae0/snowcrab: Snow crab assessment suite using aegis* and associated projects and data streams.

  proportion.legal = function(odb, region, year) {
   # sex codes
    male = 0
    female = 1
    sex.unknown = 2


    out= NULL

# Remove CW's outside norms and remove production (pre-sorted) samples
    i = which( odb$sex==male & odb$prodcd_id=="0" & odb$cw > 50 & odb$cw < 170 & odb$fishyr==year )
    j = which( odb$sex==male & odb$prodcd_id=="0" & odb$cw >= 95 & (odb$durometer >= 68 | odb$shell >=2) )

    r = polygon_inside(x=odb, region=aegis.polygons::polygon_internal_code(region), planar=F)
    total.n = intersect (r, i)
    total.com = intersect(total.n, j)

    out = c( length(total.n), length(total.com), round( length(total.com)/length(total.n) * 100) )

    return(out)
  }