0

我试图同时从列表 x 中的 5 个数据帧中获取所有五个输出,但我被困在这里。另外,如果我希望函数从每个输出返回特定的系数,比如说X24hr_PM25,该怎么办?谢谢!

这是代码:

q1<-cvd[cvd$quintile=='1',]
q2<-cvd[cvd$quintile=='2',]
q3<-cvd[cvd$quintile=='3',]
q4<-cvd[cvd$quintile=='4',]
q5<-cvd[cvd$quintile=='5',]

x<-list(q1,q2,q3,q4,q5)

analyze <- function(data, df) {
  for (i in 1:length(data)) {
    fit <- glm(cvd_o ~ X24hr_PM25 + ns(c7, 6*df) + ns(temp, 4) + ns(rh, 4) + dow,
               family=poisson, data=as.data.frame(data[i]))
    sum <- summary(fit)
    return(sum)
  }
}
4

0 回答 0