特定の範囲内の複数の列に対して線形回帰をプロットする

Question

私はいくつかの列を持つリストを持っています。左端の列はx軸のデータであり、残りはさまざまなyのデータです。

これらの列のそれぞれに対して線形回帰を生成したいが、xの特定の範囲についてのみ、線形回帰を0,0に強制したい。最後に、前記線形回帰の線勾配を抽出する。以下は私のデータの例です。ここ

x  y1  y2  y3  y4 
1 2.495 -1e-04 -1e-03 -1e-04 0e+00 
2 2.995 -2e-04 -7e-04 -2e-04 0e+00 
3 3.495 2e-04 -2e-04 1e-04 2e-04 
4 3.995 0e+00 -3e-04 0e+00 1e-04 
5 4.495 0e+00 -3e-04 0e+00 0e+00 
6 4.995 0e+00 -3e-04 0e+00 1e-04 
7 5.495 1e-04 -2e-04 0e+00 1e-04 
8 5.995 -1e-04 -1e-04 0e+00 3e-04 
9 6.495 1e-04 1e-04 2e-04 3e-04 
10 6.995 1e-04 0e+00 0e+00 3e-04 

は、例えば、私は0,0と交差する回帰を強制する、例えば行のxの3〜10の範囲のための線形回帰を得ることができます。最後に、y1、y2等に対する線形回帰からの傾きが得られる。

私はabind関数とlm関数を認識していますが、0,0の交差を強制せずに単一の回帰をプロットすることしかできませんでした。ありがとう

Answer 1

インターセプトのない回帰モデルは、回帰直線を原点に強制します。これは、Rにモデルに-1を加えることによって達成される。あなたが望むxの範囲を取得するには、単にその範囲に観測をサブセットと異なる独立変数をループにlapplyを使用します。

dat2 <- dat[3:9,] 
fits <- lapply(dat2[,-1], function(x) lm(x ~ dat2[,1] - 1)) 

ggplot2とデータと回帰ラインをプロットします。 reshape2パッケージのmeltを使用すると、データを長くすることができます。これは、カテゴリ内のデータを壊すggplot2役立ちます

library(ggplot2); library(reshape2) 
mdat <- melt(dat2,id="x") 
ggplot(data = mdat) + 
    geom_point(aes(y = value, x = x, color=variable)) + 
    geom_smooth(data=mdat,aes(x=x,y=value,color = variable),formula = y ~ x - 1, 
       method = "lm", se = F)