如何从 R 中的 princomp 对象获取系数?

机器算法验证 r 主成分分析
2022-03-22 08:21:37

基本上,我想手动计算predict(princomp(somedata))产生的结果。

为了做到这一点,我需要有系数矩阵,以便

newdata %*% coefficients

将产生每个组件的预测值。

任何人都可以给我一个线索来获得系数吗?

2个回答

如果您阅读有关该princomp函数的文档,?princomp您会看到有一个标题为“值”的部分描述了返回的对象的组件。

您正在寻找的部分可能是负载。

此外,如果您stats:::predict.princomp在控制台上键入,您可以在调用时准确地看到 R 正在做什么predict

predict在这种情况下,由 完成的特定计算是,返回的对象scale(newdata, object$center, object$scale) %*% object$loadings在哪里没有办法将其减少为单个矩阵乘法,因为在一般情况下必须重新缩放以匹配原始数据。objectprincompnewdata