0

我正在尝试将回归线拟合到这个关系 angc~ext。变量 pch 将数据分为两组，我想分别拟合回归线

及其置信区间。这是我的数据框（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 )

我有两个问题：

  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 回答