0

我有一个数据集，我从中进行 pls 然后进行预测。现在我需要评估预测数据是否适合使用 Hotelling 的 T2 和 Q-Residuals 的模型（如果您将 Hotelling 的 T2 放在 x 上，将 Q-Residuals 放在 y 轴上，您会得到一个所谓的“影响图”，这有助于评估样本数据 (plsdat$X1) 是否适合模型数据 (plsdat$X)。

我无法从预测中计算 Hotelling 的 T2 和 Q 残差，有人建议吗？

require(pls)
plsdat <- list()

plsdat$X <- matrix(rnorm(5000), ncol = 50)
plsdat$Y <-  rnorm(nrow(plsdat$X))

plsdat$plsmod <- plsr(plsdat$Y ~ plsdat$X
                      , ncomp = 4
                      , validation = "none"
                      , segments = 5
                      , segment.type="random"
                      , length.seg=2)

plsdat$X1 <- matrix(rnorm(5000), ncol = 50)

plsdat$prediction <- predict(object = plsdat$plsmod
                          , newdata = plsdat$X1
                          , ncomp = 4
                          , se.fit = T, interval = "confidence", level = 0.95)
4

0 回答 0