我想使用函数使用 PLSR 的载荷(权重)计算一组变量的分量plsr
。
我认为这些组件是通过将每个变量的值乘以估计的负载(重量)来计算的。
但是,使用输出plsr
并这样做并没有给我预期的值:
例子:
library("pls")
data(oliveoil)
sens.pcr <- plsr(sensory ~ chemical, ncomp = 4, scale = F, data = oliveoil)
提取载荷/重量:
df <- cbind(sens.pcr$loadings[,1],sens.pcr$loadings[,2],sens.pcr$loadings[,3],sens.pcr$loadings[,4])
一项测试观察:
firstrow <- oliveoil$chemical[1,]
提取组件(分数):
scores <- sens.pcr$scores
做线性组合:
sum(firstrow*df[,1])
[1] -12.81924
哪个不等于第一个分数scores[1,1]
=0.5100166
我错过了什么?
使用sense.pcr$loadings.weigths
也没有什么大的不同。