<https://archive.ics.uci.edu/dataset/235/individual+household+electric+power+consumption> 上記のサイトから取得できるデータセットを用いて電力データの「時系列予測」を行っていて、 'Global_active_power','Sub_metering_1','Sub_metering_2','Sub_metering_3'の4つのデータから 'Global_active_power'の予測を行いたくPython環境でKerasのライブラリのLSTMを使い「時系列予測」を行っているのですが、4つのデータから1つのデータの予測のやり方がわからないです。 現在のコードが以下のようになっています。4つのデータから4つをそのまま予測しているようです。 どのようにしたら4つのデータから1つのデータを予測できるでしょうか。ご教授してもらえると幸いです。 追記 参考にしているサイトが以下のサイトです <https://www.kaggle.com/code/yassinesfaihi/lstm-time-series-household-power-consumption> 追追記 時刻t-1からt-5の間の'Global_active_power','Sub_metering_1','Sub_metering_2','Sub_metering_3'の計4つのデータから 時刻tの'Global_active_power'を予測したいと思っています。 ``` # Import necessary libraries and packages import numpy as np import matplotlib.pyplot as plt import pandas as pd # Set floating point precision option for pandas pd.set_option('display.float_format', lambda x: '%.4f' % x) # Import seaborn library and set context and style import seaborn as sns sns.set_context("paper", font_scale=1.3) sns.set_style('white') # Import warnings and set filter to ignore warnings import warnings warnings.filterwarnings('ignore') # Import time library from time import time # Import matplotlib ticker and scipy stats import matplotlib.ticker as tkr from scipy import stats # Import statistical tools for time series analysis from statsmodels.tsa.stattools import adfuller # Import preprocessing from sklearn from sklearn import preprocessing # Import partial autocorrelation function from statsmodels from statsmodels.tsa.stattools import pacf # Enable inline plotting in Jupyter Notebook %matplotlib inline # Import math library import math # Import necessary functions from keras import keras from keras.models import Sequential from keras.layers import Dense from keras.layers import LSTM from keras.layers import Dropout from keras.layers import * # Import MinMaxScaler from sklearn from sklearn.preprocessing import MinMaxScaler # Import mean squared error and mean absolute error from sklearn from sklearn.metrics import mean_squared_error from sklearn.metrics import mean_absolute_error # Import early stopping from keras callbacks from keras.callbacks import EarlyStopping #結果の読み出し(前回の続きから出力) data = pd.read_csv('/content/drive/MyDrive/Research_folder/Data_set/Data_set_from_UCI/household_power_consumption_Edited.csv') data.head(5) dataset = data.loc[:,['Global_active_power','Sub_metering_1','Sub_metering_2','Sub_metering_3']] dataset = dataset.values.astype('float32') #Create an instance of the MinMaxScaler class to scale the values between 0 and 1 scaler = MinMaxScaler(feature_range=(0, 1)) #Fit the MinMaxScaler to the transformed data and transform the values dataset = scaler.fit_transform(dataset) train_size = int(len(dataset) * 0.80) test_size = len(dataset) - train_size train, test = dataset[0:train_size,:], dataset[train_size:len(dataset),:] # convert an array of values into a dataset matrix def create_dataset(dataset, look_back=1): X, Y = [], [] for i in range(len(dataset)-look_back-1): a = dataset[i:(i+look_back), 0] X.append(a) Y.append(dataset[i + look_back, 0]) return np.array(X), np.array(Y) # reshape into X=t and Y=t+1 look_back = 30 X_train, Y_train = create_dataset(train, look_back) X_test, Y_test = create_dataset(test, look_back) # reshape input to be [samples, time steps, features] X_train = np.reshape(X_train, (X_train.shape[0], 1, X_train.shape[1])) X_test = np.reshape(X_test, (X_test.shape[0], 1, X_test.shape[1])) # Defining the LSTM model model = Sequential() # Adding the first layer with 100 LSTM units and input shape of the data model.add(LSTM(100, input_shape=(X_train.shape[1], X_train.shape[2]))) # Adding a dropout layer to avoid overfitting model.add(Dropout(0.2)) # Adding a dense layer with 1 unit to make predictions model.add(Dense(1)) # Compiling the model with mean squared error as the loss function and using Adam optimizer model.compile(loss='mean_squared_error', optimizer='adam') # Fitting the model on training data and using early stopping to avoid overfitting history = model.fit(X_train, Y_train, epochs=20, batch_size=1240, validation_data=(X_test, Y_test), callbacks=[EarlyStopping(monitor='val_loss', patience=4)], verbose=1, shuffle=False) # Displaying a summary of the model model.summary() # make predictions train_predict = model.predict(X_train) test_predict = model.predict(X_test) pad_col = np.zeros(dataset.shape[1]-1) def pad_array(val): return np.array([np.insert(pad_col, 0, x) for x in val]) train_predict = scaler.inverse_transform(pad_array(train_predict)) Y_train = scaler.inverse_transform(pad_array(Y_train)) test_predict = scaler.inverse_transform(pad_array(test_predict)) Y_test = scaler.inverse_transform(pad_array(Y_test)) aa=[x for x in range(200)] # Creating a figure object with desired figure size plt.figure(figsize=(20,6)) plt.tick_params(labelsize=20) # Plotting the actual values in blue with a dot marker plt.plot(aa, Y_test[0][:200], marker='.', label="actual", color='purple') # Plotting the predicted values in green with a solid line plt.plot(aa, test_predict[:,0][:200], '-', label="prediction", color='red') # Removing the top spines sns.despine(top=True) # Adjusting the subplot location plt.subplots_adjust(left=0.07) # Labeling the y-axis plt.ylabel('Global_active_power', size=22) # Labeling the x-axis plt.xlabel('Time step', size=22) # Adding a legend with font size of 15 plt.legend(fontsize=22) # Display the plot plt.show() ```