# NOT RUN {
model <- lm(Sepal.Length~., data=iris)
new_observation <- iris[1,]
br <- broken(model, new_observation)
plot(br)
# works for interactions as well
model <- lm(Sepal.Length ~ Petal.Width*Species, data = iris)
summary(model)
new_observation <- iris[1,]
br <- broken(model, new_observation)
br
plot(br)
br2 <- broken(model, new_observation, predict.function = betas)
br2
plot(br2)
# }
Run the code above in your browser using DataCamp Workspace