如何增加一行

1. mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv")
   
2. library(foreign)
   
3. write.dta(mydata,"e:/abc.dta")
   
4. head(mydata)
   
5. mydata$rank <- factor(mydata$rank)
   
6. regress.out<- glm(admit ~ gre + gpa + rank, data = mydata, family = "binomial")
   
7. usual.output <-summary(regress.out)
   
8. usual.output
   
9. 
   
10. 输出的结果中Coefficients:
    
11.              Estimate Std. Error z value Pr(>|z|)   
    
12. (Intercept) -3.989979   1.139951  -3.500 0.000465 ***
    
13. gre          0.002264   0.001094   2.070 0.038465 *  
    
14. gpa          0.804038   0.331819   2.423 0.015388 *  
    
15. rank2       -0.675443   0.316490  -2.134 0.032829 *  
    
16. rank3       -1.340204   0.345306  -3.881 0.000104 ***
    
17. rank4       -1.551464   0.417832  -3.713 0.000205 ***
    
18. 
    
19. 
    
20. 把系数转成成数据框后，如何为多分类自变量多增加一行
    
21. Coefficients:
    
22.              Estimate Std. Error z value Pr(>|z|)   
    
23. (Intercept) -3.989979   1.139951  -3.500 0.000465 ***
    
24. gre          0.002264   0.001094   2.070 0.038465 *  
    
25. gpa          0.804038   0.331819   2.423 0.015388 *  
    
26. rank1        0
    
27. rank2       -0.675443   0.316490  -2.134 0.032829 *  
    
28. rank3       -1.340204   0.345306  -3.881 0.000104 ***
    
29. rank4       -1.551464   0.417832  -3.713 0.000205 ***
    
30. 谢谢！

复制代码

x <- rbind(usual.output$coefficients[1:3, ], c(0, rep(NA, 3)), usual.output$coefficients[-(1:3), ])
                rownames(x)[4] <- "rank1"
               
                printCoefmat(x, na.print="")

ntsean 发表于 2013-4-27 10:49
x

如果有多个多分类的自变量，如何自动实现呢？