Advertisement
gronke

Untitled

Jul 3rd, 2014
183
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. > summary(mydata)
  2.       Won              Club          Healthcare         Other        
  3.  Min.   :0.0000   Min.   :0.0000   Min.   :0.0000   Min.   :0.00000  
  4.  1st Qu.:0.0000   1st Qu.:0.0000   1st Qu.:0.0000   1st Qu.:0.00000  
  5.  Median :0.0000   Median :0.0000   Median :0.0000   Median :0.00000  
  6.  Mean   :0.2731   Mean   :0.1101   Mean   :0.1359   Mean   :0.09723  
  7.  3rd Qu.:1.0000   3rd Qu.:0.0000   3rd Qu.:0.0000   3rd Qu.:0.00000  
  8.  Max.   :1.0000   Max.   :1.0000   Max.   :1.0000   Max.   :1.00000  
  9.    Government    
  10.  Min.   :0.0000  
  11.  1st Qu.:0.0000  
  12.  Median :1.0000  
  13.  Mean   :0.6567  
  14.  3rd Qu.:1.0000  
  15.  Max.   :1.0000
  16.  
  17. > str(mydata)
  18. 'data.frame':   7981 obs. of  5 variables:
  19.  $ Won       : int  1 0 0 0 0 0 0 0 0 0 ...
  20.  $ Club      : int  0 0 0 0 0 0 1 1 0 0 ...
  21.  $ Healthcare: int  0 0 0 0 0 0 0 0 0 0 ...
  22.  $ Other     : int  0 0 0 0 0 0 0 0 0 0 ...
  23.  $ Government: int  1 1 1 1 1 1 0 0 1 1 ...
  24.  
  25. > glm.out = glm(w ~ c+h+o, family=binomial(logit), data=mydata)
  26. > summary(glm.out)
  27.  
  28. Call:
  29. glm(formula = w ~ c + h + o, family = binomial(logit), data = mydata)
  30.  
  31. Deviance Residuals:
  32.     Min       1Q   Median       3Q      Max  
  33. -1.1649  -0.7063  -0.7063   1.1900   1.7382  
  34.  
  35. Coefficients:
  36.             Estimate Std. Error z value Pr(>|z|)    
  37. (Intercept) -0.02958    0.06747  -0.438    0.661    
  38. c           -0.79804    0.09437  -8.456  < 2e-16 ***
  39. h           -0.65584    0.10166  -6.451 1.11e-10 ***
  40. o           -1.23167    0.07524 -16.370  < 2e-16 ***
  41. ---
  42. Signif. codes:  0***0.001**0.01*0.05 ‘.’ 0.1 ‘ ’ 1
  43.  
  44. (Dispersion parameter for binomial family taken to be 1)
  45.  
  46.     Null deviance: 9359.6  on 7980  degrees of freedom
  47. Residual deviance: 9074.3  on 7977  degrees of freedom
  48. AIC: 9082.3
  49.  
  50. Number of Fisher Scoring iterations: 4
  51.  
  52. > anova(glm.out, test="Chisq")
  53. Analysis of Deviance Table
  54.  
  55. Model: binomial, link: logit
  56.  
  57. Response: w
  58.  
  59. Terms added sequentially (first to last)
  60.  
  61.  
  62.      Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
  63. NULL                  7980     9359.6              
  64. c     1    5.965      7979     9353.6   0.01459 *  
  65. h     1   19.076      7978     9334.5 1.256e-05 ***
  66. o     1  260.238      7977     9074.3 < 2.2e-16 ***
  67. ---
  68. Signif. codes:  0***0.001**0.01*0.05 ‘.’ 0.1 ‘ ’ 1
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement