"
餅を1億口食べると7人が窒息死するという
https://youtu.be/sSCM9mRAsyc?t=315

某AIの答
>>
日本でワクチンが原因と認定された死亡例は、約79件でした。この中で、ワクチン接種後に死亡が確認されたケースは51件でした1。
一方、日本での総ワクチン接種回数は約4億3961万回に上ります。
<<

【問題】
（１）餅による窒息死とワクチン原因死のリスク比（もしくはオッズ比）の95%信頼区間を算出せよ。

【博物館入りの古典問題】
餅による窒息死とワクチン原因死の確率は同じを帰無仮説として有意差検定せよ。
すなわち、帰無仮説を前提に実測値とmore extremeな場合の確率の合計=p値を求めよ。
"

source('toolmini.R')
r1=7
r2=51
n1=1e8
n2=4.3961e8
Fisher.test(c(r1,r2),c(n1,n2))
prop.test(c(r1,r2),c(n1,n2))

hit=c(r1,r2)
shot=c(n1,n2)
mat=cbind(hit,shot-hit)
Epi::twoby2(mat)

mochi=rbeta(1e9,7+0.5,1e8-7+0.5)
vacci=rbeta(1e9,51+0.5,4.3961e-51+0.5)
d=mochi-vacci
hist(d)
summary(d)
HDInterval::hdi(d)
mean(mochi>vacci)

re=beta.diff(r1,r2,n1,n2)
re$cdf(0)