Julia, vc pode usar, por exemplo:
m1<-lm(Y~-1+cbind(x1,x2,x3,I(x1^2) ,I(x2^2),I(x3^2)
))
Luzia
Em 18/04/2012 17:47, Julia Maria escreveu:
m1<-lm(Y~-1+x1+x2+x3+I(x1^2) +I(x2^2)+I(x3^2)