# ある仮想の難治疾患患者25人従来薬を投与して3人治癒した。
# 新薬が登場して3人に投与したところ治癒した人はいなかった。
# この新薬を継続して使う価値があるかどうか検討せよ。
別バージョン
# 巨乳女子大で25人に声をかけたら3人が誘いにのった。
# 桃尻女子大で3人に声をかけたら誰も誘いにのらなかった。
# どちらが口説きやすいか検討せよ。

JAGSでMCMCして治癒率の確率密度関数を描くとこうなる。
http://i.imgur.com/y49H5AK.png

治癒率差の不偏推定量は
> mean(dif)
[1] -0.05136971
54％が負
> c(mean(dif<0),mean(dif>0))
[1] 0.5395 0.4605
５％幅の違いは同等扱いにすると
> c(mean(dif<ROPE[1]),mean(ROPE[1]<dif & dif<ROPE[2]), mean(dif>ROPE[2]))
[1] 0.4834 0.1236 0.3930
と計算できる。
98％HDIは
> HDInterval::hdi(dif)
lower upper
-0.4247349 0.2535357
と０を挟む。

RのパッケージBESTを改造して、治癒率の差の確率密度をかくと
http://i.imgur.com/vdIj7ES.png

ゆえに新薬は無効とはいえないだけでなく、不偏推定量から従来薬を凌駕する可能性が54％ある。