<分区>
我有一个数据框如下:
df
Date Hour ID Par1 Par2 Par3
08-01-15 0 A 2 3 4
08-01-15 0 B 4 5 6
08-01-15 1 N 2 9 10
08-01-15 1 A 3 7 23
08-01-15 1 B 4 7 22
08-02-15 0 E 2 4 12
08-02-15 0 A 3 7 9
所以我想按小时拆分此数据框,如下所示:
splitdata<-split(df<-split(df, df$Hour)
拆分后,我想对拆分后的数据集应用线性模型。
result <- lapply(splitdata, function(df){
lm1 <-lm(Par1~Par2,data=df)
summary <- (lm1$summary)
data.frame(as.list(summary))
})
result
我的结果没有显示任何内容。尽管如果我将摘要更改为:
summary <- (lm1$coef)
data.frame(as.list(summary))
然后它会产生一个结果。
所以主要问题是,如何按小时获取每个线性模型的摘要列表,而不仅仅是系数?
谢谢!