# # Epage # # Written by: # -- # John L. Weatherwax 2009-04-21 # # email: wax@alum.mit.edu # # Please send comments and especially bug reports to the # above email address. # #----- # version$language == "R" for R version$language == NULL for SPlus if(is.null(version$language) == FALSE) require(alr3) else library(alr3) # Oehlert's cake baking data #cakes <- read.table("data/cakes.txt", header=TRUE) if(is.null(version$language) == FALSE) data(cakes) attach(cakes) m1 <- lm(Y ~ X1 + X2 + I(X1^2) + I(X2^2) + X1:X2, data=cakes) s1 <- summary(m1) m2 <- lm(Y ~ X1 + X2 + I(X1^2) + I(X2^2) , data=cakes) s2 <- summary(m2) anova(m2,m1)