search for: omega2

...t 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*tim...

Dear R-helpers, I have time series data from 16 subjects and 2 treatment groups. The seasonal variation can be best described by two harmonics, I called the frequencies omega and omega2. I now want to test whether (1) the seasonal pattern differs between the treatments (some kind of overall test). If this is the case, (2) I want to conduct tests to find out whether the amplitude of the fist harmonic, the amplitude of the second harmonic, the phase of the first harmonic and the...

...eautiful diagram on p643 (figure 18.5). The gsl package does sn, cn, dn et seq, but before I put a huge amount of effort in, has anyone else coded up any related material? I'd be interested to hear if anyone has anything on Weierstrass's P function (in particular, computing omega1 and omega2 from general g2 and g3) and his sigma functions. Jacobi's theta functions would be good too. -- Robin Hankin Uncertainty Analyst Southampton Oceanography Centre European Way, Southampton SO14 3ZH, UK tel 023-8059-7743

...(airmiles)==13),Dec02=1*(seq(airmiles)==84)),method='ML') # Additive outliers are incorporated as dummy variables in xreg. # Transfer function components are incorporated by the xtransf and transfer # arguments. # Here, the transfer function consists of two parts omega0*P(t) and # omega1/(1-omega2*B)P(t) where the inputs of the two transfer # functions are identical and equals the dummy variable that is 1 at September # 2001 (the 69th data point) and zero otherwise. # xtransf is a matrix whose columns are the input variables. # transfer is a list consisting of the pair of (MA order, AR order...

- These recipients of your message have been processed by the mail server: renato-costa.vieira-pinto at serpro.gov.br; Failed; 5.3.0 (other or undefined mail system status) -------------- next part -------------- An embedded message was scrubbed... From: unknown sender Subject: no subject Date: no date Size: 143 Url: