# NOT RUN {
data(stumplift)
library(nlme)
modConstPow<-gls(Productivity~Machine+Machine*I((Diameter-70)^2),
data=stumplift,
weights=varPower(),
corr=corAR1(form=~Stump|Machine))
# }
Run the code above in your browser using DataLab