在 R 中可视化固定效应回归

问题描述

我有一个包含各种解释变量的面板数据模型,但我主要对其中一个的影响感兴趣。我想通过将所有其他解释变量置于它们的平均值来可视化我的模型,并将我的模型绘制在我的因变量和我的主要解释变量的散点图上。

如何在 R 中编写代码

非常感谢您的回答。

解决方法

你要找的是一个变量的边际效应,一种快速的方法是使用sjPlot,虽然我不确定“面板数据模型”,这是否意味着你使用的是线性模型lm()和面板回归 plm(),下面是 lm() 的示例:

library(sjPlot)
library(ggplot2)

fit <- lm(mpg ~ hp + disp + drat,data = mtcars)
plot_model(fit,type = "pred",terms = "hp")

enter image description here