Spaces:
Running
Running
import streamlit as st | |
import pandas as pd | |
import time | |
from datetime import datetime | |
import numpy as np | |
import pmdarima as pm | |
import matplotlib.pyplot as plt | |
from pmdarima import auto_arima | |
import plotly.graph_objects as go | |
import torch | |
from transformers import pipeline, TapasTokenizer, TapasForQuestionAnswering | |
st.set_page_config( | |
page_title="Sales Forecasting System", | |
page_icon="π", | |
layout="wide", | |
initial_sidebar_state="expanded", | |
) | |
# Preprocessing | |
def merge(B, C, A): | |
i = j = k = 0 | |
# Convert 'Date' columns to datetime.date objects | |
B['Date'] = pd.to_datetime(B['Date']).dt.date | |
C['Date'] = pd.to_datetime(C['Date']).dt.date | |
A['Date'] = pd.to_datetime(A['Date']).dt.date | |
while i < len(B) and j < len(C): | |
if B['Date'].iloc[i] <= C['Date'].iloc[j]: | |
A['Date'].iloc[k] = B['Date'].iloc[i] | |
A['Sales'].iloc[k] = B['Sales'].iloc[i] | |
i += 1 | |
else: | |
A['Date'].iloc[k] = C['Date'].iloc[j] | |
A['Sales'].iloc[k] = C['Sales'].iloc[j] | |
j += 1 | |
k += 1 | |
while i < len(B): | |
A['Date'].iloc[k] = B['Date'].iloc[i] | |
A['Sales'].iloc[k] = B['Sales'].iloc[i] | |
i += 1 | |
k += 1 | |
while j < len(C): | |
A['Date'].iloc[k] = C['Date'].iloc[j] | |
A['Sales'].iloc[k] = C['Sales'].iloc[j] | |
j += 1 | |
k += 1 | |
return A | |
def merge_sort(dataframe): | |
if len(dataframe) > 1: | |
center = len(dataframe) // 2 | |
left = dataframe.iloc[:center] | |
right = dataframe.iloc[center:] | |
merge_sort(left) | |
merge_sort(right) | |
return merge(left, right, dataframe) | |
else: | |
return dataframe | |
def drop (dataframe): | |
def get_columns_containing(dataframe, substrings): | |
return [col for col in dataframe.columns if any(substring.lower() in col.lower() for substring in substrings)] | |
columns_to_keep = get_columns_containing(dataframe, ["date", "sale"]) | |
dataframe = dataframe.drop(columns=dataframe.columns.difference(columns_to_keep)) | |
dataframe = dataframe.dropna() | |
return dataframe | |
def date_format(dataframe): | |
for i, d, s in dataframe.itertuples(): | |
dataframe['Date'][i] = dataframe['Date'][i].strip() | |
for i, d, s in dataframe.itertuples(): | |
new_date = datetime.strptime(dataframe['Date'][i], "%m/%d/%Y").date() | |
dataframe['Date'][i] = new_date | |
return dataframe | |
def group_to_three(dataframe): | |
dataframe['Date'] = pd.to_datetime(dataframe['Date']) | |
dataframe = dataframe.groupby([pd.Grouper(key='Date', freq='3D')])['Sales'].mean().round(2) | |
dataframe = dataframe.replace(0, np.nan).dropna() | |
return dataframe | |
def series_to_df_exogenous(series): | |
dataframe = series.to_frame() | |
dataframe = dataframe.reset_index() | |
dataframe = dataframe.set_index('Date') | |
dataframe = dataframe.dropna() | |
# Create the eXogenous values | |
dataframe['Sales First Difference'] = dataframe['Sales'] - dataframe['Sales'].shift(1) | |
dataframe['Seasonal First Difference'] = dataframe['Sales'] - dataframe['Sales'].shift(12) | |
dataframe = dataframe.dropna() | |
return dataframe | |
def dates_df(dataframe): | |
dataframe = dataframe.reset_index() | |
dataframe['Date'] = dataframe['Date'].dt.strftime('%B %d, %Y') | |
dataframe[dataframe.columns] = dataframe[dataframe.columns].astype(str) | |
return dataframe | |
def get_forecast_period(period): | |
return round(period / 3) | |
# SARIMAX Model | |
def train_test(dataframe, n): | |
training_y = dataframe.iloc[:-n,0] | |
test_y = dataframe.iloc[-n:,0] | |
test_y_series = pd.Series(test_y, index=dataframe.iloc[-n:, 0].index) | |
training_X = dataframe.iloc[:-n,1:] | |
test_X = dataframe.iloc[-n:,1:] | |
future_X = dataframe.iloc[0:,1:] | |
return (training_y, test_y, test_y_series, training_X, test_X, future_X) | |
def test_fitting(dataframe, Exo, trainY): | |
trainTestModel = auto_arima(X = Exo, y = trainY, start_p=1, start_q=1, | |
test='adf',min_p=1,min_q=1, | |
max_p=3, max_q=3, m=12, | |
start_P=2, start_Q=2, seasonal=True, | |
d=None, D=1, trace=True, | |
error_action='ignore', | |
suppress_warnings=True, | |
stepwise=True, maxiter = 50) | |
model = trainTestModel | |
return model | |
def forecast_accuracy(forecast, actual): | |
mape = np.mean(np.abs(forecast - actual)/np.abs(actual)).round(4) # MAPE | |
rmse = (np.mean((forecast - actual)**2)**.5).round(2) # RMSE | |
corr = np.corrcoef(forecast, actual)[0,1] # corr | |
mins = np.amin(np.hstack([forecast[:,None], | |
actual[:,None]]), axis=1) | |
maxs = np.amax(np.hstack([forecast[:,None], | |
actual[:,None]]), axis=1) | |
minmax = 1 - np.mean(mins/maxs) # minmax | |
return({'mape':mape, 'rmse':rmse, 'corr':corr, 'min-max':minmax}) | |
def sales_growth(dataframe, fittedValues): | |
sales_growth = fittedValues.to_frame() | |
sales_growth = sales_growth.reset_index() | |
sales_growth.columns = ("Date", "Sales") | |
sales_growth = sales_growth.set_index('Date') | |
sales_growth['Sales'] = (sales_growth['Sales']).round(2) | |
# Calculate and create the column for sales difference and growth | |
sales_growth['Forecasted Sales First Difference']=(sales_growth['Sales']-sales_growth['Sales'].shift(1)).round(2) | |
sales_growth['Forecasted Sales Growth']=(((sales_growth['Sales']-sales_growth['Sales'].shift(1))/sales_growth['Sales'].shift(1))*100).round(2) | |
# Calculate and create the first row for sales difference and growth | |
sales_growth['Forecasted Sales First Difference'].iloc[0] = (dataframe['Sales'].iloc[-1]-dataframe['Sales'].iloc[-2]).round(2) | |
sales_growth['Forecasted Sales Growth'].iloc[0]=(((dataframe['Sales'].iloc[-1]-dataframe['Sales'].iloc[-2])/dataframe['Sales'].iloc[-1])*100).round(2) | |
return sales_growth | |
def merge_forecast_data(actual, predicted, future): # debug | |
actual = actual.to_frame() | |
print("BEFORE RENAME ACTUAL") | |
print(actual) | |
actual.rename(columns={actual.columns[0]: "Actual Sales"}, inplace=True) | |
print("ACTUAL") | |
print(actual) | |
predicted = predicted.to_frame() | |
predicted.rename(columns={predicted.columns[0]: "Predicted Sales"}, inplace=True) | |
print("PREDICTED") | |
print(predicted) | |
future = future.to_frame() | |
future = future.rename_axis('Date') | |
future.rename(columns={future.columns[0]: "Forecasted Future Sales"}, inplace=True) | |
print("FUTURE") | |
print(future) | |
merged_dataframe = pd.concat([actual, predicted, future], axis=1) | |
print("MERGED DATAFRAME") | |
print(merged_dataframe) | |
merged_dataframe = merged_dataframe.reset_index() | |
print("MERGED DATAFRAME RESET INDEX") | |
print(merged_dataframe) | |
return merged_dataframe | |
def interpret_mape(mape_score): | |
score = (mape_score * 100).round(2) | |
if score < 10: | |
interpretation = "Great" | |
color = "green" | |
elif score < 20: | |
interpretation = "Good" | |
color = "seagreen" | |
elif score < 50: | |
interpretation = "Relatively good" | |
color = "orange" | |
else: | |
interpretation = "Poor" | |
color = "red" | |
return score, interpretation, color | |
# TAPAS Model | |
def load_tapas_model(): | |
model_name = "google/tapas-large-finetuned-wtq" | |
tokenizer = TapasTokenizer.from_pretrained(model_name) | |
model = TapasForQuestionAnswering.from_pretrained(model_name, local_files_only=False) | |
pipe = pipeline("table-question-answering", model=model, tokenizer=tokenizer) | |
return pipe | |
pipe = load_tapas_model() | |
def get_answer(table, query): | |
answers = pipe(table=table, query=query) | |
return answers | |
def convert_answer(answer): | |
if answer['aggregator'] == 'SUM': | |
cells = answer['cells'] | |
converted = sum(float(value.replace(',', '')) for value in cells) | |
return converted | |
if answer['aggregator'] == 'AVERAGE': | |
cells = answer['cells'] | |
values = [float(value.replace(',', '')) for value in cells] | |
converted = sum(values) / len(values) | |
return converted | |
if answer['aggregator'] == 'COUNT': | |
cells = answer['cells'] | |
converted = sum(int(value.replace(',', '')) for value in cells) | |
return converted | |
else: | |
return answer['answer'] | |
def get_converted_answer(table, query): | |
converted_answer = convert_answer(get_answer(table, query)) | |
return converted_answer | |
# Session States | |
if 'uploaded' not in st.session_state: | |
st.session_state.uploaded = False | |
if 'forecasted' not in st.session_state: | |
st.session_state.forecasted = False | |
# Web Application | |
st.title("Forecasting Dashboard π") | |
if not st.session_state.uploaded: | |
st.subheader("Welcome User, get started forecasting by uploading your file in the sidebar!") | |
# Sidebar Menu | |
with st.sidebar: | |
# TODO Name for product | |
st.title("MLCast v1.1") | |
st.subheader("An intelligent sales forecasting system") | |
uploaded_file = st.file_uploader("Upload your store data here to proceed (must atleast contain Date and Sales)", type=["csv"]) | |
if uploaded_file is not None: | |
date_found = False | |
sales_found = False | |
df = pd.read_csv(uploaded_file, parse_dates=True) | |
for column in df.columns: | |
if 'Date' in column: | |
date_found = True | |
if 'Sales' in column: | |
sales_found = True | |
if(date_found == False or sales_found == False): | |
st.error('Please upload a csv containing both Date and Sales...') | |
st.stop() | |
st.success("File uploaded successfully!") | |
st.write("Your uploaded data:") | |
st.write(df) | |
df = drop(df) | |
df = date_format(df) | |
merge_sort(df) | |
series = group_to_three(df) | |
st.session_state.uploaded = True | |
with open('sample.csv', 'rb') as f: | |
st.download_button("Download our sample CSV", f, file_name='sample.csv') | |
if (st.session_state.uploaded): | |
st.subheader("Sales History") | |
st.line_chart(series) | |
MIN_DAYS = 30 | |
MAX_DAYS = 90 | |
period = st.slider('How many days would you like to forecast?', min_value=MIN_DAYS, max_value=MAX_DAYS) | |
forecast_period = get_forecast_period(period) | |
forecast_button = st.button( | |
'Start Forecasting', | |
key='forecast_button', | |
type="primary", | |
) | |
if (forecast_button or st.session_state.forecasted): | |
df = series_to_df_exogenous(series) | |
n_periods = round(len(df) * 0.2) | |
print(n_periods) # debug | |
train = train_test(df, n_periods) | |
training_y, test_y, test_y_series, training_X, test_X, future_X = train | |
train_test_model = test_fitting(df, training_X, training_y) | |
print(df) # debug | |
print(len(df)) # debug | |
future_n_periods = forecast_period | |
fitted, confint = train_test_model.predict(X=test_X, n_periods=n_periods, return_conf_int=True) | |
index_of_fc = test_y_series.index | |
# make series for plotting purpose | |
fitted_series = pd.Series(fitted) | |
fitted_series.index = index_of_fc | |
lower_series = pd.Series(confint[:, 0], index=index_of_fc) | |
upper_series = pd.Series(confint[:, 1], index=index_of_fc) | |
#Future predictions | |
frequency = '3D' | |
future_fitted, confint = train_test_model.predict(X=df.iloc[-future_n_periods:,1:], n_periods=future_n_periods, return_conf_int=True, freq=frequency) | |
future_index_of_fc = pd.date_range(df['Sales'].index[-1], periods = future_n_periods, freq=frequency) | |
# make series for future plotting purpose | |
future_fitted_series = pd.Series(future_fitted) | |
future_fitted_series.index = future_index_of_fc | |
# future_lower_series = pd.Series(confint[:, 0], index=future_index_of_fc) | |
# future_upper_series = pd.Series(confint[:, 1], index=future_index_of_fc) | |
future_sales_growth = sales_growth(df, future_fitted_series) | |
test_y, predictions = np.array(test_y), np.array(fitted) | |
print("Test Y:", test_y) # debug | |
print("Prediction:", fitted) # debug | |
score = forecast_accuracy(predictions, test_y) | |
print("Score:", score) # debug | |
mape, interpretation, mape_color = interpret_mape(score['mape']) | |
print(df) | |
print(df['Sales']) | |
merged_data = merge_forecast_data(df['Sales'], fitted_series, future_fitted_series) | |
col_charts = st.columns(2) | |
print(merged_data) # debug | |
print(merged_data.info) | |
print(merged_data.dtypes) | |
with col_charts[0]: | |
fig_compare = go.Figure() | |
fig_compare.add_trace(go.Scatter(x=merged_data[merged_data.columns[0]], y=merged_data['Actual Sales'], mode='lines', name='Actual Sales')) | |
fig_compare.add_trace(go.Scatter(x=merged_data[merged_data.columns[0]], y=merged_data['Predicted Sales'], mode='lines', name='Predicted Sales', line=dict(color='#006400'))) | |
fig_compare.update_layout(title='Historical Sales Data', xaxis_title='Date', yaxis_title='Sales') | |
st.plotly_chart(fig_compare, use_container_width=True) | |
with col_charts[1]: | |
fig_forecast = go.Figure() | |
fig_forecast.add_trace(go.Scatter(x=merged_data[merged_data.columns[0]], y=merged_data['Actual Sales'], mode='lines', name='Actual Sales')) | |
fig_forecast.add_trace(go.Scatter(x=merged_data[merged_data.columns[0]], y=merged_data['Forecasted Future Sales'], mode='lines', name='Future Forecasted Sales', line=dict(color=mape_color))) | |
fig_forecast.update_layout(title='Forecasted Sales Data', xaxis_title='Date', yaxis_title='Sales') | |
st.plotly_chart(fig_forecast, use_container_width=True) | |
st.write(f"MAPE score: {mape}% - {interpretation}") | |
df = dates_df(future_sales_growth) | |
col_table = st.columns(2) | |
with col_table[0]: | |
col_table[0].subheader(f"Forecasted sales in the next {period} days") | |
col_table[0].write(df) | |
with col_table[1]: | |
col_table[1] = st.subheader("Question-Answering") | |
with st.form("question_form"): | |
question = st.text_input('Ask a Question about the Forecasted Data', placeholder="What is the total sales in the month of December?") | |
query_button = st.form_submit_button(label='Generate Answer') | |
if query_button or question: | |
answer = get_converted_answer(df, question) | |
if answer is not None: | |
st.write("The answer is:", answer) | |
else: | |
st.write("Answer is not found in table") | |
st.session_state.forecasted = True | |
# Hide Streamlit default style | |
hide_st_style = """ | |
<style> | |
footer {visibility: hidden;} | |
</style> | |
""" | |
st.markdown(hide_st_style, unsafe_allow_html=True) |