如何merge

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 ~ rank, data = mydata, family = "binomial")
   
7. usual.output <-summary(regress.out)
   
8. usual.output
   
9. 
   
10. data<-data.frame(variable=c("rank","","",""),x=c(1,2,3,4),x1=c(10,10,10,10))
    
11. 希望得到的结果
    
12. variable x x1   coef        se        z
    
13. rank    1 10     0             -         -
    
14.            2 10   -0.7500 0.3080 -2.435
    
15.            3 10   -1.3647 0.3354 -4.069
    
16.            4 10   -1.6867 0.4093 -4.121

复制代码

如何将logistic回归的结果merge到数据框data，谢谢！

data[,c("coef")]<-c("0",format(usual.output$coefficients[2:4,1],digits=4,trim=F,nsmall=4))
data[,c("se")]<-c("-",format(usual.output$coefficients[2:4,2],digits=4,trim=F,nsmall=4))
data[,c("z")]<-c("-",format(usual.output$coefficients[2:4,3],digits=3,trim=F,nsmall=3))

R这么复杂啊？

直接cbind不就行了吗？

qoiqpwqr 发表于 2013-4-29 04:36
直接cbind不就行了吗？

rank    1 10     0             -         -

楼主，你的结果我弄出来了，为何不兑现１００金帀啊