bstraub
calculates credibility premiums in the B�hlmann-Straub
credibility model.bstraub(ratios, weights,
heterogeneity = c("iterative", "unbiased"),
TOL = 1e-06, echo = FALSE)
"iterative"
for the
Bischel-Straub estimator; "unbiased"
for the usual
B�hlmann-Straub estimator (see below)heterogeneity =
"unbiased"
) of the between contracts heterogeneity parameter is
$$\hat{a} = c \left( \sum_{i = 1}^I w_{i\cdot} (X_{iw} - X_{ww})^2 -
(I - 1)\hat{s}^2 \right),$$
where $c = w_{\cdot\cdot}/(w_{\cdot\cdot}^2 - \sum_{i = 1}^I
w_{i\cdot}^2)$ and $I$ is the
number of contracts.
The Bishel-Straub pseudo-estimator (heterogeneity =
"iterative"
) is obtained recursively as the solution of
$$\hat{a} = \frac{1}{I - 1} \sum_{i=1}^I z_i (X_{iw} - X_{zw})^2.$$
The fixed point algorithm is used up, with a relative error of
TOL
stopping criteria.NA
in both the matrix of ratios and
the matrix of weights. The function can cope with complete lines
of NA
in case a contract has no experience.data(hachemeister)
## Credibility premiums calculated with the iterative estimator
bstraub(hachemeister$claims, hachemeister$weights)
## Credibility premiums calculated with the unbiased estimator
bstraub(hachemeister$claims, hachemeister$weights, heterogeneity = "unbiased")
Run the code above in your browser using DataLab