Function to estimate the seasons

Given the fraction of the year when each of the indices are observed, this function estimates the number of seasons, the fraction of the year of each season and the season when each of the indices are observed.

periods(findicesB = NULL, findicesP = NULL)

Arguments

findicesB

A named vector with fraction of the year when each of the total biomass indices are observed.

findicesP

A named vector with fraction of the year when each of the proportion of recruits indices are observed.

Value

A list with four numeric elements: nper, with the number of seasons; f with the fraction of the year corresponding to each of the seasons; perindicesB with the season in which the index in total biomass is observed, this object is a named vector with one element per index; and perindicesP with the season in which the index of proportion of recruits is observed, this object is a named vector with one element per index.

See also

Author

Leire Ibaibarriaga & Sonia Sanchez.

Examples

# Load data
data(ane)

# Generate population estimates, given some estimated parameters
per <- periods( findicesB=unlist(lapply( indicesB.ane, function(x) mean(range(x)[c('startf','endf')]))), 
                findicesP=unlist(lapply( indicesP.ane, function(x) mean(range(x)[c('startf','endf')]))))
class(per)
#> [1] "list"

per$nper
#> [1] 2
per$f
#> [1] 0.375 0.625
per$perindicesB
#>     depm acoustic 
#>        2        2 
per$perindicesB
#>     depm acoustic 
#>        2        2