sp.index
Sparial indices
Description
Estimation of abundance, biomass, mean individual weight (MIW), and sex ratio indices over a statistical grid
Usage
sp.index(TA, TB, sspp, type, range=NA, GSA=NA, country=NA, threshold=NA, grid.sf=cgpmgrid,land=countries)
Arguments
TA |
data frame containing the hauls data (TA, table A). |
TB |
data frame containing the catches data (TB, table B). |
sspp |
MEDITS code of the reference species |
type |
type of index to estimate: "abundance", "invCV", "biomass", "MIW", "sex ratio" |
range |
coordinate range of the study area in the form of c(xmin,xmax, ymin,ymax) |
GSA |
integer value corresponding to the GSA number |
country |
string value indicating the selected country for the analysis in case the analysis should be performed by country |
threshold |
numeric value indicating the minimum number of specimens per haul to be considered in the sex ratio estimation. If the threshold parameter is not defined all the available specimens are considered in the estimation of the sex ratio |
grid.sf |
... |
land |
... |
Details
Mean biomass and mean abundance
The mean abundance (likewise the mean biomass) in the GFCM grid cells (Dcell) is calculated as the average of the standardized numbers of individuals (number/km2) over the most recent 10 years of the time series (if the time series is shorter than 10 years, all the available year data are considered):
where n is the count of the combinations year-haul in all the last 10 years. The variance of the mean
abundance in the cells is calculated as:
The CV is calculated as the ratio between the standard deviation of the mean annual value by haul and year (numerator) and the mean biomass (or abundance) in the cell (denominator).
For each GFCM cell the mean individual weight is calculated by year y and haul h as ratio between the total weight W in the haul and the total number N in the haul (from the MEDITS samples data) as follows:
where n is the count of the combinations year-haul in all the last 10 years.
The variance of the MIW in the cells is calculated using the following formula:
Then, the Coefficient of Variation (CV) of the MIW is calculated as:
Sex ratio
The sex ratio in each GFCM cell is calculated as the ratio between the sum of the standardized number of females and the sum of the standardized number of males and females over the hauls of the last 10 years:
where NF and NM are the standardized number of the females and of males in the haul h and year y.
The variance of the sex ratio in the cell is calculated using the following formula:
where n is the count of the combinations year-haul in all the last pooled 10 years.
The CV is calculated as the ratio between the standard deviation of the sex ratio by haul and year to
the sex ratio in the cell.
Inverse of mean abundance Coefficient of Variation (CV)
The inverse of the coefficient of variation of the mean abundance by GFCM cell is plotted.
Value
the function returns the plot of the selected indices over the GFCM (General Fisheries Commission for the Mediterranean) statistical grid. Moreover, the data frame containing the values of the estimated indices is returned
cgpmgridlevel |
identification number value of the corresponding grid cell. |
GSA |
the corresponding geographical sub-area (GSA) of the relative cell. |
meanNkm2 |
mean value of the indices. |
sdNkm2 |
standard deviation of the mean value of the indices. |
cvNkm2 |
Coefficient of variation of the mean indices. |
inverse_cvNkm2 |
Inverse value of the mean indices coefficient of variation. |
nhauls |
number of hauls in any given cell. |
positive_hauls |
number of positive hauls used for the estimation of the indices. |
lon |
longitude coordinate of the cell's centroid |
lat |
latitude coordinate of the cell's centroid |
Author(s)
Walter Zupa
Examples
library(MEDITS) sp.index(TA,TB,sspp="ARISFOL",type="abundance",range = c(10, 20, 38, 42))