我在弄清楚如何解决此问题时遇到了一些困难。 我有一个数据框,我将其拆分为不同的站点(link5)。 拆分后,我基本上想对子集运行线性回归模型。 这是我正在使用的代码,但是绝对不正确。 另外,如果我可以将模型结果输出到新的数据框中,使每个站点都具有一行带有模型参数估计值的行,那将非常好,这只是一个愿望,而现在不是必须的。 感谢您的任何帮助!
les_events <- split(les, les$link5)
result <- lapply(les_events) {
lm1 <-lm(cpe~K,data=les_events)
coef <- coef(lm1)
q.hat <- -coef(lm1)[2]
les_events$N0.hat <- coef(lm1[1]/q.hat)
}