Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- data1errors<-function(model)
- {
- mse <- 0
- mae <- 0
- mape <- 0
- amape <- 0
- forecast = forecast.Arima(model, 5)
- forecast.check <- cbind(data1.xts$index[846:850],forecast$mean)
- for (i in 1:5) {
- mse = mse + (forecast.check[i,1] - forecast.check[i,2])^2
- mae = mae + abs(forecast.check[i,1] - forecast.check[i,2])
- mape = mape + abs((forecast.check[i,1] - forecast.check[i,2])/forecast.check[i,1])
- amape = abs((forecast.check[i,1] - forecast.check[i,2])/(forecast.check[i,1]+forecast.check[i,2]))
- }
- mse = mse/5
- mae = mae/5
- mape = mape/5
- amape = amape/5
- fore_check1 = cbind(mse,mae,mape,amape)
- return(fore_check1)
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement