Question

我正在尝试将回归线拟合到这个关系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 )

我有两个问题：

如何为圆圈获得所需的回归线？它应该是一条几乎垂直的线（检查c1）
我没有得到正确的置信区间

感谢您的帮助，

桑蒂

Answer 1

在ggplot2中，您可以相当有效地执行此操作：

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

这将创建geom_point()与angc的散点图（ext），其中点的形状基于pch。此外，在pch中的每个唯一元素的绘图中绘制了一条回归线。名称geom_smooth()来自于它绘制了数据的平滑版本，在这种情况下是线性回归。

将多个回归线拟合到R中的散点图

1 个答案: