I want to be able to index my model after having fit the model. Say I have
df <- data.frame(a = c(1,2,3), b = c(2,3,1000),