控制 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)