GMR restrictions? [Conferences / Workshops]

posted by mittyri – Russia, 2016-09-21 02:32 (2745 d 17:04 ago) – Posting: # 16654
Views: 5,737

(edited by mittyri on 2016-09-21 13:41)

Hi Helmut!

could you please explain this:

❝ Statistically the GMR-restriction makes no sense indeed but I guess it will stay for a long time due to political reasons.


My little undestanding:
library(PowerTOST)
library(ggplot2)
rs <- reg_const("EMA")
rs$pe_constr <- FALSE
GMR <- 0.7999
y1 <- data.frame(
  N=integer(), Power=double(), LowerCI=double(), UpperCI=double()
  )

n1 <- 80:150
for(i in seq_along(n1))
{
 
  y1[nrow(y1)+1,] <-c(N=n1[i],
                      Power = power.scABEL(CV=0.5, n=n1[i], theta0=GMR, regulator=rs),
                      LowerCI=CI.BE(pe=GMR, CV=0.5, n=n1[i], design = "2x2x4")[1],
                      UpperCI=CI.BE(pe=GMR, CV=0.5, n=n1[i], design = "2x2x4")[2]   )
}
ggplot(y1)+
  geom_point(aes(N, Power, colour = Power))+
  geom_line(aes(N, LowerCI ), colour = "red")+
  geom_line(aes(N, UpperCI), colour = "red")+
  geom_hline(yintercept = 0.6984, color="red", size=3)+
  ylab("The Power and CIs")+
  ggtitle("CV = 50%, GMR = 79.99%")

[image]

So not necessary that (GMR is inside 0.8-1.25)&(CIs are inside the limits)
May be I understtod something wrong, please feel free to correct me and criticize my code :-D

Kind regards,
Mittyri

Complete thread:

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

Nothing shows a lack of mathematical education more
than an overly precise calculation.    Carl Friedrich Gauß

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