Linear RDD绘图仅显示两个数据点的问题求助

问题描述 #

运行RDD相关代码时，plot(reg_para0)仅显示两个数据点，同时收到警告：[1] "Mass points detected in the running variable."，已尝试dev.off()无效。

相关代码：

# 初始代码（输出符合预期）
data(house)  
house_rdd = rdd_data(x=x, y=y, data=house, cutpoint=0)  
summary(house_rdd) 
plot(house_rdd) 

# 后续代码（绘图异常）
reg_para0 = rdd_reg_lm(rdd_object=house_rdd, slope =("same"))
summary(reg_para0)
plot(reg_para0,
     ylim=c(0,1))

可行解决方案 #

1. 调整绘图的type参数 #

默认绘图逻辑在存在密集点时仅展示断点两侧均值，指定type="fit"可强制显示拟合线与原始数据：

plot(reg_para0,
     ylim=c(0,1),
     type="fit")

2. 处理运行变量的密集点 #

警告提示运行变量x存在大量重复值，可添加微小随机抖动分散数据，同时保留断点分组逻辑：

# 给x添加小幅度抖动
house$x_jitter = jitter(house$x, amount = 0.01)
# 重建RDD对象并重新拟合
house_rdd_jitter = rdd_data(x=x_jitter, y=y, data=house, cutpoint=0)
reg_para0_jitter = rdd_reg_lm(rdd_object=house_rdd_jitter, slope =("same"))
# 绘图
plot(reg_para0_jitter, ylim=c(0,1))

注：抖动幅度需根据x的取值范围调整，避免破坏断点两侧的分组。

3. 用ggplot2自定义绘图 #

如果默认绘图函数不满足需求，可使用ggplot2手动构建更灵活的可视化：

library(ggplot2)
# 生成拟合预测值
house$pred = predict(reg_para0, newdata=house)
# 绘制散点+拟合线+断点
ggplot(house, aes(x=x, y=y)) +
  geom_point(alpha=0.3) +
  geom_line(aes(y=pred, color=factor(x>=0)), linewidth=1) +
  geom_vline(xintercept=0, linetype="dashed") +
  ylim(0,1) +
  theme_minimal()

4. 对运行变量进行分箱处理 #

创建rdd_data时指定分箱参数，让绘图展示分箱均值与拟合线：

# 创建带分箱的RDD对象
house_rdd_binned = rdd_data(x=x, y=y, data=house, cutpoint=0, binning=TRUE, bins=50)
# 重新拟合并绘图
reg_para0_binned = rdd_reg_lm(rdd_object=house_rdd_binned, slope =("same"))
plot(reg_para0_binned, ylim=c(0,1))

注：bins参数可调整分箱数量，平衡展示效果与数据细节。

内容的提问来源于stack exchange，提问作者Andy_H