Technicality: Weigths for the inverse normal approach [Two-Stage / GS Designs]

posted by d_labes  – Berlin, Germany, 2018-04-25 16:19 (1979 d 19:46 ago) – Posting: # 18729
Views: 15,491

Dear Helmut,

great post:clap:.

Only one remark about the weights you choose for the maximum combination test in your R code.


❝     for (j in seq_along(CV)) {
❝       n[j] <- sampleN.TOST(alpha=alpha, CV=CV[j], theta0=GMR, theta1=theta1,
❝                            theta2=theta2, targetpower=targetpower,
❝                            print=FALSE, details=FALSE)[["Sample size"]]
❝       if (n[j] < 12) n[j] <- 12
❝       for (k in seq_along(n1)) {
❝         # median of expected total sample size as a 'best guess'
❝         n.tot <-, CV=CV[j], n1=n1[k], GMR=GMR,
❝                               usePE=usePE, theta1=theta1, theta2=theta2,
❝                               targetpower=targetpower, fCrit=fCrit,
❝                               fClower=fClower, fCNmax=fCNmax, pmethod=pmethod,
❝                               npct=0.5)$nperc[["50%"]]
❝         w     <- c(n1[k], n.tot - n1[k]) / n.tot
❝         # force extreme weights if expected to stop in stage 1 with n1
❝         if (w[1] == 1) w <- w + c(-1, +1) * 1e-6
❝     ...

Defining the weights that way is IMHO not what you intended. Or I don't understand what you intended.
It is correct if you think in terms of the standard combination test and think further that you have to specify two weights for that. But since the two weights are connected by w, 1-w the second is calculated within the function automatically. You only need to define w[1] in the input argument.

The idea behind the maximum combination test now is:
If our first pair of weights w, 1-w (chosen anyhow) is not "optimal", choose a second pair of weights w*, 1-w* which is more adapted to the real situation.
If you were too optimistic in your planing of n2, i.e. have chosen n2 too low compared to what really happens in the sample size adaption, it would be wise to define w* lower than w.
You do that, but your choice (w in w[1]=0.999999, w* in w[2]=1e-6) is too extreme I think and not your intention I further think. The second pair of weights w*=1e-6, 1-w*=0.999999 here is for a situation were the p-values from the second stage nearly exclusively determine the overall outcome of the maximum combination test. The p-values from the first stage data are down-weighted with w*=1e-6.

Hope this sermon is not too confusing.

BTW: Choosing the weights "optimal" is for me a mystery. To do that, we had to know the outcomes of the two stages, but we don't have them until the study has been done. On the other hand we have to predefine them to gain strict TIE control. Hier beißt sich die Katze in den Schwanz.



Complete thread:

UA Flag
 Admin contact
22,759 posts in 4,775 threads, 1,628 registered users;
19 visitors (0 registered, 19 guests [including 8 identified bots]).
Forum time: 12:06 CEST (Europe/Vienna)

Whenever a theory appears to you as the only possible one,
take this as a sign that you have neither understood the theory
nor the problem which it was intended to solve.    Karl R. Popper

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz