Displaying 1 result from an estimated 1 matches for "eigvnocent".
Did you mean:
eigvcent
2009 Mar 08
2
prcomp(X,center=F) ??
...not really a PCA.
See the following example, in which the second column of
data matrix X is linearly correlated to the first column:
> X <- cbind(rnorm(100,100,50),rnorm(100,100,50))
> X[,2] <- X[,1]*1.5-50 +runif(100,-70,70)
> plot(X)
> cor(X[,1],X[,2])
[1] 0.903597
> eigvnocent <- prcomp(X,center=F,scaling=F)[[1]]
> eigvcent <- prcomp(X,center=T,scaling=F)[[1]]
> eigvecnocent <- prcomp(X,center=F,scaling=F)[[2]]
> eigveccent <- prcomp(X,center=T,scaling=F)[[2]]
> PCnocent <- X%*%eigvecnocent
> PCcent <- X%*%eigveccent
> par(mfrow...