Caros colegas,

não estou conseguindo estimar os betas dos modelos quadrático e cúbico da regressão através da função reg.poly() do pacote ExpDes.    

x<-c(1000000, 10000000, 50000000, 100000000)

Glresíduo<- 72

Sqresíduo<- 1159.9

Gltratamento<- 3

Sqtratamento<- 167.6

reg.poly(y,x,72,1159.9,3,167.6)

Agradeço pela atenção

Att.

Tiago.