Hello! I am having trouble including covariates in a fuzzy RDD model in R with rddtools. I run the following commands:
library(rddtools) data.complete2 <- as.data.frame(cbind(GDPpc.rel, GR.rate, fitval.firstStage)) colnames(data.complete2) <- c("GDPpc.rel", "GR.rate", "FitTreat") data.RDD<-rdd_data(x=data.complete2$GDPpc.rel, y=data.complete2$GR.rate, z=data.complete2$FitTreat, covar=dummytime, cutpoint = 0.75) data.RDD1<-rdd_data(x=data.complete2$GDPpc.rel, y=data.complete2$GR.rate, z=data.complete2$FitTreat, covar=cbind(dummytime, dummyreg), cutpoint = 0.75) reg_para <- rdd_reg_lm(rdd_object=data.RDD, covar.opt = list(strategy= "include"), order = 4) reg_para1 <- rdd_reg_lm(rdd_object=data.RDD1, covar.opt = list(strategy = "include"), order = 4) where dummytime and dummyreg are two dummy variables considering time and regional effects. The problem is that looking at the summary of the two regressions, they seem to be exactly the same and the coefficient which should refer to the dummy variables is not reported. My guess is that I did not actually implement the covariates in the regression, but I cannot understand where is the problem. Any suggestion? Thank you in advance. Best regards, Chiara Malavasi [[alternative HTML version deleted]] ______________________________________________ Remail@example.com mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.