绘制ggplot2中的多元回归线

时间:2017-02-14 03:58:24

标签: r ggplot2 linear-regression

这里有点新手问题。我想用ggplot2来绘制我的多元回归模型,但是我遇到了一个错误。在决定询问社区之前,我在论坛上做了一些研究。 original question错误地发布在错误的论坛中。

我理解如何绘制简单的线性回归:

ggplot(data, aes(x=X, y=Y))+geom_point()+ 
  geom_smooth(method='lm',formula=Y~X)

但是当我尝试进行多元回归时:

Model<-lm(Y~x1*x2*x3*x4*x5, data, na.action=na.omit)

ggplot(data, aes(x=X, y=Model))+geom_point()+ 
  geom_smooth(method='lm',formula=Y~x1*x2*x3*x4*x5)

我收到此错误:

  

不知道如何自动为lm类型的对象选择比例。违约持续。
      错误:美学必须是长度1或与数据(209)相同:x,y

2 个答案:

答案 0 :(得分:0)

您是否尝试过使用abline()函数?

尝试使用

abline(Model,lwd=3,col="red")

答案 1 :(得分:0)

如果你有两个独立的变量,你可以尝试多次回归

set.seed(1) 
data <- data.frame(x1=rnorm(100, 5, 2), x2=rnorm(100, 10, 5)) # some random data
data$y <- 3*data$x1 + 2*data$x2 + rnorm(100, 2, 5)

lmFit <- lm(y ~ x1 + x2, data = data)

x1 <- seq(from = min(data$x1), to =  max(data$x1), by = 0.1)
x2 <- seq(from = min(data$x2), to = max(data$x2), by = 0.1)
y <- outer(X = x1, Y = x2, FUN = function(x1, x2) {
  predict(lmFit, newdata = data.frame(x1 = x1, x2 = x2))
})

open3d()
bg3d("white")
material3d(col="black")
persp3d(x1, x2, y,  col = "green3", xlab = "x1", ylab = "x2", zlab = "y",theta=50, phi=25, expand=0.75, ticktype="detailed")
points3d(data$x1, data$x2, data$y, col = "red", size=5)

enter image description here