制約条件を持つ最適化問題 scipy

制約条件を持つ最適化問題

制約条件を持つ最適化問題を解くためには、scipy.optimizeモジュールを使用することが一般的です。

以下は、scipyを使用して制約条件のある最適化問題を解くサンプルコードです。

この例では、制約条件付きの2次元目的関数を最小化する問題を考えます。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import minimize

# 目的関数の定義
def objective(x):
    return x[0]**2 + x[1]**2

# 制約条件
def constraint(x):
    return x[0] + x[1] - 1

# 初期値
x0 = np.array([0.5, 0.5])

# 制約条件付き最適化
constraints = {'type': 'eq', 'fun': constraint}
result = minimize(objective, x0, constraints=constraints)

# 結果の表示
print("最適な解:", result.x)
print("最小値:", result.fun)

[実行結果]

グラフ化

次に、結果をグラフ化してみましょう。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
# 目的関数を等高線プロット
x = np.linspace(-2, 2, 400)
y = np.linspace(-2, 2, 400)
X, Y = np.meshgrid(x, y)
Z = X**2 + Y**2

plt.contour(X, Y, Z, levels=np.logspace(-1, 2, 10), cmap='viridis')

# 制約条件をプロット
plt.plot(x, 1 - x, '-r', label='x + y = 1')

# 最適解をプロット
plt.plot(result.x[0], result.x[1], 'ro', label='最適解')

plt.xlabel('x')
plt.ylabel('y')
plt.title('Constrained Optimization')
plt.legend()
plt.show()

このコードでは、2次元の目的関数を最小化しながら、制約条件 x + y = 1 を満たす最適解を求めています。

結果は等高線プロットとして表示され、最適解は赤い点で示されます。

なお、制約条件のある最適化問題は様々な種類があり、それぞれ異なるアプローチが必要です。

上記の例は一般的なサンプルであり、実際の問題に応じて適切なアルゴリズムや制約条件の設計が必要です。

[実行結果]

異常検出 scikit-learn

異常検出

異常検出の例として、シンプルなユースケースとして1次元のデータで異常を検出することを考えてみましょう。

以下は、正規分布に従うデータに対して異常を検出する例です。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
import numpy as np
import matplotlib.pyplot as plt
from sklearn.ensemble import IsolationForest

# データ生成
np.random.seed(42)
normal_data = np.random.normal(loc=0, scale=1, size=1000)
anomalies = np.array([4, 5, -3])  # 異常値の追加

# データと異常値を結合
data = np.concatenate([normal_data, anomalies])

# データのプロット
plt.figure(figsize=(10, 6))
plt.scatter(range(len(data)), data, marker='x', color='blue', label='Normal')
plt.scatter(range(len(normal_data), len(data)), anomalies, marker='o', color='red', label='Anomaly')
plt.xlabel('Sample Index')
plt.ylabel('Value')
plt.title('Anomaly Detection Example')
plt.legend()
plt.show()

# Isolation Forestモデルの構築と異常検出
model = IsolationForest(contamination=0.05, random_state=42)  # contaminationは異常値の割合
model.fit(data.reshape(-1, 1))

# データの異常スコアを予測
anomaly_scores = model.decision_function(data.reshape(-1, 1))

# 異常スコアのプロット
plt.figure(figsize=(10, 6))
plt.plot(anomaly_scores, marker='o', linestyle='-', color='green')
plt.xlabel('Sample Index')
plt.ylabel('Anomaly Score')
plt.title('Anomaly Scores')
plt.show()

このコードでは、正規分布に従うデータに3つの異常値を追加しています。

Isolation Forestアルゴリズムを使用して異常検出を行い、異常スコアをプロットしています。

異常スコアが閾値を超えるデータポイントは、異常として検出されたとみなすことができます。


なお、実際のデータではより多くの特徴量や高次元データを扱うことになるため、適切なデータの前処理やモデルの調整が必要です。

また、異常の定義や異常値の割合に応じてcontaminationパラメータを調整する必要があります。

ソースコード解説

以下にコードの各部分の詳細な説明を行います。

1. モジュールのインポート:

最初の数行では、必要なライブラリとモジュールをインポートしています。

  • numpy（npとしてエイリアス）: 数値計算用のライブラリ。
  • matplotlib.pyplot（pltとしてエイリアス）: データの可視化用のライブラリ。
  • IsolationForestクラス: Scikit-learnの異常検出アルゴリズムであるIsolation Forestを提供するクラス。

2. データの生成と準備:

  • normal_data: 正規分布に従うランダムなデータを生成しています。
    平均 (loc) は0、標準偏差 (scale) は1で、サイズは1000です。
  • anomalies: 異常値の配列を手動で定義しています。

3. データと異常値のプロット:

  • plt.scatter: データと異常値をプロットしています。
    青い「x」マーカーは正常なデータポイントを示し、赤い「o」マーカーは異常値を示します。
    X軸はサンプルのインデックス、Y軸は値です。

4. Isolation Forestモデルの構築と異常検出:

  • IsolationForestのインスタンスを作成しています。
    contaminationパラメータは異常値の割合を指定します。
    random_stateは乱数のシードを固定します。
  • model.fit: モデルをデータに適合させます。

5. データの異常スコアを予測:

  • model.decision_function: データポイントの異常スコアを計算します。
    異常スコアは、モデルがデータポイントを異常とみなす確信度の尺度です。

6. 異常スコアのプロット:

  • plt.plot: データポイントの異常スコアをプロットしています。
    異常スコアが低いほど正常、高いほど異常と判断される可能性が高くなります。
    X軸はサンプルのインデックス、Y軸は異常スコアです。

このコードは、Isolation Forestアルゴリズムを使用して異常検出を行い、異常スコアをプロットして視覚化するサンプルです。

ただし、実際のケースではデータの前処理、ハイパーパラメータの調整、適切なモデルの選択などが必要です。

結果解説

異常検出の結果として表示される2つのグラフについて、詳しく説明します。

1. データと異常値のプロット:

このグラフは、元のデータと追加した異常値を含むデータポイントを表示しています。
以下の要素に注意してください：

  • 青い「x」マーカー:
    正常なデータポイントを表しています。これらは正規分布に従うランダムなデータです。
  • 赤い「o」マーカー:
    異常値を表しています。これらはわざとデータに追加された異常な値です。

このプロットは、異常検出のタスクを視覚的に理解するのに役立ちます。
異常値が通常のデータポイントからどれだけ外れているかが分かります。

2. 異常スコアのプロット:

このグラフは、Isolation Forestアルゴリズムによって計算された異常スコアを示しています。
以下の要素に注意してください：

  • 縦軸 (Y軸): 異常スコアを表しています。
    スコアが低いほどデータポイントは正常、高いほど異常と判断された可能性が高いです。
  • 横軸 (X軸): データポイントのインデックスを表しています。

グリーンの線で表されたプロットは、各データポイントの異常スコアを示しています。
異常スコアが閾値を超えるデータポイントは、異常として検出されたとみなされます。

これらのグラフは、異常検出のプロセスを理解し、異常検出アルゴリズムのパフォーマンスを評価する際に役立ちます。

異常スコアを基にして異常を検出するため、グラフから異常値がどれだけ明確に特定されているかを確認することが重要です。

飲食業の需要予測 Prophet

飲食業の需要予測

飲食業の需要予測の例題をProphetを使用して解いてみましょう。

まず、仮想のデータセットを用意しましょう。

例えば、過去4年間のあるレストランの月間売上データを持っているとします。

以下にデータの例を示します：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
日付,売上
2022-01-01,1000
2022-02-01,1200
2022-03-01,1100
2022-04-01,1300
2022-05-01,1400
2022-06-01,1500
2022-07-01,1300
2022-08-01,1200
2022-09-01,1100
2022-10-01,1000
2022-11-01,1200
2022-12-01,1300
2023-01-01,1400
2023-02-01,1500
2023-03-01,1600
2023-04-01,1400
2023-05-01,1300
2023-06-01,1200
2023-07-01,1100
2023-08-01,1000
2023-09-01,1200
2023-10-01,1300
2023-11-01,1400
2023-12-01,1500
2024-01-01,1600
2024-02-01,1700
2024-03-01,1500
2024-04-01,1400
2024-05-01,1300
2024-06-01,1200
2024-07-01,1100
2024-08-01,1000
2024-09-01,1200
2024-10-01,1300
2024-11-01,1400
2024-12-01,1500
2025-01-01,1600
2025-02-01,1700
2025-03-01,1800
2025-04-01,1600
2025-05-01,1500
2025-06-01,1400
2025-07-01,1300
2025-08-01,1200
2025-09-01,1100
2025-10-01,1000
2025-11-01,1200
2025-12-01,1300

それでは、PythonのProphetライブラリを使用してこのデータセットで需要予測を行ってみましょう。

まず、Prophetライブラリがインストールされていない場合は、次のコマンドを使用してインストールします：

1
!pip install prophet

インストールが完了したら、以下のコードを使って進めます。

以下は、提供されたデータセットを使用してProphetを使った需要予測を行う例です：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
import pandas as pd
from prophet import Prophet
import matplotlib.pyplot as plt

# データセットの定義
data = pd.read_csv('データセットのパス.csv')

# Prophet用にデータを準備する
data = data.rename(columns={'日付': 'ds', '売上': 'y'})

# Prophetモデルの作成
model = Prophet()

# データにモデルをフィットさせる
model.fit(data)

# 予測のための未来の日付を生成する
future_dates = model.make_future_dataframe(periods=12, freq='M')

# 予測を実行する
forecast = model.predict(future_dates)

# 予測結果をグラフ化する
model.plot(forecast, xlabel='Date', ylabel='Sales')
plt.title('Restaurant Demand Forecast')
plt.show()

この例では、まず必要なライブラリをインポートし、Prophetもインポートしています。

次にデータセットを読み込み、列名をProphetの予測モデルに合わせて変更しています。

その後、Prophetモデルを作成し、データにフィットさせ、予測のための未来の日付を生成します。

最後に、モデルを使用して未来の日付の売上を予測し、予測結果をグラフ化しています。

'データセットのパス.csv'の部分は、実際のデータセットファイルのパスに置き換えてください。

また、make_future_dataframeのperiodsパラメータを調整して、予測する未来の期間の数を指定することもできます。

ソースコード解説

以下、コードの各部分の詳細な説明です。

1. ライブラリのインポート:

1
2
3
import pandas as pd
from prophet import Prophet
import matplotlib.pyplot as plt
  • pandas: データ操作や分析のためのPythonライブラリ。
  • Prophet: Facebookが開発した時系列予測ライブラリ。このライブラリを使用して時系列モデルを構築し、予測を行う。
  • matplotlib.pyplot: グラフを描画するためのPythonライブラリ。

2. データセットの読み込み:

1
data = pd.read_csv('データセットのパス.csv')

CSV形式のデータセットを読み込み、dataという名前のDataFrame（表形式のデータ構造）に格納しています。

3. データの準備:

1
data = data.rename(columns={'日付': 'ds', '売上': 'y'})

データのカラム名をProphetの予測モデルに合わせるために変更しています。
'日付'カラムは時系列の日付を表し、'売上'カラムは予測する対象の数値データ（売上）を表しています。

4. Prophetモデルの作成:

1
model = Prophet()

Prophetのモデルをインスタンス化します。
このモデルは、時系列データからトレンド、季節性、ホリデーエフェクトなどを学習し、予測を行うために使用されます。

5. モデルのデータへのフィッティング（適合）:

1
model.fit(data)

dataに格納されたデータを使用して、Prophetモデルに適合させます。
モデルはここでデータのトレンドや季節性を学習します。

6. 未来の日付の生成:

1
future_dates = model.make_future_dataframe(periods=12, freq='M')

予測を行うための未来の日付を生成します。
periodsで指定した期間（ここでは12ヶ月）分の日付を生成し、freqで頻度を’M’（月ごと）として指定しています。

7. 予測の実行:

1
forecast = model.predict(future_dates)

未来の日付に対して予測を実行し、forecastというDataFrameに予測結果を格納します。

8. 予測結果のグラフ化:

1
2
3
model.plot(forecast, xlabel='Date', ylabel='Sales')
plt.title('Restaurant Demand Forecast')
plt.show()

予測結果をグラフ化して表示します。
x軸には日付を、y軸には売上を表示します。
また、グラフのタイトルも設定されています。

このコードは、Prophetを使用して時系列データの売上予測を行い、その結果を可視化する一連の手順を示しています。

結果解説

グラフの内容を詳しく説明します。

1. データポイント（黒い点）:

グラフ上に表示される黒い点は、実際の売上データポイントを表しています。

これらの点は、元のデータセットから取得され、予測モデルの作成に使用されました。

2. 予測値（青い線）:

グラフ上に表示される青い線は、Prophetモデルによって予測された売上の値を表しています。
この線は、過去のデータからのトレンドや季節性を考慮して、将来の売上を予測したものです。

3. 不確実性範囲（青い領域）:

グラフ上に表示される青い領域は、予測値の不確実性範囲を表しています。
Prophetは、予測の信頼性を示すために、予測値の上下に不確実性範囲を示します。
この範囲は、予測の信頼性が高いほど狭くなります。

このグラフを通じて、過去のデータからのトレンドや季節性を考慮した予測値が表示されます。

また、不確実性範囲も示されるため、予測の信頼性を評価することができます。

経済予測 Prophet

経済予測

Prophetを使用して経済予測を行う一般的な手順とコードの例を提供します。

以下はPythonコードの例です。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
import pandas as pd
from prophet import Prophet
import matplotlib.pyplot as plt

# 入力データの準備 (仮想的な経済データ)
data = pd.DataFrame({
    'ds': pd.date_range(start='2023-01-01', periods=24, freq='M'),
    'y': [100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 200, 210,
          220, 230, 240, 250, 260, 270, 280, 290, 300, 310, 320, 330]
})

# Prophetモデルの初期化
model = Prophet()

# モデルへのデータの適用
model.fit(data)

# 予測期間の準備
future = model.make_future_dataframe(periods=12, freq='M')

# 予測の実行
forecast = model.predict(future)

# 予測結果のグラフ化
fig = model.plot(forecast)
plt.title('Economic Forecast with Prophet')
plt.xlabel('Date')
plt.ylabel('Economic Indicator')
plt.show()

この例では、仮想的な経済データを使ってProphetモデルを構築し、将来の12ヶ月間の経済予測を行います。

データと予測結果はグラフで表示されます。

実際の経済データを使用する場合、適切なデータの前処理やハイパーパラメータの調整が必要です。

また、実際の経済データを使用する場合は、データを適切な形式で読み込んで上記の例に組み込む必要があります。

ソースコード解説

以下にコードの詳細な説明を提供します。

1. モジュールのインポート:

必要なライブラリやモジュールをインポートします。
pandasはデータ操作のために、ProphetはProphetモデルの構築と予測のために、matplotlib.pyplotはグラフの描画のために使用されます。

2. 入力データの準備:

仮想的な経済データを作成してpandasのDataFrameに格納します。
'ds'列には日付が、'y'列には対応する経済指標の値が入ります
ここでは24ヶ月分のデータを作成しています。

3. Prophetモデルの初期化:

Prophetクラスのインスタンスを作成して、Prophetモデルを初期化します。
このモデルは時系列データのトレンドや季節性をモデリングするために使用されます。

4. モデルへのデータの適用:

model.fit(data)を呼び出すことで、準備した経済データをモデルに適用して学習させます。
これにより、トレンドや季節性がモデルに学習されます。

5. 予測期間の準備:

model.make_future_dataframe()を使って、未来の予測期間の日付を含むDataFrameを作成します。
ここでは12ヶ月分の予測期間を設定しています。

6. 予測の実行:

model.predict(future)を呼び出すことで、予測期間における経済指標の予測値を計算します。
結果はforecastというDataFrameとして格納されます。

7. 予測結果のグラフ化:

model.plot(forecast)を呼び出すことで、予測結果をグラフとして表示します。
グラフには元データ、予測区間、トレンド、季節性などが含まれます。
matplotlibを使用してタイトルや軸ラベルを設定し、最終的にplt.show()でグラフを表示します。

このコードは、Prophetを使用して簡単に経済予測を行い、その結果を視覚化するための手順を示しています。

実際のデータを使用する場合は、データの前処理やハイパーパラメータの調整が必要です。

結果解説

経済予測のグラフの主要な要素の説明します。

1. 元データ (Observed Data):

グラフの黒い点線は、元のデータポイントを示しています。
これらのデータポイントは、予測モデルの構築に使用された実際の経済データです。
点線は元のデータポイントがある日付と対応する経済指標の値を示します。

2. 予測区間 (Prediction Interval):

グラフの青い領域は予測区間を示しています。
これは、予測モデルによって予測された値の範囲です。
通常、Prophetは不確実性を考慮して、予測値の上下に予測区間を描画します。
一般的に、将来のデータポイントがこの範囲内に入る確率は高いです。

3. トレンド (Trend):

グラフの青い線はトレンドを示しています。
Prophetは、時間の経過に伴う長期的な変動やトレンドをモデル化し、この青い線で表現しています。
トレンドは、データが上昇するか下降するかを示すものです。

Prophetの予測結果のグラフは、元のデータと予測されたトレンド、季節性、不確実性の情報を視覚的に表現するものです。

データの性質や用途に応じて、これらの要素がどのように現れるかが異なる場合があります。

製品売上予測 Prophet

製品売上予測

Prophetを使った製品売上予測の例を示します。

入力データ(日付と売上)と祝日情報をソースコード内に定義し、予測結果をグラフ化します。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
import pandas as pd
from prophet import Prophet
import matplotlib.pyplot as plt

# 入力データの定義（日付と売上）
data = pd.DataFrame({
    'ds': ['2023-01-01', '2023-01-02', '2023-01-03', '2023-01-04', '2023-01-05', '2023-01-06', '2023-01-07'],
    'y': [100, 150, 120, 130, 140, 160, 155]
})
data['ds'] = pd.to_datetime(data['ds'])

# 祝日情報の定義
holidays = pd.DataFrame({
    'holiday': 'holiday',
    'ds': pd.to_datetime(['2023-01-01', '2023-01-07']),
    'lower_window': 0,
    'upper_window': 0,
})

# モデルの作成
model = Prophet(holidays=holidays)
model.fit(data)

# 未来の日付を生成（30日分の予測）
future = model.make_future_dataframe(periods=30)

# 予測の実行
forecast = model.predict(future)

# 予測結果のグラフ化
fig = model.plot(forecast)
plt.title('Product Lifetime Prediction')
plt.xlabel('Date')
plt.ylabel('Sales')
plt.show()

このコードでは、入力データとして日付と売上を持つDataFrameを定義し、祝日情報も定義しています。

Prophetモデルを作成し、未来の日付を生成して予測を実行し、結果をグラフ化しています。

この例では7日間のデータを元に、さらに30日分の予測を行っています。

ご自身のデータに合わせて、日付と売上の値を適切に入力してください。

また、祝日情報もデータに合わせて調整してください。

ソースコード解説

それぞれのステップを詳しく説明します。

1. ライブラリのインポート:

  • pandas: データの操作と解析のためのライブラリ。
  • Prophet: Facebookが提供する時系列予測ライブラリ。
  • matplotlib.pyplot: データ可視化のためのライブラリ。

2. 入力データの定義:

  • data: 製品の売上データを持つDataFrame。
    ‘ds’列は日付を、’y’列は売上を表しています。
    日付は文字列型からpd.to_datetimeを使用して日付型に変換しています。

3. 祝日情報の定義:

  • holidays: 祝日の情報を持つDataFrame。
    ‘ds’列には祝日の日付を、’lower_window’と’upper_window’は祝日の前後の日数範囲を指定しています。
    ここでは祝日日付に影響を与えるために0を指定しています。

4. モデルの作成:

  • Prophetクラスのインスタンスを作成し、holidays情報を渡してモデルを初期化します。
  • fitメソッドでモデルをデータに適合させます。

5. 未来の日付の生成:

  • 予測結果を得るために、未来の日付を生成するためのDataFrame future を作成します。
    ここでは30日分の予測を行うため、periods引数に30を指定しています。

6. 予測の実行:

  • predictメソッドでモデルを用いて未来の売上を予測します。
    予測結果はforecastに格納されます。

7. 予測結果のグラフ化:

  • model.plot(forecast)で予測結果をグラフ化します。
  • plt.title、plt.xlabel、plt.ylabelでグラフのタイトルと軸ラベルを設定します。
  • plt.show()でグラフを表示します。

このコードは、元データの売上トレンド、未来の予測、および祝日の影響を一つのグラフで視覚的に表現するために使用されます。

データや祝日情報を適切に調整することで、自身のデータに対して予測を行うことができます。

結果解説

表示されるグラフは、Prophetによる製品寿命の売上予測の結果を表しています。

グラフには以下の要素が含まれています。

1. 黒い点（元データ）:

入力データにおける実際の売上を表す黒い点がプロットされています。
これらの点を通じて、元データの売上パターンが視覚的に確認できます。

2. 青いと線（予測データ）:

未来の予測日付に対する売上を表す青い点がプロットされています。
青い線は、Prophetによる予測のトレンドを示しています。

3. 青い領域（不確実性範囲）:

予測データの線の周りには、青い領域が描かれています。
これは予測の不確実性を示しており、売上がこれらの領域内に収まる確率が高いと言えます。

4. 祝日の影響:

グラフ上に縦線が描かれている部分は、祝日情報に基づいて設定された日付です。
これらの日付の周りには、祝日が予測にどのような影響を与えるかを示すエリアが表示されています。

このグラフを通じて、元データの売上トレンドや季節性、Prophetによる予測のトレンドと不確実性、祝日の影響などを視覚的に確認することができます。

予測結果がどの程度元データと一致するか、また祝日が売上にどのような影響を持つかを把握するのに役立ちます。

トラック運送の予測 Prophet

トラック運送の予測

トラック運送の予測を行うためのProphetのサンプルコードを提供します。

以下の例では、過去のトラック運送データと祝日の情報を使用して、未来のトラック運送量を予測します。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
import pandas as pd
from fbprophet import Prophet
import matplotlib.pyplot as plt

# 過去のトラック運送データ
data = {
    'ds': ['2023-01-01', '2023-01-02', '2023-01-03', '2023-01-04', '2023-01-05', '2023-01-06', '2023-01-07', '2023-01-08', '2023-01-09', '2023-01-10'],
    'y': [10, 20, 15, 25, 30, 35, 40, 45, 50, 55]
}

# Prophet用のデータフレームに変換
df = pd.DataFrame(data)
df['ds'] = pd.to_datetime(df['ds'])

# 祝日の情報を追加
holidays = pd.DataFrame({
  'holiday': 'holiday',
  'ds': pd.to_datetime(['2023-01-01', '2023-01-05']),
  'lower_window': 0,
  'upper_window': 1,
})

# モデルの作成と学習
model = Prophet(holidays=holidays)
model.fit(df)

# 予測期間の設定
future = model.make_future_dataframe(periods=365)

# 予測の実行
forecast = model.predict(future)

# 予測結果のグラフ化
fig, ax = plt.subplots(figsize=(10, 6))
model.plot(forecast, ax=ax)
model.plot_components(forecast, ax=ax)
plt.title('Truck Delivery Forecast')
plt.xlabel('Date')
plt.ylabel('Delivery Volume')
plt.show()

上記のコードでは、過去のトラック運送データをdataとして定義し、Prophet用のデータフレームに変換しています。

また、holidaysというデータフレームを作成し、祝日の情報を追加しています。

モデルの作成時にholidaysを指定し、make_future_dataframeで予測期間を設定しています。

予測の実行後、plotメソッドを使用して予測結果をグラフ化しています。

さらに、plot_componentsメソッドを使用してトレンドや季節性の成分をグラフ化しています。

このサンプルコードを実行すると、トラック運送の予測結果がグラフとして表示されます。





グラフを通じて、予測値と実測値の関係やトレンドの変化、祝日の影響などを視覚的に確認することができます。

ご自身のデータに合わせて適切にカスタマイズし、トラック運送の予測結果をグラフ化してみてください。

ソースコード解説

コードの各部分について説明します。

1. ライブラリのインポート:

  • import pandas as pd: Pandasライブラリをインポートして、データ処理に使用します。
  • from prophet import Prophet: Prophetクラスをインポートして、予測モデルを構築・訓練するために使用します。
  • import matplotlib.pyplot as plt: Matplotlibライブラリをインポートして、グラフの描画に使用します。

2. 日次のデータ生成:

  • pd.date_range: 2020年1月1日から2023年1月1日までの日次の日付データを生成します。
  • pd.DataFrame: 生成した日付データとそのインデックスに対応する連番を含むDataFrameを作成します。
    列名は ‘ds’ と ‘y’ です。’ds’ は日付を、’y’ はその日の連番を表します。

3. モデルのインスタンス化:

  • model = Prophet(): Prophetモデルのインスタンスを作成します。

4. 季節性の追加:

  • model.add_seasonality: モデルに週次と年次の季節性を追加します。
    週次季節性の周期は7日で、フーリエオーダーが3です。
    年次季節性の周期は365.25日（うるう年を考慮）で、フーリエオーダーが10です。

5. 祝日の指定:

  • pd.DataFrame: ‘special_day’ という名前の祝日を指定します。
    祝日の日付は、2020年から2022年までの12月25日です。
    lower_window と upper_window は祝日の前後に適用するウィンドウの日数を指定します。
  • model.add_country_holidays: モデルに国に基づく祝日（この場合はアメリカ）を追加します。

6. データへの適合:

  • model.fit(data): モデルを提供されたデータに適合させます。

7. 予測期間の作成:

  • model.make_future_dataframe: 予測期間のDataFrameを生成します。
    ここでは365日（1年間）の期間を指定しています。

8. 予測の実行:

  • model.predict(future): 予測期間内の予測を実行し、予測結果を得ます。

9. 予測結果の可視化:

  • 予測結果を含むグラフを描画し、plt.title や plt.xlabel でタイトルと横軸ラベルを設定します。

10. 季節性の可視化:

  • model.plot_components: 季節性を構成要素ごとに可視化するグラフを作成します。

以上の手順によって、Prophetモデルを使用して時系列データの予測と可視化が行われます。

結果解説

表示される5つのグラフについて説明します。

1. 予測結果の折れ線グラフ（Forecast Plot）：

このグラフは、実測値と予測値を折れ線で表示します。
実測値は黒い点で表され、予測値は青い折れ線で表されます。
このグラフを通じて、予測値が実測値にどれくらい近いかを確認することができます。

2. トレンドのグラフ（Trend Plot）：

このグラフは、時間の経過に伴うデータのトレンドを示します。
予測値のトレンドは青い折れ線で表され、過去のデータから未来の予測までの傾向を示します。
このグラフを通じて、データの増減や変化の傾向を把握することができます。

3. 季節性のグラフ（Seasonality Plot）：

このグラフは、データの季節性パターンを示します。
季節性は、周期的な変動やパターンを持つデータの特徴です。
Prophetは、季節性を捉えるためにフーリエ変換を使用しています。
このグラフを通じて、データの季節性の変化や周期性を確認することができます。

4. 週次季節性のグラフ（Weekly Seasonality Plot）：

このグラフは、週ごとの季節性パターンを示します。
週次季節性は、週単位でのデータの変動やパターンを表します。
Prophetは、週次季節性を捉えるためにフーリエ変換を使用しています。
このグラフを通じて、週ごとのデータの変動や季節性のパターンを確認することができます。

5. 月次季節性のグラフ（Monthly Seasonality Plot）：

このグラフは、月ごとの季節性パターンを示します。
月次季節性は、月単位でのデータの変動やパターンを表します。
Prophetは、月次季節性を捉えるためにフーリエ変換を使用しています。
このグラフを通じて、月ごとのデータの変動や季節性のパターンを確認することができます。

これらのグラフを組み合わせることで、データのトレンドや季節性、予測値の信頼性などを詳しく分析することができます。

また、グラフの表示によって、データの特徴や傾向を視覚的に把握することができます。

オンライン教育の利用予測 Prophet

オンライン教育の利用予測

オンライン教育の利用予測の例題を解決するために、Facebookが開発したProphetというツールを使用します。

Prophetは、時系列データの予測に特化したオープンソースの予測モデルです。

以下の例では、”Date”と”Usage”の2つの列を持つデータセットを想定しています。

“Date”列は各データポイントのタイムスタンプを表し、”Usage”列はその日付のオンライン教育利用者数を表します。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
import pandas as pd
from prophet import Prophet

# 入力データを定義する
data = pd.DataFrame({
    'ds': ['2023-01-01', '2023-01-08', '2023-01-15', '2023-01-22', '2023-01-29', '2023-02-05', '2023-02-12', '2023-02-19', '2023-02-26', '2023-03-05'],
    'y': [100, 120, 110, 130, 150, 140, 160, 180, 170, 190]
})

# 'ds'列を日付型に変換する
data['ds'] = pd.to_datetime(data['ds'])

# Prophetモデルを初期化してフィットさせる
model = Prophet(weekly_seasonality=True)
model.fit(data)

# 予測用の未来のデータフレームを作成する
future = model.make_future_dataframe(periods=10, freq='W')

# 予測を行う
forecast = model.predict(future)

# 予測結果をプロットする
fig = model.plot(forecast)

この例では、2023年1月1日から2023年3月5日までの10個のデータポイントを持つ入力データを定義しています。

weekly_seasonality=Trueパラメータは、モデルがデータの週次周期性を考慮するようにします。

make_future_dataframe関数は、10週間の追加予測用データフレームを作成するために使用されます。

freq='W'パラメータは、頻度を週次に設定します。

最後に、plot関数を使用して予測結果をグラフで可視化します。

これは単なる基本的な例であり、具体的な要件に応じて入力データやパラメータをカスタマイズすることができます。

ソースコード解説

提供されたコードは、Prophetという時系列予測モデルを用いて、一連の時系列データの予測を行う基本的な手順を示しています。以下に各部分の詳細説明を示します。

1. 必要なライブラリのインポート:

まず、必要なライブラリをインポートしています。
pandas ライブラリはデータの操作や整形、Prophet ライブラリは時系列予測モデルを構築するために使用されます。

2. 入力データの定義:

data という名前のデータフレームを作成しています。
このデータフレームには、’ds’ 列に日付情報、’y’ 列に対応する時系列データ（ここでは例として利用回数）を含むデータが格納されています。

3. 日付データの整形:

‘ds’ 列は日付を文字列として格納しています。
そのため、pd.to_datetime を使用してこの列を日付型に変換しています。
これにより、日付データを効果的に処理できるようになります。

4. Prophetモデルの初期化とフィッティング:

Prophet クラスのインスタンスを作成し、weekly_seasonality=True を設定して週次の周期性を考慮するようにしています。
その後、このモデルに data データフレームを適合させています。

5. 予測用の未来データの作成:

予測のために、元のデータ期間に加えて未来の期間を含むデータフレーム future を作成しています。periods=10 と freq='W' を指定して、未来の10週間分の日付データを生成しています。

6. 予測の実行:

model.predict(future) を呼び出すことで、未来の期間に対する予測を行います。
予測結果は forecast というデータフレームとして返されます。

7. 予測結果のプロット:

model.plot(forecast) を呼び出すことで、元のデータと予測結果を含むグラフが表示されます。

このグラフを通じて、予測の範囲、予測された時系列データの変化、週次の季節性などが視覚的に確認できます。

このコードは、Prophetモデルを使って時系列データの予測を行う基本的な手法を示しており、実際のデータに適用する際には、データの特性に合わせて適切な調整が必要です。

結果解説

グラフは、以下の要素を含んでいます。

1. 黒い点:

元のデータポイントを表しており、各データポイントの x 軸には日付（’ds’ 列の値）、y 軸には対応する時系列データの値（’y’ 列の値）が配置されています。

2. 青い線と青い範囲:

予測の結果として、黒い線が描かれ、これは予測された時系列データの値を示しています。
青い範囲は、予測の不確実性を示すもので、予測された値の上下に広がっています。

3. 縦軸 (y 軸):

利用回数などの時系列データの値を表します。

4. 横軸 (x 軸):

日付を示し、時間の経過に伴うデータの変化を示しています。
青い点と黒い線が連続している範囲が、過去から未来への時系列の推移を表しています。

5. 週次の周期性:

グラフ上に縦線が表示され、週ごとの周期性を示しています。
この週次の周期性は、Prophetモデルによって考慮されています。

このグラフを通じて、元のデータと予測されたデータの変化、および週ごとの季節性を視覚的に理解することができます。

飲食業の需要予測 Prophet

飲食業の需要予測

飲食業の需要予測の例を示し、Prophetを使用して予測し、結果をグラフ化する手順を説明します。

まず、必要なライブラリをインポートします。

1
2
3
4
import pandas as pd
import numpy as np
from prophet import Prophet
import matplotlib.pyplot as plt

次に、仮想的な飲食業の需要データを生成します。

1
2
3
4
5
6
7
8
9
10
np.random.seed(10)  # 乱数生成のシードを固定して再現性を保つ

# 週次周期の日付データを生成
dates = pd.date_range(start='2023-01-01', periods=104, freq='W-SAT')

# ランダムな需要データを生成
demand = np.random.randint(100, 1000, size=len(dates))

# データフレームに格納
data = pd.DataFrame({'ds': dates, 'y': demand})

Prophetを使用して予測を行い、結果をグラフ化します。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
# Prophetモデルの初期化
model = Prophet(weekly_seasonality=True)

# モデルにデータをフィット
model.fit(data)

# 予測用データフレームの作成
future = model.make_future_dataframe(periods=12, freq='W-SAT')  # 予測を12週先まで行う

# 予測を実行
forecast = model.predict(future)

# グラフ化
fig = model.plot(forecast)
plt.xlabel('Date')
plt.ylabel('Demand')
plt.title('Restaurant Demand Forecast with Weekly Seasonality')
plt.show()

上記のコードでは、仮想的な飲食業の需要データを生成し、Prophetモデルを使用して予測を行いました。

生成されたグラフは、日付に対する需要の予測を示しています。

実際のデータやニーズに合わせてコードを調整してください。

ソースコード解説

ソースコードの各部分の詳細を説明します。

1. ライブラリのインポート:

  • pandasはデータ操作のためのライブラリです。
  • numpyは数値計算のためのライブラリです。
  • ProphetはFacebookが開発した時系列予測ライブラリです。
  • matplotlib.pyplotはデータの可視化のためのライブラリです。

2. 乱数生成のシードの固定:

  • np.random.seed(10)は、乱数生成のシードを固定して再現性を保つためのコードです。

3. 週次周期の日付データの生成:

  • pd.date_rangeを使用して、’2023-01-01’から104週間分の週次周期の日付データを生成します。

4. ランダムな需要データの生成:

  • np.random.randintを使用して、100から1000までのランダムな需要データを生成します。

5. データフレームへの格納:

  • pd.DataFrameを使用して、日付データと需要データを組み合わせてデータフレームを作成します。

6. Prophetモデルの初期化:

  • Prophet(weekly_seasonality=True)を使用して、Prophetモデルを初期化します。
    weekly_seasonality=Trueは週次周期性を考慮するためのパラメータです。

7. モデルへのデータのフィット:

  • model.fit(data)を使用して、データをモデルにフィットさせます。

8. 予測用データフレームの作成:

  • model.make_future_dataframeを使用して、12週先までの予測用のデータフレームを作成します。

9. 予測の実行:

  • model.predictを使用して、予測を実行します。

10. グラフ化:

  • model.plotを使用して、予測結果をグラフ化します。
  • plt.xlabelとplt.ylabelを使用して、x軸とy軸のラベルを設定します。
  • plt.titleを使用して、グラフのタイトルを設定します。
  • plt.showを使用して、グラフを表示します。

このソースコードは、ランダムな需要データを使用してProphetモデルを構築し、週次周期性を考慮した需要の予測を行います。

予測結果はグラフとして表示されます。

ただし、具体的なデータやパラメータによって結果は異なるため、実際のデータに合わせて調整してください。

結果解説

グラフの特徴や内容について説明いたします。

Prophetによる需要予測のグラフは、横軸が日付（週次で表示されている場合は週ごとの日付）、縦軸が需要（または売上など）を表しています。

以下に、典型的なProphetの予測グラフで見られる要素について説明します。

1. 黒い点線:

実測データの点線が表示されます。
これは、過去の実際の需要データを示しています。

2. 青い線と範囲:

黒い線は、Prophetモデルによって予測された需要の中央値（予測値）を表しています。
また、周囲に淡い色で描かれた領域は、不確実性を示しており、予測の信頼区間を表しています。

3. 週次周期性の影響:

グラフに週次周期性の影響が現れます。
通常、週末に需要が上昇する業界では、週末にピークが現れる傾向があるでしょう。

4. 未来の予測:

青い線が実測データの最後から伸びている部分は、未来の予測を示しています。
これにより、将来の需要の予測値と信頼区間を確認することができます。

5. トレンド:

グラフ内のデータポイントや線の動きに沿ったトレンドが観察できるかもしれません。

需要が増加している、減少している、あるいは安定しているかを確認できます。

グラフを通じて、過去の需要データからProphetモデルが予測した将来の需要の変動パターンと、それに伴う不確実性を視覚的に理解することができます。

あなたのデータに合わせてこのグラフを解釈し、適切なビジネス判断を行うことが重要です。

自然災害（地震）の発生予測 Prophet

自然災害（地震）の発生予測

自然災害（地震）の発生予測の例題を示します。

こちらはダミーデータを使用しています。

また、年次周期性を考慮してProphetを適用し、予測結果をグラフで表示します。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
import pandas as pd
from prophet import Prophet
import matplotlib.pyplot as plt

# ダミーデータを作成
data = {'ds': pd.date_range(start='2023-01-01', periods=365, freq='D'),
        'y': [0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0,
              1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0,
              0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0,
              0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0,
              0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0,
              0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0,
              0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0,
              0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
              0, 0, 0, 0, 0]}  # 365日間の地震発生（1）・非発生（0）
print(len(data['y']))
df = pd.DataFrame(data)

# Prophetモデルの作成
model = Prophet(
    changepoint_prior_scale=0.5,  # 変化点の影響度
    yearly_seasonality=True,      # 年次周期性を考慮
    weekly_seasonality=False      # 週次周期性は考慮しない
)
model.add_seasonality(name='monthly', period=30.5, fourier_order=5)  # 月次周期性を考慮

model.fit(df)

# 予測用データフレームを作成
future = model.make_future_dataframe(periods=30)  # 30日間の予測
forecast = model.predict(future)

# グラフで表示
fig = model.plot(forecast)
plt.title('Earthquake Occurrence Prediction')
plt.xlabel('Date')
plt.ylabel('Occurrence (1) / Non-occurrence (0)')
plt.ylim([-0.1, 1.1])
plt.show()

この例では、365日間の地震発生データ（1：地震発生、0：地震非発生）を使用してProphetモデルを作成し、その後30日間の地震発生を予測しています。

年次周期性を考慮し、月次の周期性も加えています。

予測結果のグラフは、年次の周期性や月次の変動を反映しています。

地震の発生・非発生のパターンを視覚的に理解することができます。

ただし、この例では実際の地震発生データではなく、ダミーデータを使用しています。

実際の地震発生予測には多くの要因が影響し、複雑なモデリングが必要です。

ソースコード解説

このソースコードは、地震発生予測の例題を示しています。

プログラムはPythonで記述されており、Prophetライブラリを使用して地震の発生予測を行い、結果をグラフで表示しています。

以下にコードの各部分を説明します。

1. import文:

必要なライブラリをインポートしています。
pandasはデータの操作、Prophetは時系列予測モデルの作成、matplotlib.pyplotはグラフの描画に使用されます。

2. ダミーデータの作成:

data辞書を用いて、日付と地震の発生（1）・非発生（0）のデータを作成しています。

3. dfの作成:

data辞書を使用してpandasのデータフレーム df を作成しています。
これには日付（’ds’）と地震の発生状況（’y’）が含まれています。

4. Prophetモデルの作成:

Prophetモデルを作成します。
changepoint_prior_scaleは変化点の影響度、yearly_seasonalityは年次周期性、weekly_seasonalityは週次周期性を指定しています。

5. 月次周期性の追加:

add_seasonalityメソッドを使用して、月次周期性を追加しています。
periodは周期の長さ（30.5日）、fourier_orderは周期性の調整パラメータです。

6. モデルの学習:

fitメソッドを使用して、モデルをデータに適合させています。

7. 予測用データフレームの作成:

make_future_dataframeメソッドを使用して、将来の予測データを含むデータフレーム future を作成しています。

8. 予測の実行:

predictメソッドを使用して、予測を行い、forecastに結果を格納しています。

9. グラフの描画:

model.plot(forecast)で予測結果のグラフを描画しています。
plt.title、plt.xlabel、plt.ylabelなどを使用して、グラフのタイトルと軸ラベルを設定しています。
また、plt.ylimで縦軸の範囲を設定しています。

10. plt.show():

グラフを表示します。

このプログラムは、Prophetモデルを使用して地震発生予測を行い、その結果を日付ごとの地震発生（1）・非発生（0）の状況としてグラフで可視化しています。

結果解説

地震発生予測のグラフについて説明します。

このグラフは、Prophetモデルを使用して予測された地震の発生（1）または非発生（0）を表しています。

横軸は日付を示し、縦軸は地震の発生を表しています。

以下にグラフ内の主な要素を説明します。

1. 黒い点線:

実際のデータポイントを示します。
黒い点線の位置が、過去のデータに基づいてモデルが予測した地震発生の状態を表します。

2. 青い線（予測）:

予測された地震発生の状態を示す黒い線です。
この線は、過去のデータからモデルが学習したパターンと、その後の予測期間に対する予測を基にして描かれています。

3. 青い領域:

青い領域は、予測の信頼区間を表しています。
予測は確実性を持って行えるわけではないため、信頼区間内での予測が行われていることを示します。

4. 横軸（日付）:

日付が横軸に表示されており、過去のデータと将来の予測期間が表示されています。

5. 縦軸（地震発生）:

縦軸は地震の発生（1）または非発生（0）を表しています。
縦軸の値が0または1の間で変動します。

このグラフを通じて、過去の地震発生データに基づいてモデルが学習し、将来の予測期間における地震発生の傾向や確からしさを示すことができます。

信頼区間内の予測結果がどれだけ合致しているかを確認し、地震発生のパターンを把握する助けとなります。

人材需要の予測 Prophet

人材需要の予測

人材需要の予測を行いましょう。

以下の例をご覧ください。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
import pandas as pd
from prophet import Prophet
import numpy as np

# ランダムな入力データの生成
np.random.seed(0)
dates = pd.date_range(start='2021-01-01', periods=24, freq='M')
data = pd.DataFrame({
    'ds': dates,
    'y': np.random.randint(low=50, high=200, size=24)
})

# Prophetモデルの作成と学習
model = Prophet()
model.fit(data)

# 予測期間の設定
future = model.make_future_dataframe(periods=12, freq='M')

# 予測の実行
forecast = model.predict(future)

# グラフ化
fig = model.plot(forecast)
fig.show()

この例では、24ヶ月分のランダムな人材需要のデータを生成し、次の12ヶ月間の予測を行います。

ds列には日付、y列にはランダムな人材需要の指標が入っています。

Prophetモデルを使ってデータを学習し、make_future_dataframe関数で予測期間を設定します。

そして、predict関数を使って予測を実行します。

最後に、plot関数を使って予測結果をグラフ化し、show関数でグラフを表示します。

ソースコード解説

このソースコードは、Prophetを使用してランダムな人材需要の予測を行うための手順を示しています。

以下、各行の詳細な説明です。

1. import pandas as pd:

pandasライブラリをインポートし、データの操作と解析のための機能を提供します。

2. from prophet import Prophet:

ProphetライブラリからProphetクラスをインポートします。
Prophetは、時系列データの予測モデリングに特化したライブラリです。

3. import numpy as np:

numpyライブラリをインポートし、数値計算やランダムなデータの生成などの機能を提供します。

4. np.random.seed(0):

乱数のシードを設定します。これにより、再現性のあるランダムなデータが生成されます。

5. dates = pd.date_range(start='2021-01-01', periods=24, freq='M'):

pandasのdate_range関数を使って、24ヶ月分の日付データを生成します。
開始日は’2021-01-01’で、1ヶ月ごとにデータが生成されます。

6. data = pd.DataFrame({'ds': dates, 'y': np.random.randint(low=50, high=200, size=24)}):

pandasのDataFrameを使って、日付データとランダムな人材需要のデータを組み合わせてデータフレームを作成します。
ds列には日付データが、y列には50から200の範囲でランダムに生成された人材需要のデータが格納されます。

7. model = Prophet():

Prophetクラスのインスタンスを作成し、予測モデルを初期化します。

8. model.fit(data):

Prophetモデルにデータを適合させて学習を行います
。ここで、過去のデータを元にトレンドや季節性を学習します。

9. future = model.make_future_dataframe(periods=12, freq='M'):

予測期間を設定するために、将来の日付データを生成します。
periods引数には予測する期間の長さを指定し、freq引数には期間の間隔を指定します。

10. forecast = model.predict(future):

予測を実行し、将来の人材需要を予測します。
predict関数は、学習済みのモデルと将来の日付データを入力として、予測結果を返します。

11. fig = model.plot(forecast):

予測結果をグラフ化します。
plot関数は、予測結果を可視化するためのグラフを作成します。

12. fig.show(): グラフを表示します。

以上がソースコードの詳細な説明です。

このコードを実行することで、ランダムな人材需要の予測結果をグラフで表示することができます。

結果解説

結果のグラフについて説明します。

Prophetによって生成されたグラフには、以下の要素が含まれています。

1. 元のデータポイント（Observed Data Points）:

グラフ上の黒い点が元のデータポイントを表しています。
これは、予測モデルの学習に使用された実際の人材需要のデータです。

2. 予測値（Trend Forecast）:

グラフ上の黒い線が予測値を表しています。
この線は、Prophetモデルによって予測された人材需要のトレンドを示しています。
トレンドは、データの中で見られる長期的な変化や傾向を表しています。

3. 予測の不確実性（Uncertainty Intervals）:

グラフ上の青い領域が予測の不確実性を示しています。
これは、予測値の周りに描かれた領域で、将来のデータポイントがどの程度のばらつきを持つかを示しています。
一般的に、不確実性が高いほど予測の信頼性が低くなります。

これらの要素を組み合わせることで、Prophetによる人材需要の予測結果が視覚化されます。
グラフを通じて、トレンドや季節性のパターン、予測の不確実性などを把握することができます。

以上が結果のグラフの詳細な説明です。

予測結果を分析する際には、これらの要素を考慮に入れてください