Estou estudando soluções para realizar uma predição de um produto que depende de outras variáveis.
Nesse meu estudo estou usando como minha base de dados Seatbelts
, uma Série Temporal que já é nativa do R. Que se trata de uma série histórica de acidentes de carros com morte de 1969 a 1983.
Nessa base possui oito variáveis, sendo que meu objetivo é usar a variável DriversKilled
(Motoristas Mortos) para criar um modelo preditivo que preveja a quantidade de Motoristas Mortos nos próximos 5 anos.
colnames(Seatbelts)
[1] "DriversKilled" "drivers" "front" "rear" "kms"
[6] "PetrolPrice" "VanKilled" "law"
Usando o modelo de Regressão Linear tslm
e incluindo as variáveis trend
(tendência) e season
(sazonalidade), consegui fazer a previsão usando a função Forecast
com sucesso.
mortos = window(Seatbelts[,c("DriversKilled")], start = c(1975,1), end= c(1984,12))
treino = window(mortos, start=c(1975,1), end=c(1979,12))
teste = window(mortos, start=c(1980,1), end=c(1984,12))
modelo_1 = tslm(treino ~ trend + season, data = treino)
Prev1 = forecast(modelo_1, h = 60)
plot(mortos)
lines(Prev1$mean, col="red")
Como visto na imagem acima, a cor vermelha indica a previsão usando a Regressão Linear, porém quero melhorá-la, considerando outras variáveis como por exemplo a variável drivers
(motoristas), com toda certeza esta variável influencia no número de mortos e se ela for incluída em meu modelo posso melhorar a acurácia da minha previsão.
É aí que começa meu problema, venho tentando acrescentar outras variáveis no meu modelo de Regressão Linear, mas não consigo fazer que a função Forecast
reconheça essas novas variáveis. Como posso fazer isso? É possível?
Como posso realizar uma predição usando um modelo de Regressão Linear Múltipla ou Multivariada no R?