Sou novo na página e bastante novo em estatística e R. Estou trabalhando em um projeto para a faculdade com o objetivo de encontrar a correlação entre a chuva e o nível do fluxo de água nos rios. Depois que a correlação for comprovada, quero prever / prever.
Os dados Eu tenho um conjunto de dados de vários anos (coletados a cada 5 minutos) para um determinado rio que contém:
- Precipitação em milímetros
- Fluxo do rio em metros cúbicos por segundo
Como o rio não tem neve, o modelo é baseado apenas na chuva e no tempo. Ocasionalmente, há temperaturas congelantes, mas estou pensando em remover esses períodos dos dados, de modo tão externo quanto essa situação está fora do escopo do meu projeto.
Exemplos Aqui você tem alguns gráficos de dados da chuva e do aumento da água algumas horas depois.
A linha vermelha é o fluxo do rio. A laranja é a chuva. Você pode ver que sempre chove antes que a água suba no rio. Há chuva começando novamente no final da série cronológica, mas isso afetará o fluxo do rio mais tarde.
A correlação está lá. Aqui está o que eu fiz no R para provar a correlação usando o ccf no R:
- a correlação cruzada
- a variável principal
- o atraso
Esta é minha linha R usada para o segundo exemplo (um período de precipitação):
ccf(arnoiaex1$Caudal, arnoiaex1$Precip, lag.max=1000, plot=TRUE, main="Flow & Rain")
Minha interpretação é:
- que a chuva leva (acontece primeiro),
ccf
Estou certo?
Sobre as séries temporais . Esta série cronológica não possui periodicidade ou sazonalidade. A chuva pode chegar a qualquer momento e causar um efeito. Reduz no verão, mas ainda acontece, é uma área com muita chuva o ano todo.
Modelo e previsão.
Não sei como criar um modelo para poder fazer uma previsão que me diga quanto um rio vai aumentar o volume após um período de chuva. Eu tenho tentado alguns arima
, auto arima
mas não tive muito sucesso. Devo usar Arima
, vars
ou outro modelo multivariado diferente? Qualquer link para um exemplo seria de grande ajuda.
Por favor, deixe-me saber se você sabe a melhor maneira de criar essa previsão, que modelo devo usar. Há algumas outras coisas que estou pensando em fazer, mas as tirei dessa explicação por simplicidade. Posso compartilhar alguns dados, se necessário.