R中回归不连续的固定效应的控制

问题描述

我正在使用Rdrobust软件包来评估国家政策对县级成果的影响。在我的协变量中,我包括了指示状态的虚拟变量,以控制状态级别的固定效果。但是,当我们运行代码时,出现以下错误消息:

chol.default(ZWZ)错误

33号命令的未成年人不是肯定的。

其中Z是包含我的钴铁矿的矩阵,第33个变量是状态1的虚拟。

我的代码是:

out = rdrobust(y,x,covs=z,kernel = "triangular",p=2,bwselect="mserd",cluster= cluster)

我不确定如何清除此错误消息。更重要的是,我想知道Rdrobust软件包中是否还有其他方法可以控制固定效果(例如状态级别)。

预先感谢您的帮助

解决方法

最近的 rdrobust 版本 1.0.1 应该处理这个错误。要从分类变量“状态”生成固定效果,您可以执行以下操作:

state.f = factor(state)

state.d = model.matrix(~state.f+0)

out = rdrobust(y,x,covs = state.d)