0

この関係angc〜extに回帰直線を当てはめようとしています。変数pchは、データを2つのセットに分割し、それぞれに回帰直線を当てはめます。

その信頼区間で。これが私のデータフレーム（C）です：

"ext" "angc" "pch"
25 3.76288002820208 0
29 4.44255895177431 0
21 2.45214044383301 0
35 4.01334352881766 0
35 9.86225452423762 0
28 19.9304126868056 1
32 25.6984064030981 1
20 5.10582966112880 0
36 5.75603291081328 0
11 4.62311785943305 0
33 4.94401591414043 0
27 8.10039123328465 0
29 16.3882499757369 1
30 29.3492784626796 1
29 3.85960848290140 0
32 5.35857680326963 0
26 4.86451443776053 0
16 8.22008387344697 0
30 10.2212259432413 0
32 17.2519440101067 1
29 27.5011256290209 1

私のコード：

c0 <-  C[C$pch == 0, ]
c1 <- C[C$pch == 1, ]
prd0 <- as.data.frame( predict( lm(c0$angc ~ c0$ext), interval = c("confidence") ) )
prd1 <- as.data.frame( predict( lm(c1$angc ~ c1$ext), interval = c("confidence") ) )
dev.new()
plot( C$angc ~ C$ext, type = 'n' )
points( c0$angc ~ c0$ext, pch = 17 ) # triangles
abline(lm(c0$angc ~ c0$ext)) # regression line
lines(prd0$lwr) # lower CI
lines(prd0$upr) # upper CI
points( c1$angc ~ c1$ext, pch = 1 ) # circles
abline(lm(c1$angc ~ c1$ext))
lines(prd1$lwr, type = 'l', lty = 3 )
lines(prd1$upr, type = 'l', lty = 3 )

私には2つの問題があります：

  1. 円に必要な回帰直線を取得するにはどうすればよいですか？ほぼ垂直な線である必要があります（c1を確認してください）
  2. 正しい信頼区間が得られません

ご協力ありがとうございました、

サンティ

4

1 に答える 1

2

ggplot2あなたはこれをかなり効率的に行うことができます：

ggplot(C, aes(x = ext, y = angc, shape = pch)) + geom_point() + 
     geom_smooth(method = "lm")

これにより、 vsの散布図（geom_point()）が作成されます。ここで、点の形状はに基づいています。さらに、の一意の要素ごとに回帰直線がプロットに描画されます。この名前は、データの平滑化されたバージョン、この場合は線形回帰を描画するという事実に由来しています。angcextpchpchgeom_smooth()

于 2013-03-11T10:21:56.433 に答える