我试图同时从列表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)
}
}