我想使用" MA"中的lmodel2()
函数计算线性回归。 (长轴)方法(使用lmodel2
包)。另外,我想获得回归的斜率,其中我明确地将截距给lmodel2()
作为0.(因为数据不能取负值,我知道截距应为0的先验)。使用函数lm()
,以下语法通过添加 0 + :
lin <- data.frame(x = c(0:6), y = c(0.3, 0.1, 0.9, 3.1, 5, 4.9, 6.2))
plot (lin$x, lin$y)
fit <- lm(lin$y ~ 0 + lin$x)
summary(fit)
如果我尝试使用与lmodel2函数相同的方法(添加 0 + ),它就不起作用了。两者都给出相同的结果并估计非零截距。
library(lmodel2)
fit_1 <- lmodel2(lin$y ~ lin$x)
fit_2 <- lmodel2(lin$y ~ 0 + lin$x)
此外,如果可能,了解如何使用abline()
绘制结果将是有价值的。
(注意这篇文章使用帖子中的示例代码&#34;线性回归与R&#34中的显式截距;这解释了如何使用普通最小二乘回归抑制截距lm()
)