helper functions to do p value adjustment
gof_adj_C1(dta, rnull, vals, TS, typeTS, TSextra, B = 1000L)
A matrix of powers
data set
R function (generate data under null hypothesis)
values of discrete random variable
function to calculate test statistics
integer indicating type of test statistic
list to pass to TS
=1000 Number of simulation runs