PopulationCov <- function(n)
{
a <- matrix(NA, n, n);
a[] <- phi^abs(row(a) - col(a));
for (i in 1:n)
a[i, ] <- a[i, ] - mean(a[i, baseline]);
for (j in 1:n)
a[ , j] <- a[ , j] - mean(a[baseline, j]);
invisible(a);
}
e <- eigen(PopulationCov(n));
lines(e$vectors[,1], col = 2, lwd = 2);
The National Academies | 500 Fifth St. N.W. | Washington, D.C. 20001
Copyright © National Academy of Sciences. All rights reserved.
Terms of Use and Privacy Statement
"
Appendix B: R Code for Figure 9-2 ."
Surface Temperature Reconstructions for the Last 2,000 Years . Washington, DC: The National Academies Press,
2006 .
Please select a format: