การประมาณค่าและการอนุมานแบบ Posterior

การสร้างแบบจำลอง Bayesian ด้วย RJAGS

Alicia Johnson

Associate Professor, Macalester College

โมเดลการถดถอยแบบ Bayesian

$Y_i$ = น้ำหนักของผู้ใหญ่คนที่ $i$ (กก.)
$X_i$ = ส่วนสูงของผู้ใหญ่คนที่ $i$ (ซม.)

$\;$

โมเดล

$Y_i \sim N(m_i, s^2)$
$m_i = a + b X_i$

$a \sim N(0, 200^2)$
$b \sim N(1, 0.5^2)$
$s \sim \text{Unif}(0, 20)$

การประมาณค่าจุดของ Posterior

summary(weight_sim_big)

1. Empirical mean and standard deviation for each variable,
   plus standard error of the mean:
      Mean      SD  Naive SE Time-series SE
a -104.038 7.85296 0.0248332       0.661515
b    1.012 0.04581 0.0001449       0.003849
s    9.331 0.29495 0.0009327       0.001216

2. Quantiles for each variable:
       2.5%       25%      50%     75%   97.5%
a -118.6843 -109.5171 -104.365 -99.036 -87.470
b    0.9152    0.9828    1.014   1.044   1.098
s    8.7764    9.1284    9.322   9.524   9.933

ค่าเฉลี่ย Posterior ของ $a$ $\approx$ -104.038

ค่าเฉลี่ย Posterior ของ $b$ $\approx$ 1.012

การประมาณค่าจุดของ Posterior

แนวโน้มค่าเฉลี่ย Posterior:

$m_i = -104.038 + 1.012 X_i$

ผลลัพธ์ Markov chain:

head(weight_chains)

             a        b        s
[1,] -113.9029 1.072505 8.772007
[2,] -115.0644 1.077914 8.986393
[3,] -114.6958 1.077130 9.679812
[4,] -115.0568 1.072668 8.814403
[5,] -114.0782 1.071775 8.895299
[6,] -114.3271 1.069477 9.016185

ความไม่แน่นอนของ Posterior

แนวโน้มค่าเฉลี่ย Posterior:

$m_i = -104.038 + 1.012 X_i$

ผลลัพธ์ Markov chain:

head(weight_chains)

             a        b        s
[1,] -113.9029 1.072505 8.772007
[2,] -115.0644 1.077914 8.986393
[3,] -114.6958 1.077130 9.679812
[4,] -115.0568 1.072668 8.814403
[5,] -114.0782 1.071775 8.895299
[6,] -114.3271 1.069477 9.016185

ช่วงความน่าเชื่อถือของ Posterior

summary(weight_sim_big)

1. Empirical mean and standard deviation for each variable,
   plus standard error of the mean:
      Mean      SD  Naive SE Time-series SE
a -104.038 7.85296 0.0248332       0.661515
b    1.012 0.04581 0.0001449       0.003849
s    9.331 0.29495 0.0009327       0.001216

2. Quantiles for each variable:
       2.5%       25%      50%     75%   97.5%
a -118.6843 -109.5171 -104.365 -99.036 -87.470
b    0.9152    0.9828    1.014   1.044   1.098
s    8.7764    9.1284    9.322   9.524   9.933

ช่วงความน่าเชื่อถือ Posterior 95% ของ $a$: (-118.6843, -87.470)

ช่วงความน่าเชื่อถือ Posterior 95% ของ $b$: (0.9152, 1.098)

ช่วงความน่าเชื่อถือของ Posterior

การตีความ

เมื่อพิจารณาจาก Prior และข้อมูลที่สังเกตได้ มีโอกาส 95% (Posterior) ที่ $b$ จะอยู่ระหว่าง 0.9152 ถึง 1.098 กก./ซม.

ความน่าจะเป็นแบบ Posterior

table(weight_chains$b > 1.1)

FALSE  TRUE 
97835  2165

mean(weight_chains$b > 1.1)

0.02165

การตีความ:
มีโอกาส Posterior 2.165% ที่ $b$ จะมากกว่า 1.1 กก./ซม.

มาฝึกกันเถอะ!

การสร้างแบบจำลอง Bayesian ด้วย RJAGS