import pandas as pd
import numpy as np
from sklearn.preprocessing import StandardScaler
import tensorflow as tf
from tensorflow.keras.models import Sequential
from tensorflow.keras.layers import SimpleRNN, Dense, Input
# Load the Air Quality dataset from the URL
url = "https://codesignal-staging-assets.s3.amazonaws.com/uploads/1742293523899/AirQualityUCI.csv"
df = pd.read_csv(url, sep=';', decimal=',')
# Replace -200 with NaN
df.replace(-200, np.nan, inplace=True)
# Combine 'Date' and 'Time' into a single 'DateTime' column
df['DateTime'] = pd.to_datetime(df['Date'] + ' ' + df['Time'], format="%d/%m/%Y %H.%M.%S")
# Drop the original 'Date' and 'Time' columns
df.drop(columns=['Date', 'Time'], inplace=True)
# Drop unnecessary columns
df.drop(columns=['Unnamed: 15', 'Unnamed: 16'], inplace=True)
# Drop rows where essential features are missing
df.dropna(subset=['CO(GT)', 'NO2(GT)', 'T', 'RH'], inplace=True)
# Fill missing values using forward-fill & backward-fill
df.ffill(inplace=True)
df.bfill(inplace=True)
# Set DateTime as index
df.set_index('DateTime', inplace=True)
# Select relevant features for RNN
features = ['CO(GT)', 'NO2(GT)', 'PT08.S5(O3)', 'RH', 'T']
df_selected = df[features]
# Normalize the data
scaler = StandardScaler()
df_scaled = scaler.fit_transform(df_selected)
# Function to create sequences for multi-input RNN
def create_multivariate_sequences(data, seq_length=10):
X, y = [], []
for i in range(len(data) - seq_length):
X.append(data[i:i+seq_length]) # Multi-feature input
y.append(data[i+seq_length, -1]) # Predicting Temperature (T)
return np.array(X), np.array(y)
# Define sequence length
sequence_length = 10
# Create sequences
X, y = create_multivariate_sequences(df_scaled, sequence_length)
# Reshape input for RNN
X = X.reshape((X.shape[0], sequence_length, len(features))) # (samples, timesteps, features)
# Split into training and test sets
train_size = int(0.8 * len(X))
X_train, X_test = X[:train_size], X[train_size:]
y_train, y_test = y[:train_size], y[train_size:]
# Define a modified RNN model for multivariate time series
model = Sequential([
Input(shape=(sequence_length, len(features))),
SimpleRNN(30, activation='relu'),
Dense(10, activation='relu'),
Dense(1)
])
# Compile model
model.compile(optimizer='adam', loss='mse')
# Train model and obtain the history object
history = model.fit(X_train, y_train, epochs=5, batch_size=16, verbose=1)
# Print model summary
model.summary()