强制固定效应的总和为零

问题描述

我正在运行具有固定效应 (id) 的回归模型。我不想使用默认值，而是将固定效果的总和强制为零。换句话说，下面最后一行代码的结果应该为零。我相信有一种方法可以使用对比参数来实现这一点，但我一直无法确切地找出方法。有什么想法吗？

library(lfe)
m1 <- felm(y ~
x +
z | id,data = df)
fx_m1 <- getfe(m1)
sum(fx_m1$effects)

解决方法

如何scale处理数据？ 示例：

library(plm)
data("Produc",package="plm")
Produc[-(1:3)] <- scale(Produc[-(1:3)])
zz <- plm(gsp ~ pcap + pc + emp + unemp,data=Produc,index=c("state","year"),model="within")
round(sum(fixef(zz)),14)
# [1] 0