Pattern? [Bioanalytics]

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

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,675 registered users;
70 visitors (0 registered, 70 guests [including 12 identified bots]).
Forum time: 00:19 CEST (Europe/Vienna)

Every man gets a narrower and narrower field of knowledge
in which he must be an expert in order to compete with other people.
The specialist knows more and more about less and less
and finally knows everything about nothing.    Konrad Lorenz

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