0

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

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

其中 covar 在帮助文件中仅表示"Exogeneous variables"。但是什么类型？数据框？一个列表？

第二：函数 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)

其中，根据帮助文件， covariates 参数仅表示"Formula to include covariates"。同样，我不清楚应用这些协变量的正确方法是什么。

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

我很感激这里的一些帮助。我已经一次又一次地阅读了帮助和小插图文件，在许多博客中进行了搜索，但仍然一无所获。

我已经在下面检查过这个主题

如何使用 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"))

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

谢谢！

4

1 回答 1

0

您可以使用协变量创建数据框，然后将其包含在 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"))
于 2020-05-19T16:28:17.763 回答