如何在rddtools rdd_reg_lm函数中使用协变量？

Question

我正在尝试使用rddtools R软件包运行参数RD回归。但是，包装文档对我来说不是很清楚。

First：定义RD对象的函数是：rdd_data(y, x, covar, cutpoint, z, labels, data)

帮助文件中的covar仅表示“外生变量”]。但是什么类型的呢？一个数据框？清单？

第二：函数rdd_reg_lm再次要求以这种方式通知协变量：

rdd_reg_lm(rdd_object, covariates = NULL, order = 1, bw = NULL,
  slope = c("separate", "same"), covar.opt = list(strategy = c("include",
  "residual"), slope = c("same", "separate"), bw = NULL),
  covar.strat = c("include", "residual"), weights)

在这里，根据帮助文件，协变量参数仅表示“包含协变量的公式”。同样，我不清楚应用这些协变量的正确方法是什么。

此外，是否可以在此函数rdd_data()和rdd_reg_lm()中包含多个协变量？

我非常感谢您的帮助。我已经读过一遍又一遍的帮助和小插图文件，在许多博客中搜索过，仍然一无所获。

我已经在下面检查了此主题

How to include a linear trend in a regression discontinuity design using rddtools

向我显示了以下示例：

rd.medic <- rdd_data(y = er ,x = ageyrs, covar = ageyrs, cutpoint=65, data = medicare)

rd.reg <- rdd_reg_lm(rdd_object=rd.medic, covariates = 'ageyrs', slope =
("same"), covar.opt = list("include"))

即使如此，语法仍然不清晰，因为我试图添加多个协变量而没有成功

谢谢！

Answer 1

您可以使用协变量创建一个数据框，然后将其包含在rdd_data中。

covariates<-data.frame(z1=ageyrs, z2=ageyrs2)

rd.medic <- rdd_data(y = er ,x = ageyrs, covar = covariates, cutpoint=65, data = medicare)

rd.reg <- rdd_reg_lm(rdd_object=rd.medic, covariates =TRUE, slope =("same"))