MackChainladder Example

import pandas as pd
import chainladder as cl

This example demonstrates how you can can use the Mack Chainladder method.

# Load the data
data = cl.load_sample('raa')

# Compute Mack Chainladder ultimates and Std Err using 'volume' average
mack = cl.MackChainladder()
dev = cl.Development(average='volume')
mack.fit(dev.fit_transform(data))

plot_data = mack.summary_.to_frame(origin_as_datetime=False)
plot_data

	Latest	IBNR	Ultimate	Mack Std Err
1981	18834.0	NaN	18834.000000	NaN
1982	16704.0	153.953917	16857.953917	0.761504
1983	23466.0	617.370924	24083.370924	356.573425
1984	27067.0	1636.142163	28703.142163	446.086164
1985	26180.0	2746.736343	28926.736343	771.404425
1986	15852.0	3649.103184	19501.103184	809.026792
1987	12314.0	5435.302590	17749.302590	830.803503
1988	13112.0	10907.192510	24019.192510	2047.154824
1989	5395.0	10649.984101	16044.984101	1917.116582
1990	2063.0	16339.442529	18402.442529	7288.933411

Show code cell source

Hide code cell source

import matplotlib.pyplot as plt
plt.style.use('ggplot')
%config InlineBackend.figure_format = 'retina'

ax = plot_data[['Latest', 'IBNR']].plot(
    kind='bar', stacked=True, ylim=(0, None), grid=True,
    yerr=pd.DataFrame({'latest': plot_data['Mack Std Err']*0,
                       'IBNR': plot_data['Mack Std Err']}),
    title='Mack Chainladder Ultimate',
    xlabel='Accident Year', ylabel='Loss');