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