Dear 'R' Users,

I have CO2 timeseries which I want to fit with fourth harmonic function and
then would like to compute residuals.
I tried with two option but not sure which one is correct, kindly any one
can help correcting me:

#
#I would like to fit a model with Fourth harmonic function as:

#m(t)=a+bt+ct^2+c1sin(omega1*t)+d1cos(omega1.t)+c2sin(omega2.t)+d2cos(omega2.t)+c3sin(omega3.t)+d3cos(omega3.t)+c4sin(omega4.t)+d4cos(omega4.t)
# In my data 'file' I have two variable 'co2obs' and 'time'
# SO above function in R will look like EITHER as:

testfit<-lm(co2obs~1+time+(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)
#RESIDUALS COMPUTING for above data fit
testfit$residuals

OR as:

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)
 #RESIDUALS COMPUTING for above data fit
testfit$residuals

NOTE: The 'third' term (ct^2) in both above is wrriten differently, rest
terms are same.

Kindly advise which one in above is correct, OR both are wrong ?

Great Thanks,

Regards,
Yogesh

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to