次の形式のデータセット(このリンクにあります https://drive.google.com/open?id=0B2Iv8dfU4fTUY2ltNGVkMG05V )があります。
_ time X Y
0.000543 0 10
0.000575 0 10
0.041324 1 10
0.041331 2 10
0.041336 3 10
0.04134 4 10
...
9.987735 55 239
9.987739 56 239
9.987744 57 239
9.987749 58 239
9.987938 59 239
_
データセットの3番目の列(Y)は私の真の値です。それが私が予測(推定)したかったものです。 Y
の予測を実行したい(つまり、Y
の前の100個のローリング値に従ってX
の現在の値を予測します。このために、次のpython
スクリプトは_random forest regression model
_を使用して動作します。
_#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
@author: deshag
"""
import pandas as pd
import numpy as np
from io import StringIO
from sklearn.ensemble import RandomForestRegressor
from sklearn.metrics import mean_squared_error
from math import sqrt
df = pd.read_csv('estimated_pred.csv')
for i in range(1,100):
df['X_t'+str(i)] = df['X'].shift(i)
print(df)
df.dropna(inplace=True)
X=pd.DataFrame({ 'X_%d'%i : df['X'].shift(i) for i in range(100)}).apply(np.nan_to_num, axis=0).values
y = df['Y'].values
reg = RandomForestRegressor(criterion='mse')
reg.fit(X,y)
modelPred = reg.predict(X)
print(modelPred)
print("Number of predictions:",len(modelPred))
meanSquaredError=mean_squared_error(y, modelPred)
print("MSE:", meanSquaredError)
rootMeanSquaredError = sqrt(meanSquaredError)
print("RMSE:", rootMeanSquaredError)
_
最後に、二乗平均平方根誤差(RMSE)を測定し、RMSE
の_19.57
_を取得しました。私がドキュメントから読んだことから、二乗誤差は応答と同じ単位を持つと言います。 RMSE
の値をパーセントで表示する方法はありますか?たとえば、予測のこのパーセントは正しいが、これはかなり間違っていると言います。
sklearn
の最近のバージョンには、mean absolute percentage error (MAPE)
を計算するための_check_array
_関数がありますが、以前のバージョンと同じように動作しないようです。以下。
_import numpy as np
from sklearn.utils import check_array
def calculate_mape(y_true, y_pred):
y_true, y_pred = check_array(y_true, y_pred)
return np.mean(np.abs((y_true - y_pred) / y_true)) * 100
calculate_mape(y, modelPred)
_
これはエラーを返します:ValueError: not enough values to unpack (expected 2, got 1)
。また、以前のバージョンとは異なり、最近のバージョンの_check_array
_関数は 単一の値 のみを返すようです。
RMSE
をパーセントで表示したり、MAPE
をsklearn
に対してPython
を使用して計算したりする方法はありますか?
calculate_mape
で削除されたcheck_arrays
関数が予期されているため、sklearn 0.16
の実装は機能していません。 check_array
はあなたが望むものではありません。
This StackOverflowの回答は、有効な実装を提供します。