问题描述
我正在使用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)