后验的更新

关于您选举支持率参数 \(p\) 的后验模型，同时综合了 \(p\) 的先验模型与民调数据 \(X\) 的信息。请运行右侧脚本，回顾在最初先验（Beta(45, 55)）与最初民调数据（在抽样的 \(n = 10\) 名选民中有 \(X = 6\) 支持您）下得到的后验。工作区中已提供 vote_model。

接下来通过 3 个步骤，您将探索采用「不同的先验模型」或观测到「新的数据」（或二者结合）会如何影响后验。

本练习是课程的一部分

使用 RJAGS 的贝叶斯建模

查看课程

交互式实操练习

通过完成这段示例代码来试试这个练习。

# COMPILE the model    
vote_jags <- jags.model(textConnection(vote_model), 
    data = list(a = 45, b = 55, X = 6, n = 10),
    inits = list(.RNG.name = "base::Wichmann-Hill", .RNG.seed = 100))

# SIMULATE the posterior
vote_sim <- coda.samples(model = vote_jags, variable.names = c("p"), n.iter = 10000)

# PLOT the posterior
plot(vote_sim, trace = FALSE, xlim = c(0,1), ylim = c(0,18))

编辑并运行代码