2013-04-17 7 views
11

Ho eseguito una regressione logistica, il cui riepilogo è il nome. "Punteggio" Di conseguenza, summary(score) mi dà il seguenteEstrai solo i coefficienti i cui valori p sono significativi da un modello logistico

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 

Quello che spero di raggiungere è quello di ottenere variabili nomi e coefficienti di quelle variabili che hanno un *, **, o *** accanto al loro valore di Pr(>|z|). In altre parole, voglio le variabili e i coefficienti sopra citati con un Pr(>|z|) < .05.

Idealmente, mi piacerebbe averli in un frame di dati. Sfortunatamente, il seguente codice che ho provato non funziona.

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

risposta

18

Che dire di questo:

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

Grazie! Molto utile. –

+2

@JonathanRossCharlton Re la ridenominazione della domanda; per favore ** non farlo **. A cosa serve questo a nessuno? –

+3

@JonathanRossCharlton A meno che un utente dichiari esplicitamente di non aver votato, non hai modo di sapere chi lo ha fatto. Le votazioni sono completamente anonime. Indipendentemente da ciò, cose come "Programmazione R" e menzioni di utenti specifici sono considerate inappropriate nei titoli. – joran