Skip to content

Select components of a list

1 message · arun

#
HI Gustav,

If you need the combined output:
res<-lapply(1:length(models),function(i) do.call(rbind,lapply(models[[i]],function(x) summary(x)$coef[row.names(summary(x)$coef)%in%c("pm10","ozone","so2"),c(1:2,4)])))
?names(res)<-1:length(res)
res1<-do.call(rbind,lapply(res,function(i) {row.names(i)<-c("pm10","ozone","so2");data.frame(i)}))
names(res1)[2:3]<- c("Std.Error","Pr(>|z|)")
res1
#??????????? Estimate??? Std.Error???? Pr(>|z|)
#1.pm10? 0.0005999185 0.0001486195 5.423004e-05
#1.ozone 0.0010117294 0.0003792739 7.640816e-03
#1.so2?? 0.0026595441 0.0009352046 4.457766e-03
#2.pm10? 0.0005720549 0.0001740368 1.012696e-03
#2.ozone 0.0009128304 0.0004364390 3.647954e-02
#2.so2?? 0.0028256121 0.0010150314 5.373144e-03
#3.pm10? 0.0005099552 0.0001559620 1.076462e-03
#3.ozone 0.0023896044 0.0004109854 6.087769e-09
#3.so2?? 0.0024097381 0.0009563814 1.174744e-02
#4.pm10? 0.0009285593 0.0001766520 1.468764e-07
#4.ozone 0.0005455392 0.0004301502 2.047076e-01
#4.so2?? 0.0017251400 0.0011635156 1.381552e-01
A.K.