如何从主成分分析中计算指数/分数?
这是一个可重现的示例
set.seed(1)
dat <- data.frame(
Diet = sample(1:2),
Outcome1 = sample(1:10),
Outcome2 = sample(11:20),
Outcome3 = sample(21:30),
Response1 = sample(31:40),
Response2 = sample(41:50),
Response3 = sample(51:60)
)
ir.pca <- prcomp(dat[,3:5], center = TRUE, scale. = TRUE)
summary(ir.pca)
loadings <- ir.pca$rotation
scores <- ir.pca$x
correlations <- t(loadings)*ir.pca$sdev
这会产生三个主要成分。我可以使用这些来计算上述数据中每一行的分数或称为“响应指数”的指数吗?