p<-runif(387)^2 #387 independent p-values, non-uniform intersection null violated
res<-BBSGoF(p)
par(mfrow=c(2,2))
plot(res) #Tarone test, within-block correlation, beta density (for automatic k),
#and decision plot (number of rejected nulls)
Run the code above in your browser using DataLab