Alternatively, the update() function can be used.

lm.fit1 <- update(lm.fit, ~. - age)

Try updating the following model by removing the rm variable and adding the lstat variable

Assume that:

• The MASS library has been loaded
• The Boston dataset has been loaded and attached