Bayesian Optimal Discovery Procedure (BODP) R code

The r code given as an additional file to the paper Bayesian optimal discovery procedure for simultaneous significance testing requires your data to be in matrix form. I spent a day trying to troubleshoot the code when all you need to do is the following:

Open R

Copy and paste the R code from Cao et al.

x<-read.table("x.txt") #your control state

y<-read.table("y.txt") #your test state

x<-as.matrix(x)

y<-as.matrix(y)

BODP(x,y,1000,10) #feel free to change 1000 and 10 to whatever you'd like