Question

我正在Card et al Medicare data上实现RDDtools包。我已经能够建立一个简单的参数回归。我现在要做的是在ageyrs中包含一个线性趋势。我输入ageyrs作为协变量，但返回错误：

library(rddtools)
load(file ="C:/Users/perdue/Desktop/Adv.MicroEconometrics/HA 9/medicare.Rdata")

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

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

    Error: $ operator is invalid for atomic vectors
4.
match.arg(covar.opt$strategy, choices = c("include", "residual"))
3.
model.matrix.rdd_data(rdd_object, covariates = covariates, order = order, bw = bw, slope = slope, covar.opt = covar.opt)
2.
model.matrix(rdd_object, covariates = covariates, order = order, bw = bw, slope = slope, covar.opt = covar.opt)
1.
rdd_reg_lm(rdd_object = rd.medic, covariates = rd.medic$ageyrs, slope = ("same"), covar.opt = ("include"))

在documentation中它说我应该作为公式输入协变量。有没有人知道那个公式是什么？

Answer 1

问题可能是covar.opt必须是一个列表。这没有产生任何错误：

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

如何使用rddtools在回归不连续性设计中包含线性趋势

1 个答案: