请问glm估计出来的模型怎么看伪R平方？

请问glm估计出来的模型怎么看伪R平方？用summary只有参数的显著性，但是没有R平方的结果。

1. Call:
   
2. glm(formula = GR ~ C + F1 + I, family = binomial(link = "probit"),
   
3. data = data)
   
4. 
   
5. Deviance Residuals:
   
6. Min 1Q Median 3Q Max
   
7. -1.2771 -0.4532 -0.3646 -0.2682 2.3106
   
8. 
   
9. Coefficients:
   
10. Estimate Std. Error z value Pr(>|z|)
    
11. (Intercept) 2.37052 3.62695 0.654 0.51338
    
12. C -0.03203 0.04870 -0.658 0.51080
    
13. F1 -1.40809 0.54465 -2.585 0.00973 **
    
14. I 6.87745 605.10217 0.011 0.99093
    
15. ---
    
16. Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
    
17. 
    
18. (Dispersion parameter for binomial family taken to be 1)
    
19. 
    
20. Null deviance: 40.472 on 40 degrees of freedom
    
21. Residual deviance: 28.160 on 37 degrees of freedom
    
22. AIC: 36.16
    
23. 
    
24. Number of Fisher Scoring iterations: 15

复制代码

1. 需要自己计算

复制代码

1. install.package("AER")
   
2. 
   
3. library("AER")
   
4. data("SwissLabor")
   
5. 
   
6. swiss_probit <- glm(participation ~ . + I(age^2),
   
7.                     data = SwissLabor, family = binomial(link = "probit"))
   
8. summary(swiss_probit)
   
9. 
   
10. swiss_probit0 <- glm(participation ~ 1 ,
    
11.                     data = SwissLabor, family = binomial(link = "probit"))
    
12. summary(swiss_probit0)
    
13. 
    
14. #Goodness of fit and prediction
    
15. pseudo_R2<-  1 - as.vector(logLik(swiss_probit)/logLik(swiss_probit0))
    
16. pseudo_R2

复制代码

这是Applied Econometrics with R 书上的例子

http://stats.stackexchange.com/questions/8511/how-to-calculate-pseudo-r2-from-rs-logistic-regression

答案里面提供了2种方法， 比较直接的是用 rms 软件包

1. require(rms)
   
2. mod1b <- lrm(y ~ x)
   
3. printed(mod1b)

复制代码

Aharach 发表于 2014-4-18 22:09
http://stats.stackexchange.com/questions/8511/how-to-calculate-pseudo-r2-from-rs-logistic-regression ...

那个网址怎么打不开？
下面那个代码怎么用呢？比如我已经用glm估计除了probit，怎么显示他的R平方呢？

qiangli 发表于 2014-4-18 20:39
这是Applied Econometrics with R 书上的例子

哦 非常感谢！

马克一下~