Bioequivalence and Bioavailability Forum

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

posted by mittyri  – Russia, 2016-02-01 23:33 (3426 d 12:44 ago) – Posting: # 15912
Views: 6,531

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:

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

---

Edit: Category changed. [Helmut]

Complete thread:

• PowerTOST CV Power EMA 3 periods discontinuitymittyri 2016-02-01 22:33 
  • Discontinuity = artifact Helmut 2016-02-02 00:23
  • 2x Laszlo ElMaestro 2016-02-02 16:28
    • 2x Laszlo Helmut 2016-02-02 16:39
      • 2x Laszlo ElMaestro 2016-02-02 17:09