如何在R中提取stan_glm模型的sigma参数值？

有两种简单方法可以获取sigma的中位数：

方法一：通过模型汇总结果提取 #

使用summary()函数生成模型的完整汇总统计，然后直接提取sigma对应的中位数（即50%分位数）：

# 生成模型汇总
fit_summary <- summary(fit)
# 提取sigma的中位数
sigma_median <- fit_summary$summary["sigma", "50%"]

方法二：提取后验样本计算中位数 #

先提取sigma的所有后验样本，再用median()函数计算中位数：

# 提取sigma的后验样本
sigma_samples <- posterior_samples(fit, pars = "sigma")$sigma
# 计算中位数
sigma_median <- median(sigma_samples)

控制台输入fit时显示的是模型的简略汇总信息，上述两种方法都能精准获取你看到的那个2.7左右的中位数。

内容的提问来源于stack exchange，提问作者slow_learner