用rdrobust软件包拟合具有协变量的R中的回归不连续模型

如何使用R中的“rdrobust”包为回归不连续模型指定协变量

问题在于下面代码中的“covs”参数。它应该允许为模型指定额外的协变量。我尝试过从数据集中提供一个向量（如下面的代码所示）和一个变量名向量，但两者都不起作用

fm = rdrobust(datasets$CONTINUANCE,
              dataset$T,
              p=3,
              covs=c(dataset$CONTENT_LENGTH, dataset$CONTENT_SNIPPET),
              kernel='uniform',
              weights=dataset$USERS,
              bwselect='msetwo',
              vce='nn',
              nnmatch=7,
              level=95)

因此，我得到以下错误：

Error in if (c <= x_min | c >= x_max) { : 
  missing value where TRUE/FALSE needed
Calls: rdrobust
Execution halted

如果（c=x_max）{ 缺少需要TRUE/FALSE的值 呼叫：rdrobust 停止执行

---

不幸的是，此包的文档没有帮助。您将协变量指定为公式（

var1+var2

）而不是向量（

c（var1，var2）

），因此类似的方法应该可以工作：

fm=rdrobust（数据集$CONTINUANCE，
数据集$T，
p=3，
covs=数据集$CONTENT\u LENGTH+数据集$CONTENT\u片段，
kernel='uniform'，
权重=数据集$USERS，
bwselect='msetwo'，
vce='nn'，
nnmatch=7，
级别=95）

在当前版本的rdrobust（0.99.4）中，应通过cbind将协变量合并为组合对象。在上述示例中：

fm = rdrobust(datasets$CONTINUANCE,
              dataset$T,
              p=3,
              covs=cbind(dataset$CONTENT_LENGTH, dataset$CONTENT_SNIPPET),
              kernel='uniform',
              weights=dataset$USERS,
              bwselect='msetwo',
              vce='nn',
              nnmatch=7,
              level=95)

我们将在即将发布的下一个版本中记住这一点。 谢谢

---

Sebastian。

谢谢！现在可以使用了。也许我会给软件包的作者发一封电子邮件来更新文档。谢谢Sebastian-非常有用！