Pattern? [Bioanalytics]

posted by Helmut Homepage – Vienna, Austria, 2014-08-15 06:16 (3938 d 18:22 ago) – Posting: # 13381
Views: 10,348

Hi Ken,

❝ In a BE study involving 26 subjects […] the last 6 subjects had many points above ULOQ.


Sometimes concentrations in a few subjects are higher than anticipated, but this should occur at random. Why in the last six subjects? What makes them so special? If I would be an inspector that would ring my alarm bell. I would ask R how likely the observed pattern could occur by pure chance (the code takes a good while to complete; be patient – 28 minutes on my machine):

set.seed(123456)
n   <- 1e8                      # no. of simulations; at least 1 mio needed
x   <- c(rep(0, 20), rep(1, 6)) # 20 "normal" subjects and 6 ">ULOQ"
obs <- 0                        # counter for the "last 6" pattern
for(j in 1:n) { if(sum(tail(sample(x), 6)) == 6) obs <- obs + 1 }
cat("p =", signif(obs/n, 4), "\n")


IMHO p = 4.49·10–6 should trigger an internal investigation. Some ideas/questions:

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,674 registered users;
55 visitors (0 registered, 55 guests [including 17 identified bots]).
Forum time: 00:38 CEST (Europe/Vienna)

Do not worry about your difficulties in mathematics.
I can assure you mine are still greater.    Albert Einstein

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