この問題にアプローチする方法を理解するのに苦労しています。個別のサイト (リンク 5) に分割しているデータ フレームがあります。分割したら、基本的にサブセットで線形回帰モデルを実行したいと思います。ここに私が取り組んでいるコードがありますが、それは間違いなく正しくありません。また、モデルの結果を新しいデータ フレームに出力して、各サイトにモデル パラメーターの推定値を含む 1 つの行が含まれるようにできれば素晴らしいと思います。これは単なる希望であり、現時点では必要ではありません。助けてくれてありがとう!
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)
}