Bioequivalence and Bioavailability Forum • Five minutes gone

Five minutes gone - power.tsd.2m() arose [Two-Stage / GS Designs]

posted by d_labes – Berlin, Germany, 2017-11-15 14:15 (2644 d 08:51 ago) – Posting: # 17982
Views: 9,932

Dear Helmut, dear All,

❝ ❝ The adjusted α showed (well, in simulations…) that the TIE is controlled if we follow the frameworks exactly, i.e., initiate the second stage only if necessary based on the interim and with the re-estimated sample size. If we increase the sample size, the chance of passing BE increases and hence, the TIE.

❝ Wait, wait ... I will be back in 5 min :cool: with some quick shoot code to verify your claim.

5 minutes are gone :PCchaos:

.
See Power2Stage on Github. Function power.tsd.2m().
It's pure Potvin B, abbreviated with implicit power monitoring, i.e. via sample size estimation, which gives n2=0 in case of 'enough' power.
What's missing at the moment is some sort of correlation between the PK metrics analogous to what is possible in power.2TOST() with the argument rho. Here I don't know at the moment at which place it comes into play. Suggestions are welcome.
The function is not public at the moment. Not in the mood to document :sleeping:

.

Example:
# install from GitHub via
# devtools::install_github("Detlew/Power2Stage")

library(Power2Stage)

# Cmax not BE, but AUC assumed BE

Power2Stage:::power.tsd.2m(CV=c(0.3, 0.2), theta0=c(1.25, 0.95), n1=28)$pBE

[1] 0.043414

# Cmax assumed BE, but AUC not BE

Power2Stage:::power.tsd.2m(CV=c(0.3, 0.2), theta0=c(0.95, 1.25), n1=28)$pBE

[1] 0.031875

# Cmax not BE, AUC also not BE

Power2Stage:::power.tsd.2m(CV=c(0.3, 0.2), theta0=c(1.25, 1.25), n1=28)$pBE

[1] 0.001452

Take the results with a grain of salt. Don't know how to validate. Suggestions here also welcome.
But I think the results look plausible.

Edited numbers after removing a little bug.

—
Regards,

Detlew

Complete thread:

Two PK metrics: Inflation of the Type I Error Helmut 2017-11-12 11:57 [Two-Stage / GS Designs]
- Two PK metrics: Inflation of the Type I Error? d_labes 2017-11-12 17:46
- A place to start ElMaestro 2017-11-12 21:43
  - A place to start? d_labes 2017-11-13 16:02
    - A place to start? Helmut 2017-11-13 16:51
      - Scientific gut feeling d_labes 2017-11-15 10:34
        
        Five minutes gone - power.tsd.2m() arosed_labes 2017-11-15 13:15
    - A better place to start. ElMaestro 2017-11-14 00:04
      - Nope Helmut 2017-11-14 00:21
        
        Nope nobody 2017-11-14 08:11
        I've meditated hard ElMaestro 2017-11-14 12:41