hope this will help
Eg:
wd<-choose.dir()
setwd(wd)

counts <- c(579, 4219, 911, 4578, 246, 3775, 361, 4532, 291, 5261, 396, 6598)
group <- gl(2, 1, 12, labels=c("cases","controls"))
blood <- gl(2, 2, 12, labels=c("A","O"))
city  <- gl(3, 4, 12, labels=c("London","Manchester","Newcastle"))
cbind(group, blood, city, counts)


z1 <- glm(counts ~ group*(city + blood), family=poisson())

names(z1)   ## type to show what items you can use

write.csv(z1$coefficients,"test.csv",row.names=F,col.names=F)

write.csv(t(z1$coefficients),"test2.csv",row.names=F,col.names=F)

I think what you need is:

model<-lm(...)

resd<-model$residuals
fitted<-model$fitted

output<-cbind("residuals"=resd,"fitted"=fitted)
write.csv(output,"outcome.csv")

为什么会出现错误呢？

警告信息：
In write.csv(z1$coefficients, "test.csv", row.names = F, col.names = F) :
  不能修改'col.names'


O(∩_∩)O谢谢

如果不仅仅输出系数，还要输出se，t或z值，p值怎么做呢？总之输出一个论文中呈现的回归报告结果该怎么做呀？

你好请问这个问题解决了嘛，我也想输出论文报告中的回归结果，另外能直接生成需要的表格形式么