这里有点新手问题。我想用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
答案 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)