powered by
head(cakes) lm(Y~X1+X2+I(X1^2)+I(X2^2)+X1:X2, data=cakes)
Run the code above in your browser using DataLab