suppressWarnings(RNGversion("3.5.0"))
set.seed(120)
x <- matrix(rnorm(1000*20),ncol=20)
y <- sample(c(1:4),size=20,replace=TRUE)
mydata <- list(x=x,y=y)
mytrain <- pamr.train(mydata)
mycv <- pamr.cv(mytrain, mydata)
pamr.plotcv(mycv)
Run the code above in your browser using DataLab