set.seed(42)
X <- matrix(rnorm(30 * 100), ncol = 30)
y <- rnorm(100)
ppls.cv <- penalized.pls.cv(X, y, lambda = c(1, 10), ncomp = 3)
myjack <- jack.ppls(ppls.cv)
Sigma <- vcov(myjack)
Sigma[1:5, 1:5]
Run the code above in your browser using DataLab