search for: testfit1

...2)/n) # 'n' is number of observations Residuals are computed by fitting a harmonic function on both the data: >testfit<-lm(co2obs~1+time+I(time^2)+sin(2*pi*time)+cos(2*pi*time)+sin(4*pi*time)+cos(4*pi*time)+sin(6*pi*time)+cos(6*pi*time)+sin(8*pi*time)+cos(8*pi*time),data=file) # >testfit1<-lm(co2model~1+time+I(time^2)+sin(2*pi*time)+cos(2*pi*time)+sin(4*pi*time)+cos(4*pi*time)+sin(6*pi*time)+cos(6*pi*time)+sin(8*pi*time)+cos(8*pi*time),data=file)# 'time' is time of observation testfit$residuals # observed.residuals # (saved in seperate file by write.table) testfit1$resid...