library(TSA)
data(wages)
lm2_wages=lm(wages~time(wages)+I(time(wages)^2))
summary(lm2_wages)
xfit<-time(wages)
yfit_2<-fitted(lm2_wages)
win.graph(width = 4.8,height = 8,pointsize = 8)
plot(wages,ylab='Time value',xlab='Time')
points(wages,,col='blue')
lines(as.vector(xfit),as.vector(yfit_2),col='red')