从逻辑模型中仅提取p值显着的系数 [英] Extract only coefficients whose p values are significant from a logistic model

问题描述

我进行了逻辑回归，并总结了我的名字.因此，summary(score)给了我以下内容

I have run a logistic regression, the summary of which I name. "score" Accordingly, summary(score) gives me the following

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.3616  -0.9806  -0.7876   1.2563   1.9246  

                       Estimate Std. Error    z value    Pr(>|z|)
(Intercept)        -4.188286233 1.94605597 -2.1521921 0.031382230 *
Overall            -0.013407201 0.06158168 -0.2177141 0.827651866
RTN                -0.052959314 0.05015013 -1.0560154 0.290961160
Recorded            0.162863294 0.07290053  2.2340482 0.025479900 *
PV                 -0.086743611 0.02950620 -2.9398438 0.003283778 **
Expire             -0.035046322 0.04577103 -0.7656878 0.443862068
Trial               0.007220173 0.03294419  0.2191637 0.826522498
Fitness             0.056135418 0.03114687  1.8022810 0.071501212 .

---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 757.25  on 572  degrees of freedom
Residual deviance: 725.66  on 565  degrees of freedom
AIC: 741.66

Number of Fisher Scoring iterations: 4

我希望实现的是获取变量名称和这些变量的系数，这些变量的Pr(>|z|)值旁边有*，**或***.换句话说，我希望上述的变量和系数的Pr(>|z|)< .05.

What I am hoping to achieve is to get variables names and coefficients of those variables which have a *, **, or *** next to their Pr(>|z|) value. In other words, I want the aforementioned variables and coefficients with a Pr(>|z|) < .05.

理想情况下，我想将它们放入数据框.不幸的是，我尝试过的以下代码不起作用.

Ideally, I'd like to get them in a data frame. Unfortunately, the following code I've tried does not work.

variable_try <-
  summary(score)$coefficients[if(summary(score)$coefficients[, 4] <= .05, 
                                 summary(score)$coefficients[, 1]),]

Error: unexpected ',' in "variable_try <-
summary(score)$coefficients[if(summary(score)$coefficients[,4] < .05,"

推荐答案

这是怎么回事:

data.frame(summary(score)$coef[summary(score)$coef[,4] <= .05, 4])

这篇关于从逻辑模型中仅提取p值显着的系数的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！