特定の地域の数年間の毎日の降雨データがあります。極端な降雨イベントの洞察を得るために、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")
