帮我理解pp- 贝叶斯 glm 中的值

机器算法验证 r 贝叶斯 p 值
2022-03-07 04:27:34

我正在尝试对此处的数据运行贝叶斯 logit 。我在 Rbayesglm()arm包中使用。编码很简单:

df = read.csv("http://dl.dropbox.com/u/1791181/bayesglm.csv", header=T)
library(arm)
model = bayesglm(PASS ~ SEX + HIGH, family=binomial(link="logit"), data=df)

summary(model)给出以下输出:

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)  0.10381    0.10240   1.014    0.311    
SEXMale      0.02408    0.09363   0.257    0.797    
HIGH        -0.27503    0.03562  -7.721 1.15e-14 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 2658.2  on 1999  degrees of freedom
Residual deviance: 2594.3  on 2000  degrees of freedom
AIC: 2600.3

请引导我完成这个。我知道这段代码使用了一个非常弱的先验(因为我没有指定先验方法),所以如果我使用glm()而不是bayesglm(). 但是输出应该仍然是贝叶斯精神,对吧?什么是p-价值观和z-这里的值?这些不是常客推理工具吗?他们在这里有不同的解释吗?

1个回答

好问题!尽管有贝叶斯 p 值,并且 arm 包的作者之一是倡导者,但您在输出中看到的并不是贝叶斯 p 值。检查类model

class(model)
"bayesglm" "glm"      "lm" 

你可以看到类bayesglm 继承自glm。此外,对 arm 包的检查显示没有针对 bayesglm 对象的特定摘要方法。所以当你这样做时

summary(model)

你实际上在做

summary.glm(model)

并获得对结果的常客解释。如果您想要更多的贝叶斯视角,则 arm 中的功能是display()