Uma explicação intuitiva por que o procedimento Benjamini-Hochberg FDR funciona?

Existe uma maneira simples de explicar por que o procedimento de Benjamini e Hochberg (1995) realmente controla a taxa de descoberta falsa (FDR)? Esse procedimento é tão elegante e compacto e, no entanto, a prova do porquê ele funciona sob independência (aparecendo no apêndice do artigo de 1995 ) não é muito acessível.

Aqui está um Rcódigo para gerar uma imagem. Ele mostrará 15 valores-p simulados plotados em relação à sua ordem. Então eles formam um padrão de ponto ascendente. Os pontos abaixo das linhas vermelha / roxa representam testes significativos no nível 0,1 ou 0,2. O FDR é o número de pontos pretos abaixo da linha dividido pelo número total de pontos abaixo da linha.

x0 <- runif(10)      #p-values of 10 true null hypotheses. They are Unif[0,1] distributed.
x1 <- rbeta(5,2,30)  # 5 false hypotheses, rather small p-values
xx <- c(x1,x0)
plot(sort(xx))
a0 <- sort(xx)
for (i in 1:length(x0)){a0[a0==x0[i]] <- NA}
points(a0,col="red")
points(c(1,15), c(1/15 * 0.1 ,0.1), type="l", col="red")
points(c(1,15), c(1/15 * 0.2 ,0.2), type="l", col="purple")

Espero que isso dê alguma idéia sobre a forma que a distribuição dos valores-p ordenados tem. O fato de as linhas estarem corretas e não, por exemplo, alguma curva em forma de parábola, tem a ver com o formato das distribuições de pedidos. Isso deve ser calculado explicitamente. De fato, a linha é apenas uma solução conservadora.