分位数回归拆分图答案

【问题标题】：Quantile regression split plots分位数回归拆分图
【发布时间】：2018-02-20 15:08:43
【问题描述】：

我有一个特定地区几年的每日降雨量数据。为了深入了解极端降雨事件，我在 R 中使用了分位数回归（quantreg 包）。全天的图如下所示。我想要的是在中间（或其他点）分割回归线并分别拟合数据的前半部分和后半部分以查看差异。

这是我使用 quantreg 的方式：

plot(data$ahmAnn~data$Days, type="p", pch=20,cex=.4, col="gray50",
                                        xlab="Days", ylab="Rainfall")
qr <- abline(rq(data$ahmAnn~data$Days,tau=.99),col="red")

【问题讨论】：

标签： r regression quantile quantreg

【解决方案1】：

如果您想在不同的数据集上运行 quantreg::rq，请替换

data$ahmAnn~data$Days

与

x <- 10
stopifnot(x <= nrow(data))

set1 <- data[1:X,]
abline(rq(set1$ahmAnn~set1$Days,tau=.99),col="red")

set2 <- data[X:nrow(data),]
abline(rq(set2$ahmAnn~set2$Days,tau=.99),col="red")

【讨论】：

我认为我的问题不清楚，现在编辑它。我在寻找什么：将回归线分成 2 或 3 个不同的部分，并放在一个图中。
查看我编辑的答案。这是你要找的吗？
感谢您的回复。是的，它为数据的前半部分和后半部分提供了两条回归线，但有什么办法可以绘制它，第一条回归线在中间停止，第二条从中间开始？现在两条线都覆盖了整个散点图。
不使用abline，它将在整个情节中画一条线。您可以使用 fit <- rq(...) predict(fit, newdata) 在某个范围内生成值，其中 newdata 具有您要查找的范围。看到这个post