Как включить линейный тренд в проект с разрывом регрессии с помощью rddtools

Question

Я внедряю пакет RDDtools для данных Card et al Medicare. Мне удалось установить простую параметрическую регрессию. То, что я хотел бы сделать сейчас, это включить линейный тренд в 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"))

В документации говорится, что я должен вводить ковариаты в виде формулы. У кого-нибудь есть представление о том, что это за формула?

Answer 1

Проблема может быть в том, что covar.opt должен быть список. Это не произвело никакой ошибки:

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