n <- 50 # number of observations
p <- 5 # number of variables
X <- matrix(rnorm(n * p), ncol = p)
y <- rnorm(n)
pls.object <- penalized.pls.cv(X, y)
my.jack <- jack.ppls(pls.object)
my.coef <- coef(my.jack)
print(my.coef)
Run the code above in your browser using DataLab