示例:
#loading data
data(iris)
#PCA
ir.pca <- prcomp(iris[,1:4],scale=T,center=F)
#date frame to check
dt <- data.frame(1,0,0,0)
names(dt) <- names(iris[,1:4])
#prediction how it is computed?
predict(ir.pca,dt)
#things of PCA
ir.pca$rotation
ir.pca$center
ir.pca$scale
如果中心为FALSE,那么方程是显而易见的:
predict(ir.pca,dt)[1] == ir.pca$rotation[1,1]/ir.pca$scale[1]
但是那个中心是什么?然后我不知道如何使用中心,旋转和比例来编写公式。