Rate a book against a market

from decimal import Decimal

from relativities import Market, Program

current = Program.load("effective://example-ca-ppa")
book = current.example_book
rates = current.rater.tables["BaseRate"]
collision = rates["coverage"] == "collision"
rates.loc[collision, "rate"] = rates.loc[collision, "rate"] * Decimal("1.05")
market = Market({"current": current, "proposed": current.with_rater(current.rater.with_tables({"BaseRate": rates}))})

results = market.rate(book)
results[["submission_id", "program", "status", "premium"]].head(4)
#   submission_id   program     status premium
# 0       S000001   current  succeeded    2636
# 1       S000001  proposed  succeeded    2694
# 2       S000002   current  succeeded    5339
# 3       S000002  proposed  succeeded    5461

premiums = results.pivot(index="submission_id", columns="program", values="premium").astype(float)
change = premiums["proposed"] / premiums["current"] - 1
change.describe()[["count", "mean", "min", "max"]].round(4)
# count    96671.0000
# mean         0.0163
# min          0.0000
# max          0.0355

matched = premiums.dropna()
round(float(matched["proposed"].sum() / matched["current"].sum() - 1), 4)  # 0.0181

results.query("status != 'succeeded'").groupby(["status", "program"]).size().unstack()
# program         current  proposed
# status
# invalid_inputs      394       394
# refused            2935      2935

steps = market.explain(book, "S000001").query("kind == 'calculation'")
steps.pivot(index="name", columns="program", values="value").query("current != proposed")
# program                                current proposed
# name
# base_rate.collision                       1519  1594.95
# premium                                   2636     2694
# vehicles[0].coverage_premium.collision     562      590
# vehicles[0].premium                       1265     1293
# vehicles[1].coverage_premium.collision     611      641
# vehicles[1].premium                       1371     1401

market.close()