HP.pri.start | R Documentation |
Estimates starting values for the priors of the modified 9-parameter Heligman-Pollard model from mortality rates by age and proportion of individuals death from non natural causes.
HP.pri.start(Si.qx, Prop.risk, Life.qx, age, rg = 0.25, control = nls.lm.control(maxiter = 100))
Si.qx |
A vector of the same length as "age" containing the total mortality rates by age (qx). |
Prop.risk |
A vector of the same length as "age" containing the proportion of death animals identified as death by external effects (other different than natural mortality). |
Life.qx |
Mortality vector qx form life table created with life.tab. |
age |
A vector containing the ages at which each age interval begins. |
rg |
The variation of the range at which the mean prior values will be multiplied to calculate low and high limits. |
control |
Allow the user to set some characteristics Levenberg-Marquardt nonlinear least squares algorithm implemented in nls.lm. see nls.lm.control |
A dataframe with low, high and mean starting values for priors, to be used for fitting a adapted 9-parametres Heligman-Pollard model
Heligman, L. and Pollard, J.H. (1980). The Age Pattern of Mortality. Journal of the Institute of Actuaries 107:49–80.
Moré, J.J. (1978). The Levenberg-Marquardt algorithm: implementation and theory, in: Lecture Notes in Mathematics 630: Numerical Analysis, G.A. Watson (Ed.), Springer-Verlag: Berlin, pp. 105-116.
HP.priors HP.mod nlsLM life.tab
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.