分位数回归分裂图

Question

我有几年特定地区的每日降雨量数据。为了深入了解极端降雨事件，我在 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")

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

分位数回归分裂图

Quantile regression split plots

r

regression

quantile

quantreg