Hi, It may be my misunderstanding, but it seems that the "na.action" in the princomp() function for principal components analysis does not work. Please see this simple example:
u <- matrix(rnorm(75), ncol=1) v <- matrix(rnorm(20), ncol=1) x <- u%*%t(v) + matrix(rnorm(20*75),ncol=20) x[1,1] <- NA pc.out <- princomp(x, na.action=na.exclude) Error in cov.wt(z) : 'x' must contain finite values only > Note, I have: > options("na.action") $na.action [1] "na.omit" Thanks, Ravi [[alternative HTML version deleted]] ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel