Learn R Programming

fishstat (version 2025.1.0.0)

fishstat-package: Global Fishery and Aquaculture Statistics

Description

The Food and Agriculture Organization of the United Nations (FAO) FishStat database is the leading source of global fishery and aquaculture statistics and provides unique information for sector analysis and monitoring.

This package provides the global production data from all fisheries and aquaculture in R format, ready for analysis.

Arguments

Author

Arni Magnusson and Rishi Sharma created this R package.

All credit for the FishStat database goes to the Statistics Team of the FAO Fisheries and Aquaculture Division, as well as national data submitters. The database terms of use are based on the CC BY-NC-SA 3.0 IGO license. The R package is released under a similar CC BY-NC-SA 4.0 license.

To cite the use of FishStat data:

FAO. [Year]. FishStat data. Fisheries and Aquaculture Division. Rome. https://www.fao.org/fishery/en/fishstat.

To cite the use of this R package to access the data, cite FishStat (above) as well as:

Magnusson, A. and R. Sharma. [Year]. fishstat: Global Fishery and Aquaculture Statistics. R package version [Version]. https://cran.r-project.org/package=fishstat.

Details

Production tables:

aquacultureaquaculture production
capturecapture production
productionaquaculture and capture production

Lookup tables:

areafishing areas
countrycountries and territories
environmentaquaculture environments
measureunits of measurement
sourcesources of production
speciestaxonomic groups
statusstatus of data entries

See Also

Examples

Run this code
head(production)

# Analyze production measured in tonnes
prod <- production[production$measure == "Q_tlw" & production$value > 0,]
prod <- merge(prod, species[c("species", "yearbook")])

# Select SOFIA species, excluding mammals, reptiles, and plants
prod <- prod[prod$yearbook ==
             "Aquatic animals (Fish, crustaceans and molluscs, etc.)",]

# Determine origin
prod <- merge(prod, area[c("area", "inlandmarine")])
prod$origin <- ifelse(prod$source == "CAPTURE", "Capture", "Aquaculture")
prod$w <- ifelse(prod$inlandmarine == "Marine areas", "marine", "inland")
prod$origin <- paste0(prod$origin, " (", prod$w, ")")

# World capture fisheries and aquaculture production
x <- xtabs(value~year+origin, aggregate(value~year+origin, prod, sum))
x <- x[,c(2,1,4,3)] / 1e6
library(areaplot)
areaplot(x, legend=TRUE, args.legend=list(x="topleft"), ylab="million tonnes")

Run the code above in your browser using DataLab