求R语言中平衡面板数据下含州固定效应的OLS虚拟变量模型命令

嘿，针对你的平衡面板数据估计含州固定效应的OLS虚拟变量模型需求，我整理了两个实用的R命令方案，都是面板分析里常用的：

如果你希望直观看到每个州的虚拟变量系数，可以直接把州变量转为因子类型，lm()会自动生成对应的虚拟变量（默认省略一个基准州）。如果需要控制年份固定效应，也可以一并加入：

# 假设你的数据集名为panel_data，因变量是y，核心自变量是x1、x2，州变量是state，年份变量是year
model_lm <- lm(y ~ x1 + x2 + factor(state) + factor(year), data = panel_data)
summary(model_lm)

小贴士：如果不需要年份固定效应，删掉+ factor(year)即可。

对于面板数据，plm包会比基础lm更高效，尤其是大样本场景，语法也更简洁，而且结果和虚拟变量法完全一致：
首先得安装并加载包：

install.packages("plm")
library(plm)

接着把数据转换为面板格式，再指定固定效应模型：

# 将普通数据框转为面板数据框，指定个体（州）和时间维度（年份）
pdata <- pdata.frame(panel_data, index = c("state", "year"))

# 仅估计州固定效应（个体固定效应）
model_plm <- plm(y ~ x1 + x2, data = pdata, model = "within", effect = "individual")

# 如果需要同时控制州和年份的双向固定效应，改成：
# model_plm <- plm(y ~ x1 + x2, data = pdata, model = "within", effect = "twoways")

# 查看模型结果
summary(model_plm)

补充：plm的within模型不会直接输出每个州的虚拟变量系数，但可以用fixef(model_plm)命令提取这些固定效应估计值。

内容的提问来源于stack exchange，提问作者QIA-man