Fit_References | R Documentation |
The function fits a collective model to a wishlist of references (expected layer losses and excess frequencies). The function allows to specify the family of the severity distribution that is used. Depending on this distribution family the function works slightly differently:
For the severity distribution PiecewisePareto
the function returns a PPP_Model
that satisfies all the references
For the severity distribution Pareto
the function returns a PPP_Model
that minimizes the squared relative deviations from the references
For the severity distribution GenPareto
the function returns a PGP_Model
that minimizes the squared relative deviations from the references
Fit_References(
Covers = NULL,
Attachment_Points = NULL,
Expected_Layer_Losses = NULL,
Thresholds = NULL,
Frequencies = NULL,
model_threshold = min(c(Attachment_Points, Thresholds)),
default_alpha = 2,
dispersion = 1,
alpha_max = 100,
severity_distribution = "PiecewisePareto",
ignore_inconsistent_references = FALSE
)
Covers |
Numeric vector. Vector containing the covers of the layers from the wishlist. |
Attachment_Points |
Numeric vector. Vector containing the attachment points of the layers from the wishlist. |
Expected_Layer_Losses |
Numeric vector. Vector containing the expected losses of the layers from the wishlist. |
Thresholds |
Numeric vector. Contains the thresholds from the wishlist for which excess frequencies are given. |
Frequencies |
Numeric vector. Expected frequencies excess the |
model_threshold |
Numerical. Lowest threshold of the fitted piecewise Pareto distribution. |
default_alpha |
Numerical. Default alpha for situations where an alpha has to be selected. |
dispersion |
Numerical. Dispersion of the claim count distribution in the resulting PPP_Model. |
alpha_max |
Numerical. Maximum alpha to be used for the matching. |
severity_distribution |
Character. Implemented distributions: "PiecewisePareto" (default), "Pareto" and "GenPareto". |
ignore_inconsistent_references |
Logical. If TRUE then inconsistent references are ignored in case of the piecewise Pareto distribution and the other references are used to fit the model |
For severity_distribution = "PiecewisePareto" or "Pareto": A PPP_Model object that contains the information about a collective model with a Panjer distributed claim count and a Piecewise Pareto distributed severity. The object contains the following elements:
FQ
Numerical. Frequency in excess of the lowest threshold of the piecewise Pareto distribution
t
Numeric vector. Vector containing the thresholds for the piecewise Pareto distribution
alpha
Numeric vector. Vector containing the Pareto alphas of the piecewise Pareto distribution
truncation
Numerical. If truncation
is not NULL
and truncation > max(t)
, then the distribution is truncated at truncation
.
truncation_type
Character. If truncation_type = "wd"
then the whole distribution is truncated. If truncation_type = "lp"
then a truncated Pareto is used for the last piece.
dispersion
Numerical. Dispersion of the Panjer distribution (i.e. variance to mean ratio).
Status
Numerical indicator: 0 = success, 1 = some information has been ignored, 2 = no solution found
Comment
Character. Information on whether the fit was successful
For severity_distribution = "GenPareto": A PGP_Model object that contains the information about a collective model with a Panjer distributed claim count and a Piecewise Pareto distributed severity. The object contains the following elements:
FQ
Expected claim count of the collective model.
t
Numeric. Threshold of the Pareto distribution.
alpha_ini
Numeric. Initial Pareto alpha (at t
).
alpha_tail
Numeric. Tail Pareto alpha.
truncation
If truncation
is not NULL
and truncation > t
, then the Pareto distribution is truncated at truncation
.
dispersion
Numerical. Dispersion of the Panjer distribution (i.e. variance to mean ratio).
Status
Numerical indicator: 0 = success, 1 = some information has been ignored, 2 = no solution found
Comment
Character. Information on whether the fit was successful
covers <- c(1000, 1000, 1000)
att_points <- c(1000, 2000, 5000)
exp_losses <- c(100, 50, 10)
thresholds <- c(4000, 10000)
fqs <- c(0.04, 0.005)
fit <- Fit_References(covers, att_points, exp_losses, thresholds, fqs)
Layer_Mean(fit, covers, att_points)
Excess_Frequency(fit, thresholds)
fit <- Fit_References(covers, att_points, exp_losses, thresholds, fqs,
severity_distribution = "Pareto")
Layer_Mean(fit, covers, att_points)
Excess_Frequency(fit, thresholds)
fit <- Fit_References(covers, att_points, exp_losses,
severity_distribution = "GenPareto")
Layer_Mean(fit, covers, att_points)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.