Bayesian Computation with R 4.2 パラメータが二つとも未知の正規データ > data(marathontimes) > attach(marathontimes) > d = myconto
階層モデリング
7.5 交換可能性を事前の確信とするモデリング 事前分布として、真の死亡率λ1,…,λ94はガンマ分布gamma(α, α/μ)からのランダムな標本だと仮定する。 g(λ|α,μ) λの事前平均と分散は、それ
マルコフ連鎖モンテカルロ法
Bayesian Comuputation with R 6.3 メトロポリス-ヘイスティング・アルゴリズム <やっとたどり着いたMCMC法!> 事後密度g(θ|y)からシミュレーションを行う。 初期値θ0か
ベイズ計算入門2:棄却サンプリングと重点サンプリング
5.8 棄却サンプリング 棄却サンプリングrejection sampling: 所与の確率分布からランダム標本をシミュレーションする汎用的なアルゴリズム <棄却サンプリングのステップ> pから標本のシミュレーションが容
ベイズ計算入門
From Baysean Computation with R 5.2 積分を計算する f(y|θ)からデータyを観測する。ここでθは、パラメータベクトルであり、その事前分布としてg(θ)を割り当てる。 θの事後密度は、
1パラメータモデル
From Bayesian Computation with R 3.2 平均が既知で分散が未知の正規分布 平均0で未知分散σ^2の正規分布標本とすると、 尤度関数:L(σ^2) ∝ (σ^2)^(-n/2)exp{-Σ
「ベイズ的思考への誘い」への考察
「ベイズ的思考への誘い」では、アメリカ人の大学生が8時間以上の睡眠を確保しているかどうかについて、ベイズ的統計学的解析が加えられた。 解析は、8時間以上睡眠を取れている場合をs、そうでない場合はfとして、全体の学生数に対
ベイズ的思考への誘い 練習問題2.9
1. > p = seq(0.0, 1.0, by = 0.125) > prior = c(0.001, 0.001, 0.95, 0.008, 0.008, 0.008, 0.008, 0.008, 0.008) >
ベイズ的思考への誘い
From Bayesian Computation with R アメリカ大学生の睡眠習慣。何割の学生が最低8時間の睡眠を取っているか?8時間の睡眠を取っている学生の割合をpとする。 ベイズでは、事前分布として主観を使え
Bayesian Computation with R Practice #1-6 その2
練習問題1-6 4. a) binomial.conf.interval y y 13 > y_90 y_90 0.5133179 0.7866821 > y y 15 12 14 12 12 1