PowerTOST CV Power EMA 3 periods discontinuity [🇷 for BE/BA]

posted by mittyri – Russia, 2016-02-01 23:33 (2978 d 06:37 ago) – Posting: # 15912
Views: 5,571

Dear All,

Could you please enlighten me about the following issue. May be it was discussed somewhere but I didn't find.
Some R code:
x1 <- seq(0.25, 0.5, length.out=100)
y1 <- vector("numeric", length=length(x1))
y2 <- vector("numeric", length=length(x1))
des = c("2x3x3", "2x2x3")
n<-48
GMR <- 0.95
for (j in seq_along(x1)) {
  if (x1[j] > 0.3) {
    y1[j] <- power.scABEL(CV=x1[j], theta0=GMR, n=n, des=des[1], reg="EMA")
    y2[j] <- power.scABEL(CV=x1[j], theta0=GMR, n=n, des=des[2], reg="EMA")
  }
  if (x1[j] <= 0.3) {
    y1[j] <- power.TOST(CV=x1[j], theta0=GMR, n=n, des= des[1])
    y2[j] <- power.TOST(CV=x1[j], theta0=GMR, n=n, des= des[2])
  }
}
plot(x1, y2,type="p",col="red", xlab="CV", ylab="Power")
lines(x1, y1,col="green")
legend('topright', c("TRT RTR", "TRR RRT RTR"), lty =1, col = c('red', 'green'))

The resulted plot:
[image]
There is a discontinuity at CV=0.3. I cannot figure out the reason. The limits are widening smoothly from this point.

PS: some intersection with Shuanghe's post :cool:


Edit: Category changed. [Helmut]

Kind regards,
Mittyri

Complete thread:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
91 visitors (0 registered, 91 guests [including 5 identified bots]).
Forum time: 06:10 CET (Europe/Vienna)

With four parameters I can fit an elephant,
and with five I can make him wiggle his trunk.    John von Neumann

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