일변량 모델로 단일 시계열 예측

데이터 세트 생성

ML 모델을 저장할 BigQuery 데이터 세트를 만듭니다.

1. Google Cloud 콘솔에서 BigQuery 페이지로 이동합니다.

   BigQuery 페이지로 이동

2. 탐색기 창에서 프로젝트 이름을 클릭합니다.

3. more_vert 작업 보기 > 데이터 세트 만들기를 클릭합니다.

   

4. 데이터 세트 만들기 페이지에서 다음을 수행합니다.

   • 데이터 세트 ID에 bqml_tutorial를 입력합니다.

   • 위치 유형에 대해 멀티 리전을 선택한 다음 US(미국 내 여러 리전)를 선택합니다.

     공개 데이터 세트는 US 멀티 리전에 저장됩니다. 편의상 같은 위치에 데이터 세트를 저장합니다.

   • 나머지 기본 설정은 그대로 두고 데이터 세트 만들기를 클릭합니다.

     

입력 데이터 시각화

모델을 만들기 전에 원하는 경우 입력 시계열 데이터를 시각화하여 분포를 파악할 수 있습니다. Looker Studio를 사용하여 이를 수행할 수 있습니다.

시계열 데이터를 시각화하려면 다음 단계를 따르세요.

import bigframes.pandas as bpd

# Start by loading the historical data from BigQuerythat you want to analyze and forecast.
# This clause indicates that you are querying the ga_sessions_* tables in the google_analytics_sample dataset.
# Read and visualize the time series you want to forecast.
df = bpd.read_gbq("bigquery-public-data.google_analytics_sample.ga_sessions_*")
parsed_date = bpd.to_datetime(df.date, format="%Y%m%d", utc=True)
visits = df["totals"].struct.field("visits")
total_visits = visits.groupby(parsed_date).sum()

# Expected output: total_visits.head()
# date
# 2016-08-01 00:00:00+00:00    1711
# 2016-08-02 00:00:00+00:00    2140
# 2016-08-03 00:00:00+00:00    2890
# 2016-08-04 00:00:00+00:00    3161
# 2016-08-05 00:00:00+00:00    2702
# Name: visits, dtype: Int64

total_visits.plot.line()

시계열 모델 만들기

totals.visits 열로 표시되는 총 사이트 방문을 예측하는 시계열 모델을 만들고 Google 애널리틱스 360 데이터에서 학습합니다.

다음 쿼리에서 OPTIONS(model_type='ARIMA_PLUS', time_series_timestamp_col='date', ...) 절은 ARIMA 기반 시계열 모델을 만들고 있음을 나타냅니다. CREATE MODEL 문에 있는 auto_arima 옵션은 기본적으로 TRUE이므로 auto.ARIMA 알고리즘이 모델에서 초매개변수를 자동으로 조정합니다. 이 알고리즘은 후보 모델 십여 개를 접합하고 Akaike 정보 기준 (AIC)이 가장 낮은 최적 모델을 선택합니다. CREATE MODEL 문에서의 data_frequency 옵션은 기본적으로 AUTO_FREQUENCY이므로 학습 프로세스가 입력 시계열의 데이터 빈도를 자동으로 추론합니다. CREATE MODEL 문에 있는 decompose_time_series 옵션은 기본적으로 TRUE로 설정되므로 다음 단계에서 모델을 평가할 때 시계열 데이터에 관한 정보가 반환됩니다.

모델을 만들려면 다음 단계를 따르세요.

1. Google Cloud 콘솔에서 BigQuery 페이지로 이동합니다.

   BigQuery로 이동

2. 쿼리 편집기에 다음 쿼리를 붙여넣고 실행을 클릭합니다.

   CREATE OR REPLACE MODEL `bqml_tutorial.ga_arima_model`
   OPTIONS
   (model_type = 'ARIMA_PLUS',
     time_series_timestamp_col = 'parsed_date',
     time_series_data_col = 'total_visits',
     auto_arima = TRUE,
     data_frequency = 'AUTO_FREQUENCY',
     decompose_time_series = TRUE
   ) AS
   SELECT
   PARSE_TIMESTAMP("%Y%m%d", date) AS parsed_date,
   SUM(totals.visits) AS total_visits
   FROM
   `bigquery-public-data.google_analytics_sample.ga_sessions_*`
   GROUP BY date;

   쿼리가 완료되는 데 약 4초가 걸리며 그 이후에는 ga_arima_model 모델이 탐색기 창에 표시됩니다. 이 쿼리에서는 CREATE MODEL 문을 사용하여 모델을 만들므로 쿼리 결과가 표시되지 않습니다.

후보 모델 평가

ML.ARIMA_EVALUATE 함수를 사용하여 시계열 모델을 평가합니다. ML.ARIMA_EVALUATE 함수는 자동 초매개변수 조정 과정에서 평가된 모든 후보 모델의 평가 측정항목을 보여줍니다.

모델을 평가하려면 다음 단계를 따르세요.

1. Google Cloud 콘솔에서 BigQuery 페이지로 이동합니다.

   BigQuery로 이동

2. 쿼리 편집기에 다음 쿼리를 붙여넣고 실행을 클릭합니다.

   SELECT
    *
   FROM
    ML.ARIMA_EVALUATE(MODEL `bqml_tutorial.ga_arima_model`);

   결과는 다음과 비슷하게 표시됩니다.

   

   non_seasonal_p, non_seasonal_d, non_seasonal_q, has_drift 출력 열은 학습 파이프라인에서 ARIMA 모델을 정의합니다. log_likelihood, AIC, variance 출력 열은 ARIMA 모델 접합 프로세스와 관련이 있습니다.

   auto.ARIMA 알고리즘은 KPSS 테스트를 사용하여 non_seasonal_d의 최적 값을 결정합니다. 이 경우에는 1입니다. non_seasonal_d가 1이면 auto.ARIMA 알고리즘이 42개의 서로 다른 후보 ARIMA 모델을 병렬로 학습시킵니다. 이 예시에서는 42개 후보 모델이 모두 유효하므로 출력에는 후보 ARIMA 모델마다 하나씩 42개의 행이 포함됩니다. 일부 모델이 유효하지 않은 경우 출력에서 제외됩니다. 이러한 후보 모델은 AIC에 따라 오름차순으로 반환됩니다. 첫 번째 행의 모델은 AIC가 가장 낮으며 최적 모델로 간주됩니다. 최적 모델은 최종 모델로 저장되며 모델에서 ML.FORECAST와 같은 함수를 호출할 때 사용됩니다.

   seasonal_periods 열에는 시계열 데이터에서 식별된 계절 패턴에 관한 정보가 포함됩니다. 이 열은 ARIMA 모델링과 관계가 없으므로 모든 출력 행에서 동일한 값을 갖습니다. 주간 패턴을 보고하며, 이는 입력 데이터를 시각화하도록 선택한 경우 확인한 결과와 일치합니다.

   has_holiday_effect, has_spikes_and_dips, has_step_changes 열은 decompose_time_series=TRUE인 경우에만 채워집니다. 또한 이러한 열은 입력 시계열 데이터에 관한 정보를 반영하며 ARIMA 모델링과 관련이 없습니다. 또한 이러한 열은 모든 출력 행에서 동일한 값을 갖습니다.

   error_message 열에는 auto.ARIMA 피팅 프로세스 중에 발생한 오류가 표시됩니다. 오류의 원인 중 하나는 선택한 non_seasonal_p, non_seasonal_d, non_seasonal_q, has_drift 열에서 시계열을 안정화하지 못하는 경우일 수 있습니다. 모든 후보 모델의 오류 메시지를 검색하려면 모델을 만들 때 show_all_candidate_models 옵션을 TRUE로 설정합니다.

   출력 열에 관한 자세한 내용은 ML.ARIMA_EVALUATE 함수를 참고하세요.

모델 계수 검사

ML.ARIMA_COEFFICIENTS 함수를 사용하여 시계열 모델의 계수를 검사합니다.

모델의 계수를 가져오려면 다음 단계를 따르세요.

1. Google Cloud 콘솔에서 BigQuery 페이지로 이동합니다.

   BigQuery로 이동

2. 쿼리 편집기에 다음 쿼리를 붙여넣고 실행을 클릭합니다.

   SELECT
    *
   FROM
    ML.ARIMA_COEFFICIENTS(MODEL `bqml_tutorial.ga_arima_model`);

   결과는 다음과 비슷하게 표시됩니다.

   

   ar_coefficients 출력 열에는 ARIMA 모델의 자동 회귀 (AR) 부분의 모델 계수가 표시됩니다. 마찬가지로 ma_coefficients 출력 열에는 ARIMA 모델의 이동 평균 (MA) 부분의 모델 계수가 표시됩니다. 두 열 모두 길이가 각각 non_seasonal_p 및 non_seasonal_q에 해당하는 배열 값을 포함합니다. ML.ARIMA_EVALUATE 함수의 출력에서 최적 모델의 non_seasonal_p 값은 2이고 non_seasonal_q 값은 3입니다. 따라서 ML.ARIMA_COEFFICIENTS 출력에서 ar_coefficients 값은 2개 요소 배열이고 ma_coefficients 값은 3개 요소 배열입니다. intercept_or_drift 값은 ARIMA 모델의 상수 항입니다.

   출력 열에 관한 자세한 내용은 ML.ARIMA_COEFFICIENTS 함수를 참고하세요.

모델을 사용하여 데이터 예측

ML.FORECAST 함수를 사용하여 미래 시계열 값을 예측합니다.

다음 GoogleSQL 쿼리에서 STRUCT(30 AS horizon, 0.8 AS confidence_level) 절은 쿼리가 30개의 미래 시점을 예측하고 80% 신뢰 수준의 예측 구간을 생성함을 나타냅니다.

모델을 사용하여 데이터를 예측하려면 다음 단계를 따르세요.

1. Google Cloud 콘솔에서 BigQuery 페이지로 이동합니다.

   BigQuery로 이동

2. 쿼리 편집기에 다음 쿼리를 붙여넣고 실행을 클릭합니다.

   SELECT
    *
   FROM
    ML.FORECAST(MODEL `bqml_tutorial.ga_arima_model`,
                STRUCT(30 AS horizon, 0.8 AS confidence_level));

   결과는 다음과 비슷하게 표시됩니다.

   

   출력 행은 forecast_timestamp 열 값별로 시간 순서로 정렬됩니다. 시계열 예측에서 prediction_interval_lower_bound 및 prediction_interval_upper_bound 열 값으로 표시되는 예측 구간은 forecast_value 열 값만큼 중요합니다. forecast_value 값은 예측 구간의 중간 포인트입니다. 예측 구간은 standard_error 및 confidence_level 열 값에 따라 달라집니다.

   출력 열에 관한 자세한 내용은 ML.FORECAST 함수를 참고하세요.

예측 결과 설명

ML.EXPLAIN_FORECAST 함수를 사용하여 예측 데이터 외에 설명 가능성 측정항목을 가져올 수 있습니다. ML.EXPLAIN_FORECAST 함수는 미래 시계열 값을 예측하고 시계열의 모든 개별 구성요소를 반환합니다.

ML.FORECAST 함수와 마찬가지로 ML.EXPLAIN_FORECAST 함수에 사용된 STRUCT(30 AS horizon, 0.8 AS confidence_level) 절은 쿼리가 미래 시점 30개를 예측하고 80% 신뢰도로 예측 구간을 생성함을 나타냅니다.

모델의 결과를 설명하려면 다음 단계를 따르세요.

1. Google Cloud 콘솔에서 BigQuery 페이지로 이동합니다.

   BigQuery로 이동

2. 쿼리 편집기에 다음 쿼리를 붙여넣고 실행을 클릭합니다.

   SELECT
    *
   FROM
    ML.EXPLAIN_FORECAST(MODEL `bqml_tutorial.ga_arima_model`,
      STRUCT(30 AS horizon, 0.8 AS confidence_level));

   결과는 다음과 비슷하게 표시됩니다.

   

   출력 행은 time_series_timestamp 열 값을 기준으로 시간 순으로 정렬됩니다.

   출력 열에 관한 자세한 내용은 ML.EXPLAIN_FORECAST 함수를 참고하세요.

   결과를 시각화하려면 입력 데이터 시각화 섹션에 설명된 대로 Looker Studio를 사용하여 다음 열을 측정항목으로 사용하여 차트를 만들 수 있습니다.

   • time_series_data
   • prediction_interval_lower_bound
   • prediction_interval_upper_bound
   • trend
   • seasonal_period_weekly
   • step_changes