4

以下に示すように、3 つの時系列を含む data.frame があります。より滑らかな時系列でそれらをプロットすると、プロットする線形モデルのパラメーターを取得できるようになりたいのですが、その方法がわかりませんか? 

> data
day   od series_id
1    1 0.10        A1
2    3 1.00        A1
3    5 0.50        A1
4    7 0.70        A1
5    1 1.70        B1
6    3 1.60        B1
7    5 1.75        B1
8    7 1.70        B1
9    1 2.10        C1
10   3 2.30        C1
11   5 2.50        C1
12   7 2.70        C1

data = data.frame (day = c(1,3,5,7,1,3,5,7,1,3,5,7), 
od = c(0.1,1.0,0.5,0.7 ,1.7,1.6,1.75,1.7 ,2.1,2.3,2.5,2.7), 
series_id = c("A1", "A1", "A1","A1", "B1", "B1","B1", "B1", "C1","C1", "C1", "C1"))

r <- ggplot(data = data, aes(x = day, y = od))
r + stat_smooth(aes(group = series_id, color = series_id),method="lm")
4

1 に答える 1

4

レイヤーオブジェクトが将来すべての計算結果を保存する必要があるかどうかはわかりませんが、現在、描画時にフィッティングが呼び出され、モデルはレイヤーオブジェクトに保存されません。簡単な回避策として、次の方法で同じ結果を得ることができます

 dlply(data,.(series_id),function(x)lm(od~day,data=x))
于 2010-05-01T01:53:33.317 に答える