Você não deu detalhes suficientes nem forneceu CMR. Eu fiz a mais ingenua das tentativas e deu certo para mim.

require(nlme)

m0 <- gls(dist~poly(speed, degree=2), data=cars)

summary(cars)
pred <- data.frame(speed=seq(4, 25, by=0.5))

predict(m0, newdata=cars)

À disposição.

Walmes.