function(k, n0, pstar) { # function to compute comparison with a standard critical values # for case D (mu0 unknown, CRN, no known structure) using the # tighter Bonferroni approximation # # k = number of systems, excluding standard # n0 = first-stage sample size # pstar = overall desired confidence level # # last update 4/23/97 # alpha <- 1 - pstar check <- function(x, k, n0, alpha, h) { abs(((1 - pt(x - h, n0 - 1)) + (k - 1) * (1 - pt(x, n0 - 1))) - alpha) } h <- qt(1 - alpha/k, n0 - 1) result <- optimize(check, interval = c(h, 2 * h), keep.xy = F, k = k, n0 = n0, alpha = alpha, h = h) g <- result$minimum hg <- matrix(c(h, g), ncol = 1) hg }