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

Question

我正在使用 Rdrobust 包来估计国家政策对县级结果的影响。 在我的协变量中，我包含了指示状态的虚拟变量，以控制状态级别的固定效应。 但是，当我们运行代码时，我收到以下错误消息：

chol.default(ZWZ) 中的错误：

33 阶的前导次要不是正定的。

其中 Z 是包含我的协变量的矩阵，第 33 个变量是状态 1 的虚拟变量。

我的代码是：

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

我不确定如何摆脱此错误消息。 更重要的是，我想知道是否还有其他方法可以控制 Rdrobust 包中的固定效果（例如状态级）。

在此先感谢您的帮助

Answer 1

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

state.f = factor(state)

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

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