# ロシアのワクチンは旧ソ連の人工衛星にちなみ「スプートニクＶ」と名付けられた。
# ランセットに掲載された論文によると、
# ７６人を対象に第１相および第２相の治験を行ったところ、
# 全員が新型コロナウイルスに対する抗体を獲得したという。
# 米Ｊ＆Ｊ、コロナワクチン治験中断　参加者が原因不明の病気
# 同社は９月、新型コロナワクチンの治験の最終段階に当たる第３相試験に着手したと発表。
# 最大６万人を対象に安全性や有効性を確認するとしていた。

# 治験が再開されて6万例でこの奇病1例だけと仮定して発生割合を1/60000とする。
# ロシアのワクチンではこの奇病の報告はないから発生割合は0/76

# 奇病発生率1/60000と0/76では、どちらが奇病発生が少ないといえるか？

rm(list=ls())
source('tools.R')

# Jeffereys prior
a=b=0.5
par(mfrow=c(2,1))
curve(dbeta(x,a+0,b+76),bty='l')
curve(dbeta(x,a+1,b+59999),col=2,bty='l')
k=1e6

SP=rbeta(k,a+0,b+76)
JJ=rbeta(k,a+1,b+59999)
mean(SP/JJ>1)
mean(SP-JJ>0)

jags4prop(0,1,76,60000,0.5,0.5)