1

我想使用函数使用 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也没有什么大的不同。

4

0 回答 0