使用 R 可视化多元线性回归答案

【问题标题】：Visualizing multiple linear regression using R使用 R 可视化多元线性回归
【发布时间】：2021-06-17 01:37:08
【问题描述】：

我正在尝试使用四个变量来拟合和可视化多元线性回归模型。

fit <- lm(general_data$Bacterial.Contigs ~ general_data$Soil.Temperature + general_data$Time.point + general_data$Year)

这里Time.point和Year是分类变量，其他是数值变量。

我使用以下代码创建了一个 3D 图。

library(plotly)

plot_ly(data = general_data, z = ~Bacterial.Contigs, x = ~Soil.Temperature, y = ~Time.point, color = ~Year, colors = c('#0C4B8E' ,'#BF382A'),opacity = 0.5) %>%
  add_markers( marker = list(size = 4))

情节是这样的：

如何在此图中为“拟合”模型添加回归线。我真的很感激任何帮助。谢谢

【问题讨论】：

This 可能会有所帮助。

标签： r data-visualization linear-regression r-plotly olsmultiplelinearregression

【解决方案1】：

您只用additive effects 拟合了一个模型，这意味着您的分类值只会增加或减少您的响应变量，不同类别的斜率不会改变。在 3D 绘图上很难将其可视化，建议您尝试 ggplot2。

以mtcars 为例，您基本上将拟合值放回数据框中并为拟合值调用一行：

dat = mtcars
dat$am = factor(dat$am)
dat$vs = factor(dat$vs)

fit <- lm(mpg ~ disp + am + vs,data=dat)
dat$fitted= fitted(fit)

library(ggplo2)
g = ggplot(dat) + geom_point(aes(x=disp,y=mpg)) + 
geom_line(aes(x=disp,y = fitted)) + facet_grid(am~vs)
print(g)

或者如果你需要一个情节：

library(plotly)
ggplotly(g)

【讨论】：