计算R中的PRESS值

Question

我有一个线性回归，我希望在R中计算它的PRESS值，我似乎无法找到允许我这样做的任何函数或包，我该怎么做？

由于

Answer 1

我找到了这3个解决方案：

y<-rnorm(20)
x<-rnorm(20)

mod<-lm(y~x)

qpcR::PRESS(mod) 

.........10.........20
$stat
[1] 22.87275

$residuals
 [1]  1.55408878 -0.40260735  0.06609265 -0.17136355 -0.66509310 -0.05861872  2.17710614 -0.12711591
 [9] -0.98517163 -1.30087693 -0.03429331  0.77811943  1.34654198 -0.68948005 -0.47502687  1.05951306
[17]  1.36271586  0.01317017 -0.59414220 -2.43794103

$P.square
[1] -0.1417693

MPV::PRESS(mod) # just the stat
[1] 22.87275

DAAG::press(mod) # idem
[1] 22.87275

希望其中一个能为你效劳

Answer 2

试试pls包，我用它得到了很好的结果： https://cran.r-project.org/web/packages/pls/