GMR restrictions? [Conferences / Workshops]

posted by mittyri – Russia, 2016-09-21 02:32 (2746 d 06:11 ago) – Posting: # 16654
Views: 5,743

(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;
85 visitors (0 registered, 85 guests [including 6 identified bots]).
Forum time: 07:43 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