*
* Example 3.2.9 on page 99
*
open data sunspots.dat
calendar 1770
data(format=free,org=columns) 1770:1 1869:1 sunspots
*
@bjident(number=40) sunspots
*
* You can get the autocovariances by using the instruction CORRELATE with the
* option COVARIANCE. METHOD=YULE gives the same method of computing these as
* described in the book.
*
corr(covariance,number=4,method=yule) sunspots