suppressWarnings(RNGversion("3.5.0"))
set.seed(120)
x <- matrix(rnorm(1000*20),ncol=20)
y <- sample(c(1:2),size=20,replace=TRUE)
x[1:50,y==2]=x[1:50,y==2]+3
mydata <- list(x=x,y=y)
mytrain <- pamr.train(mydata)
myfdr <- pamr.fdr(mytrain, mydata)
pamr.plotfdr(myfdr)
Run the code above in your browser using DataLab