我尝试运行以下代码:
data(house)
house_rdd = rdd_data(x=x, y=y, data=house, cutpoint=0)
summary(house_rdd)
plot(house_rdd)
当我绘制它时,我得到this,这是有道理的。
reg_para0 = rdd_reg_lm(rdd_object=house_rdd, slope =("same"))
summary(reg_para0)
plot(reg_para0,
ylim=c(0,1))
但是,当绘制回归图时,它仅显示两个观察结果,如此处
所示此外,我收到警告消息:“[1]“在运行变量中检测到质量点。””
图表应该看起来像this
有人知道这个问题的潜在解决方案吗?我尝试过“dev.off()”,但没有成功。
谢谢
原因是第二次调用
plot()
是创建一个全新的绘图,而不是覆盖在前一个绘图上。
这应该会让您接近您想要的(您需要清理 x 轴标题)。
library(rddtools)
data(house)
house_rdd = rdd_data(x=x, y=y, data=house, cutpoint=0)
plot(house_rdd, ylim=c(0,1))
par(new=TRUE)
reg_para0 = rdd_reg_lm(rdd_object=house_rdd, slope =("same"))
plot(reg_para0, ylim=c(0,1), axes = FALSE, xlab = "", ylab = "")
par(new=TRUE)
阻止创建新图。