p = bio.snowcrab::load.environment()
w = snowcrab.db( DS ="set.biologicals", p=p )
w = w[ polygon_inside( w, "cfanorth"),]
out = data.frame( cbind(
no.male = tapply( w$totno.male, w$yr, sum, na.rm=T) ,
no.female = tapply( w$totno.female, w$yr, sum, na.rm=T)
) )
out$year = rownames( out )
write.csv( out, file="cfanorth.counts.csv")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.