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))