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

posted by d_labes  – Berlin, Germany, 2017-11-15 13:15  – Posting: # 17982
Views: 3,911

(edited by d_labes on 2017-11-15 13:47)

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.043354
# 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.032048
# 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.001495


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:

Activity
 Mix view
Bioequivalence and Bioavailability Forum |  Admin contact
19,289 posts in 4,100 threads, 1,317 registered users;
online 11 (2 registered, 9 guests [including 6 identified bots]).

A refund for defective software might be nice,
except it would bankrupt the entire software industry
in the first year.    Andrew S. Tanenbaum

The BIOEQUIVALENCE / BIOAVAILABILITY FORUM is hosted by
BEBAC Ing. Helmut Schütz
HTML5