2013-04-17 48 views
11

我已經運行邏輯迴歸,其名稱的摘要。 「分數」因此,summary(score)給了我下面的僅從邏輯模型中提取其p值有意義的係數

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。

理想情況下,我想讓他們在數據框中。不幸的是,我試過的下面的代碼不起作用。

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," 

回答

18

這個怎麼樣:

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

謝謝!很有幫助。 –

+2

@JonathanRossCharlton重新命名問題;請**不要**那樣做。這對任何人都有用? –

+3

@JonathanRossCharlton除非用戶明確表示他們否決投票,否則無法知道是誰投了票。投票是完全匿名的。無論如何,諸如「R編程」和特定用戶的提及在標題中被認爲是不適當的。 – joran