高校数学の質問スレ Part438

# 一峰性非対称分布する離散量t1の95% highest probability intervalを返す関数
hpdi=\(t1){
t2=sort(table(t1),decreasing=TRUE)
t3=as.numeric(names(t2))
t4=t2/sum(t2)
t5=sum(cumsum(t4)<0.95)
ci1=range(t3[1:t5])
p1= sum(t2[1:t5]/sum(t2))
ci2=range(t3[1:(t5+1)])
p2=sum(t2[1:(t5+1)]/sum(t2))
list('ci.95-'=ci1,'p.95-'=p1,'ci.95+'=ci2 ,'p.95+'=p2)
}

ポアソン分布に従う乱数を発生させて上記に代入すれば算出してくれる。
Phimoseくんのキーキー電卓では算出できないみたいだな。