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