如何从 stan_glm 中的系数中提取标准误差

Question

我有一个看起来像这样的模型：

fit <- stan_glm(switch ~ arsenic + dist100, 
                family=binomial(link="logit"), 
                data=wells, 
                refresh=0)

输出如下所示：

stan_glm
 family:       binomial [logit]
 formula:      switch ~ arsenic + dist100
 observations: 3020
 predictors:   3
------
            Median  MAD_SD 
(Intercept)  0.0054  0.0810
arsenic      0.4618  0.0415
dist100     -0.9025  0.1040

------

我想要 MAD_SD 结果拟合的系数。我通常使用 se(fit) 并给出系数，但由于某种原因我得到错误

Error in se(fit) : argument "predicted" is missing, with no default

我该如何补救这个错误？

Answer 1

如果有人想知道这个问题的答案，stan_glm objects 基本上就像将名称映射到元素的大标题。在本例中，该元素名为“ses”，它是一个列表。所以你可以做

fit$ses[0]

对其进行索引

编辑：事实证明，我需要这样做的原因是因为我在 stan_glm 之后导入了一个包，它也有一个名为“se”的函数并且它被屏蔽了。

如何从 stan_glm 中的系数中提取标准误差

How to extract standard error from coefficients in stan_glm

r

glm

stan