分位数回归分裂图

Question

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

以下是我使用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")

Answer 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")