r - 贝叶斯可信区间的覆盖概率估计(正态分布)
问题描述
正态分布的贝叶斯推理,我使用以下r
代码获得后验分布。
install.packages(c("mvtnorm","loo","coda"), repos="https://cloud.r-project.org/",dependencies=TRUE)
options(repos=c(getOption('repos'), rethinking='http://xcelab.net/R'))
install.packages('rethinking',type='source')
library(rethinking)
set.seed(650)
x <- data.frame(x = rt(100,3))
fit <- rethinking::map(
alist(
x ~ dnorm(mu, sigma),
mu ~ dnorm(1, 10),
sigma ~ dunif(0, 50)
),
data=x)
precis(fit, corr=TRUE)
sim_post <- extract.samples(fit)
dim(sim_post)
post_mean <- apply(sim_post, 2, mean)
post_mean
quantile(sim_post$mu , c(.05, .95))
quantile(sim_post$sigma, c(0.05, 0.95))
我使用后验分布模拟了 10000 个样本。如何获得 mu 的覆盖概率?
解决方案
推荐阅读
- linux - Passbolt 安装失败:进程 875 当前正在使用 usermod www-data
- batch-file - 为什么批量将从文件加载的变量设置为0?
- javascript - 在js文件中添加sqlite3模块后刷新页面时出现空白电子窗口
- php - 为什么这个 PHP Curl PUT 请求返回代码 500?
- sed - 使用 sed 查找查询以删除除上次出现之外的所有条目
- java - 如果我的线程不改变地图的结构,我应该使用 ConcurrentHashMap 还是 HashMap?
- java - 我不断收到一种类型的错误:“表达式的类型必须是数组类型,但它解析为长”
- python - 类型别名和 NewType 的区别
- css - CSS 选择器颜色优先级
- python - hop_length、win_length、frame_length、n_fft、帧数之间的关系