# NOT RUN {
library(dplyr)
df <- MTPL2 %>%
mutate_at(vars(area), as.factor) %>%
mutate_at(vars(area), ~biggest_reference(., exposure))
mod1 <- glm(nclaims ~ area + premium, offset = log(exposure), family = poisson(), data = df)
mod2 <- glm(nclaims ~ area, offset = log(exposure), family = poisson(), data = df)
x <- rating_factors(mod1, mod2, model_data = df, exposure = exposure)
autoplot(x)
# }
Run the code above in your browser using DataLab