Can chance make you a killer?
bbc.co.uk
bbc.co.uk
If you have chance 1/100 of death then with sample 100 you have to have 0 or 1 deaths to fit in 60% but with sample 200 you can have 0,1,2,3 deaths to still fit.
It's just 60% is harsh limit for small samples. There should be some lower limit on sample size you have to gather for given probability before you start judging. Either that or threshold should depend on sample size and probability, but limiting sample size might be easier to understand for general public.
I guess the rest of us not in the hospital business can use this to be weary of statistics and data presented to us, even when provided from a reliable source. When gathering our own statistics its often necessary to do that for a prolonged period in order to let the data reveal what and what isn't pure chance.
The real take-home here is that if you set your 'warning' threshold to be too low, you'll get a lot of warnings. Here, they used an arbitrary cutoff of 8 as their warning threshold. A fairer move would have been to allow some data to come back, try to infer a distribution and it's parameters, and then set future cutoffs based on those parameters. One could even continuously recalculate based on incoming data.
Of course for you and me maybe, but many people wouldn't think of that.
Their book "The Tiger That Isn't" (http://www.amazon.co.uk/Tiger-That-Isnt-Through-Numbers/dp/1...) is excellent and includes some other examples around sample size and regression to mean.
One thing that struck me from the article was how computers and Monte Carlo simulations have lowered the barrier to entry for statistical understanding. It seems much simpler to run this sort of simulation than to understand the intricacies of different probability distributions.
So the number of deaths in a hospital is just
sum([random()<deathRate for i in range(numPeople)])
and they can do the analysis without needing to know what a Bernouilli or Binomial distribution is, even though these both feature in the problem.In the code above the random number is uniform(0,1), the standard random number in most languages. So each individual death is Bernouilli, although you don't need to know the word Bernouilli. Then this leads to a binomial distribution for a hospital as a whole.
But the binomial assumption doesn't need to be known or explicitly built into the model. It's effectively an emergent property of the individual deaths.