Programming the statistical procedures from SAS

Autoregressive model with random times

Frequent Contributor
Posts: 75

Autoregressive model with random times

Hello all,

I am using type=sp(pow)(time) for an AR(1) model with random times and I am getting two parameter estimates: the residual variance s^2 and correlation coef p

For my example s^2=8.8771 and p=.04977

Now my R matrix is a 2x2 matrix [{8.8771, 3.4493},{3.4493,8.8771}]

I can see why the diagonal elements which are the variances are 8.8771 but I dont understand where the 3.4493 is coming from. Is it the covariance between to adjacent observations?

Is cov(i,j)=s^2*p^|t(i)-t(j)| ????

Respected Advisor
Posts: 2,655

Re: Autoregressive model with random times

According to the PROC MIXED documentation,

cov(i,j) = s^2*p^(d(i,j)) where d(i,j) is the Euclidean distance between i and j

d(i,j) = sqrt( sum [over m from 1 to k](c(m,i) - c(m,j)^2).

So for one dimensional data, like time series, your expression should be correct.

Steve Denham
Ask a Question
Discussion stats
  • 1 reply
  • 2 in conversation