library(MPV)
Delivery<-p8.3

x1<-Delivery$x1
x2<-Delivery$x2
y<-Delivery$y

#Model 3 X1+X2
fit.Del12<-lm(y~x1+x2,data=Delivery)
summary(fit.Del12)
anova(fit.Del12)

#Model 2 X1
fit.Del1<-lm(y~x1,data=Delivery)
summary(fit.Del1)
anova(fit.Del1)

#Model 1 X2
fit.Del2<-lm(y~x2,data=Delivery)
summary(fit.Del2)
anova(fit.Del2)

#Model 0 1
fit.Del<-lm(y~1,data=Delivery)
summary(fit.Del)
anova(fit.Del)


drop1(fit.Del12,test='F')

#Model 3 X1+X2
fit.Del12<-lm(y~x2+x1,data=Delivery)
summary(fit.Del12)
anova(fit.Del12)