#install.packages("astsa") require(astsa) #help(package = astsa) plot(sunspotz) win.graph() acf(sunspotz, 100) ####################################################### # Figure 14.2: win.graph() seq.n = c(10,50,500) par(mfrow=c(3,1)) sigmasqd = 1/log(2+1:max(seq.n)) for (n in seq.n) { N = floor(n/2) lambda = 2*pi*(1:N)/n lambda = c(rev(-lambda), 0, lambda) lambda sig = sigmasqd[1:N] sig = c(rev(sig)/2, 1/log(2), sig/2) prob = sig/sum(sig) barplot(prob, xlab = paste("N=", N)) }