使用fixest::coefplot时keep与x参数报错的解决咨询

核心问题原因 #

你误用了x参数：该参数需要传入与模型全部系数数量一致的标签向量（包括被keep过滤掉的固定效应/虚拟变量），而不是仅保留的系数。这就是报错提示长度不匹配（1 vs 4）的原因。

解决方案（满足所有诉求） #

方法1：使用rename参数针对性重命名（推荐）

fixest::coefplot()提供了rename参数，接受命名向量，仅对指定系数重命名，无需处理其他固定效应/虚拟变量，完美适配你的需求。

步骤示例：

1. 先计算带聚类稳健标准误的回归结果（以lm+sandwich为例）：

library(lmtest)
library(sandwich)

# 拟合基础lm模型
m1 <- lm(Sepal.Width ~ Petal.Length + factor(Species), data = iris)
# 计算聚类稳健标准误（假设按Species聚类）
m1_robust_clustered <- coeftest(m1, vcov = vcovCL, cluster = ~Species)

2. 调用coefplot时结合keep和rename：

fixest::coefplot(m1_robust_clustered,
                 main = "", horiz = TRUE, 
                 keep = "Petal.Length",
                 rename = c("Petal.Length" = "Petal Length"))

方法2：手动构造匹配长度的标签向量（备选）

如果你坚持要用x参数，可以先提取全部系数名称，仅修改目标变量后传入：

# 获取模型全部系数名称
all_coef_names <- names(coef(m1))
# 替换目标变量名称
all_coef_names[all_coef_names == "Petal.Length"] <- "Petal Length"

# 调用coefplot
fixest::coefplot(m1_robust_clustered,
                 main = "", horiz = TRUE, 
                 keep = "Petal.Length",
                 x = all_coef_names)

对诉求的回应 #

1. 无需提前重命名所有变量：方法1的rename参数仅针对你需要的系数修改，完全不用处理50个国家固定效应的名称；
2. 无需切换到feols：上述方案均基于lm模型+稳健标准误结果，直接适配fixest::coefplot()；
3. 保留稳健聚类标准误：通过coeftest+vcovCL计算聚类稳健标准误，结果可直接传入fixest::coefplot()，无需更换其他包的绘图函数。

内容的提问来源于stack exchange，提问作者Ttytamaki