我正在预测每天前往医院进行MR扫描的客户数量。我有过去4年来每天到医院就诊的客户数量。但是我无法准确捕捉到不同月份客户到医院就诊数量的日变化。
我正在使用Rstudio,我已经尝试了Rob Hyndman建议的arima
。
modelfitsample<- read.csv("data_xreg_train.csv")
modeltest <- read.csv("data_xreg_test.csv")
ts_beverly_train <- ts(modelfitsample$Volume, start = c(2015,1), frequency=365.25)
ts_beverly_test <- ts(modeltest$Volume, start = c(2018,1), frequency=365)
xreg <- cbind(month=model.matrix(~as.factor(modelfitsample$Month)))
xreg1 <- cbind(month=model.matrix(~as.factor(modeltest$Month)))
modArima <- auto.arima(ts_beverly_train, xreg=xreg)
modArima
fit11 <- forecast(modArima, h=485, xreg = xreg1)
plot(fit11)
我需要一个能够捕捉日常变化并考虑到月度季节性的预测。