11. Comparison

gör regression forecast error vector regression_forecast_error <- testing$adj.close - regression_forecast

gör regression table regression_table <- data.frame(cbind(testing$adj.close, regression_forecast, regression_forecast_error))

ändra namnen i regression table names(regression_table) <- c("testing", "estimated", "residuals")

visa regression table regression_table

gör testing values vector testing_values <- testing$adj.close

gör holt forecast vector holt_forecast <- holt_model$mean

gör holt forecast error vector holt_forecast_error <- testing_values - holt_forecast

gör holt table holt_table <- data.frame(cbind(testing_values, holt_forecast, holt_forecast_error))

ändra namnen i holt table names(holt_table) <- c("testing", "estimated", "residuals")

visa holt table holt_table

gör forecast111 vector av ARIMA111 (h kan vara andra siffror och det kan va annan ARIMA) forecast111 <- forecast(ARIMA111, h=4)

gör ARIMA111 forecast table ARIMA111_forecast <- exp(forecast111$mean)

gör ARIMA111 forecast error vector ARIMA111_forecast_error <- testing$adj.close - ARIMA111_forecast

gör ARIMA111 table ARIMA111_table <- data.frame(cbind(testing$adj.close, ARIMA111_forecast, ARIMA111_forecast_error))

byt namn i ARIMA111 table names(ARIMA111_table) <- c("testing", "estimated", "residuals")

visa ARIMA111 table ARIMA111_table

gör residuals table Residuals_table <- data.frame(cbind(regression_table$residuals, holt_table$residuals, ARIMA111_table$residuals))

ändra namnen på kolumner i residuals table names(Residuals_table) <- c("regression", "holt", "ARIMA111")

visa residuals table Residuals_table

gör RMSE regression vector RMSE_regression <- rmse(regression_table$testing, regression_table$estimated)

gör RMSE holt vector RMSE_holt <- rmse(holt_table$testing, holt_table$estimated)

gör RMSE ARIMA111 vector (det måste inte vara 111 hela tiden) RMSE_ARIMA111 <- rmse(ARIMA111_table$testing, ARIMA111$estimated)

gör RMSE table RMSE_table <- data.frame(cbind(RMSE_regression, RMSE_holt, RMSE_ARIMA111), row.names="RMSE")

byt namn på kolumnerna i RMSE table names(RMSE_table) <- c("regression", "holt", "ARIMA111")

visa RMSE table RMSE_table

Click