ViroReportR (version 1.0.4)

generate_validation_metric: Compute Forecast Validation Metrics (SMAPE & MASE)

Description

This function evaluates forecast accuracy across multiple validation runs by computing two key performance metrics:

Usage

generate_validation_metric(data, validation_res)

Value

A tibble (data frame) with one row per forecast result and the following columns:

• train_period: Date range of the training period used for the forecast.

• forecast_period: Date range of the forecasted period.

• smape: Symmetric Mean Absolute Percentage Error between forecasted and actual values, rounded to two decimals.

• mase: Mean Absolute Scaled Error, rounded to two decimals.

Arguments

data

A data frame used in generate_validation(), containing the original training data for the model. It must include:

• date: Dates of the observed case data (class Date).

• confirm: Numeric values of observed confirmed cases.

validation_res

A list of forecast validation results, typically the output from generate_validation(). Each element should contain:

• forecast_res_quantiles: A data frame with columns date and p50 (median forecasted values).

• original_data: A data frame representing the training data used for that forecast, with a date column.

Details

• Symmetric Mean Absolute Percentage Error (SMAPE): Measures relative forecast accuracy while remaining robust to zero values in the actual data.

• Mean Absolute Scaled Error (MASE): Scales forecast errors relative to the in-sample one-step naïve forecast, allowing comparison across series with different scales.

For each forecast result, the function also reports the corresponding training and forecast periods. Computation stops once the forecast period reaches the maximum date in the model data.

• SMAPE is defined as: $$SMAPE = mean( |F - A| / ((|A| + |F|) / 2) )$$ where \(A\) are actual values and \(F\) are forecasts. It avoids division by zero and is suitable for count data with zeros.

• MASE compares the mean absolute forecast error against the mean absolute difference of successive actual: $$MASE = mean(|A - F|) / mean(|diff(A)|)$$

The function automatically excludes forecasts extending beyond the latest date in the observed model data.

See Also

generate_validation(), generate_forecast()

Examples

data <- simulate_data()
formatted_data <- get_aggregated_data(data,"date", "flu_a", "2024-10-16", "2024-12-31")
start_date <- ("2024-10-16")
validation_results <- generate_validation(formatted_data, start_date, type ="flu_a")
generate_validation_metric(formatted_data, validation_results)