У меня есть данные временного ряда для количества осадков с 1974 по 2018 год, однако я не смог найти правильное распределение, чтобы описать его, так как оно сильно искажено с более чем 90% ежедневного количества осадков равным 0 и случайным штормовым событием,Кроме того, существует тенденция, что каждые 10 лет будут происходить большие ливни.Что я должен сделать, чтобы найти правильное распределение в соответствии с шаблоном осадков?
Я попытался использовать пакет fitdistr, нашел гамма-распределение с, казалось бы, хорошо подходящим cdf, однако он сильно недооценил количество осадков.
library("fitdistrplus")
rain <- Rainfall$Rainfall
rain
fitdistr(rain,"normal")
plotdist(rain)
descdist(rain, discrete = FALSE)
fit.norm <- fitdist(rain, "norm")
plot(fit.norm)
fit.1 <- fitdist(rain, "gamma", method = "mle", lower = c(0,0), start =
list(scale = 1, shape = 1))
plot(fit.1)
фактические результаты значительно меньше ожидаемых результатов