calculates other sources of natural mortality other than predation. It is assumed that M1 mortality follows a beta distribution and that in general the small and large size classses suffer the greatest mortality.
calc_M1(nSize, nSpecies, lBound, mBound, alphaM1, betaM1, cM1, scLinfMat,
scLinf, phiMin, parameterValues)Number of size class intervals species can grow through
Number of species in the model
Lower bound of each size class interval
Mid point of each size class interval
Shape parameter (alpha) of the beta distribution
Shape parameter (beta) of the beta distribution
Scales the pdf ??????
The size class at which each species reaches L_inf (maximum length)
Model timestep (years). See calc_phi
Matrix of species specific parameters. See rochet_GB_parameterValues
A matrix is returned
A matrix (size nSize x nSpecies) of M1 (natural mortality) values. Note: M1 = 0 for size classes in which a species does not reach/grow in to.
Hall et al. (2006). A length-based multispecies model for evaluating community responses to fishing. Can. J. Fish. Aquat. Sci. 63:1344-1359.
Rochet et al. (2011). Does selective fishing conserve community biodiversity? Prediction from a length-based multispecies model. Can. J. Fish. Aquat. Sci. 68:469-486