Problem Statement:
Under the supposition that for j > 2,
(a) Obtain approximate standard errors for r1, r2, and rj, j > 2.
(b) Obtain the approximate correlation between r4 and r5.
Solution:
For part (a), we apply (2.1.11), p. 32:
For k = 1,
> s:=(1+2*r1^2)*(1+2*r1^2+2*r2^2)+(r1^2+2*r2)-4*r1*(2*r1*r2+2*r1); 2 2 2 2 s := (1 + 2 r1 ) (1 + 2 r1 + 2 r2 ) + r1 + 2 r2 - 4 r1 (2 r1 r2 + 2 r1) > simplify(s); 2 2 4 2 2 2 1 - 3 r1 + 2 r2 + 4 r1 + 4 r1 r2 + 2 r2 - 8 r1 r2
Now, for k=2,
> t:=(1+2*r1^2)*(1+2*r1^2+2*r2^2)+r2^2-4*r2*(2*r2+r1^2); 2 2 2 2 2 t := (1 + 2 r1 ) (1 + 2 r1 + 2 r2 ) + r2 - 4 r2 (r1 + 2 r2) > simplify(t); 2 2 4 2 2 2 1 + 4 r1 - 5 r2 + 4 r1 + 4 r1 r2 - 4 r1 r2
Finally, the last part of (a) is k > 2. By supposition,
, so we can immediately drop the terms
and
from (11). Also, for any integer v, at least
one of v+k or v-k is outside the interval [-2,2], so
the term
is also . Thus,
Finally, for part (b), we use (2.1.14), p. 34. Since
4 and 5 are bigger than 2, this formula applies (note that
it only applies once we are into the region where k > q
and q is the lag of the last nonzero autocorrelation).
We have
Remarks: The takehome message is that these formulae exist. If you are like me and aren't satisfied until you know where they come from, you can take our advanced time series class. The chances of these formulae ever coming up on an exam are two: slim and none.
Remember: these results are valid for Gaussian processes. Under a ``long range'' independence assumption, one can derive formulae for the approximate variances (variances of asymptotic normal distributions) which depend on the first four moments, basically. The particular algebraic forms are highly dependent on the moments of multivariate normal distributions.