数学 統計に詳しい人が語るコロナウイルスat MATH数学 統計に詳しい人が語るコロナウイルス - 暇つぶし2ch■コピペモード□スレを通常表示□オプションモード□このｽﾚｯﾄﾞのURL■項目テキスト650:１３２人目の素数さん 20/05/23 06:44:49 COZ69QMb.net bootstrap使うと大抵の数値の信頼区間は出せるけど1/1には無理だな。 651:１３２人目の素数さん 20/05/23 15:08:38 COZ69QMb.net JAGSを使って1/1のときの95％ Credibility Intervalを求めるスクリプト 事前分布はJeffereys library(rjags) data=list(x=1,n=1) cat(' model{ x ~ dbin(p,n) # binomial distribution p ~ dbeta(0.5,0.5) # Jeffereys prior }',file='tmp.txt') jagsModel=jags.model('tmp.txt',data=data,n.chains=4) update(jagsModel) codaSamples=coda.samples(jagsModel,n.iter=1e5,var='p') gelman.plot(codaSamples) # 収束を確認 js=as.data.frame(as.matrix(codaSamples)) HDInterval::hdi(js$p) # 95% Highest Densty Interval # 既存のパッケージで確認。 HDInterval::hdi(qbeta,shape1=0.5+1,shape2=0.5) binom::binom.bayes(1,1) 652:１３２人目の素数さん 20/05/23 17:34:11 COZ69QMb.net エクセルのツールをCDCが公開しているね 次ページ最新レス表示レスジャンプ類似スレ一覧スレッドの検索話題のニュースおまかせリストオプションしおりを挟むスレッドに書込スレッドの一覧暇つぶし2ch