head(area)
# Inland waters and marine areas
area[area$inlandmarine == "Inland waters", c("area", "area_name")]
area[area$inlandmarine == "Marine areas", c("area", "area_name")]
# Check if any area has zero production
nonzero <- unique(production$area[production$value > 0])
print(area[!(area$area %in% nonzero),], row.names=FALSE)
# Check which species groups are recorded in areas 98 and 99
species_98_99 <- unique(production$species[production$area %in% 98:99])
cbind(unique(species$isscaap[species$species %in% species_98_99]))
# Marine fishing areas in northern, central, and southern regions
area$area[area$faregion == "Northern regions"]
area$area[area$faregion == "Central regions"]
area$area[area$faregion == "Southern regions"]
# Examine one area
print.simple.list(area[area$area == 71,])
Run the code above in your browser using DataLab