p<-runif(387)^2 #387 p-values, intersection null violated
res<-BBSGoF(p)
summary(res) #automatic number of blocks, number of rejected nulls,
#estimated FDR, beta and beta-binomial parameters,
#Tarone test of no correlation
Run the code above in your browser using DataLab