Online shopping decisions rely on how consumers engage with online store content. You work for a new startup company that has just launched a new online shopping website. The marketing team asks you, a new data scientist, to review a dataset of online shoppers' purchasing intentions gathered over the last year. Specifically, the team wants you to generate some insights into customer browsing behaviors in November and December, the busiest months for shoppers. You have decided to identify two groups of customers: those with a low purchase rate and returning customers. After identifying these groups, you want to determine the probability that any of these customers will make a purchase in a new marketing campaign to help gauge potential success for next year's sales.
Data description:
You are given an online_shopping_session_data.csv that contains several columns about each shopping session. Each shopping session corresponded to a single user.
| Column | Description |
|---|---|
SessionID | unique session ID |
Administrative | number of pages visited related to the customer account |
Administrative_Duration | total amount of time spent (in seconds) on administrative pages |
Informational | number of pages visited related to the website and the company |
Informational_Duration | total amount of time spent (in seconds) on informational pages |
ProductRelated | number of pages visited related to available products |
ProductRelated_Duration | total amount of time spent (in seconds) on product-related pages |
BounceRates | average bounce rate of pages visited by the customer |
ExitRates | average exit rate of pages visited by the customer |
PageValues | average page value of pages visited by the customer |
SpecialDay | closeness of the site visiting time to a specific special day |
Weekend | indicator whether the session is on a weekend |
Month | month of the session date |
CustomerType | customer type |
Purchase | class label whether the customer make a purchase |
# Import required libraries
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import binom
# Load and view your data
shopping_data = pd.read_csv("online_shopping_session_data.csv")
shopping_data.head()# Start your code here!
nov_dec_shopping = shopping_data[shopping_data['Month'].isin(['Nov','Dec'])]
# Use as many cells as you likenov_dec_shopping['Month'].unique()# Start your code here!
purchase_rates = nov_dec_shopping.groupby(['CustomerType'])['Purchase'].mean()
# Sort the index in descending order
purchase_rates = purchase_rates.sort_index(ascending=False)
# Convert the Series to a dictionary
purchase_rates = purchase_rates.to_dict()
purchase_ratescorr_matrix = nov_dec_shopping[['Administrative_Duration', 'Informational_Duration', 'ProductRelated_Duration']].corr()
corr_matrix# Reset the index to convert the matrix into a DataFrame for easier processing
correlation_df = corr_matrix.stack().reset_index()
correlation_dfcorrelation_df.columns = ['Column1', 'Column2', 'Correlation']
correlation_df# Filter out duplicate correlations and self-correlations
correlation_df = correlation_df[correlation_df['Column1'] != correlation_df['Column2']]
correlation_df# Identify the row with the highest correlation
max_corr_row = correlation_df.loc[correlation_df['Correlation'].idxmax()]# Extract the pair and the correlation
top_correlation = {
"pair": (max_corr_row['Column1'], max_corr_row['Column2']),
"correlation": max_corr_row['Correlation']
}
print(top_correlation)returning_customer_boost = 1.15*purchase_rates['Returning_Customer']
returning_customer_boostprob_at_least_100_sales = 1 - binom.cdf(100, 500, returning_customer_boost)
prob_at_least_100_sales# Parameters
n = 500 # Number of trials
p = returning_customer_boost # Probability of success
x = np.arange(0, n + 1) # Possible values of successes
pmf = binom.pmf(x, n, p) # Binomial PMF
# Plot
plt.figure(figsize=(8, 5))
plt.bar(x, pmf, color='skyblue', alpha=0.7, label='Binomial PMF')
plt.vlines(100, 0, 0.08, color='r', linestyle='dashed', label="sales=100")
plt.xlabel('Number of Successes (k)')
plt.ylabel('Probability')
plt.title(f'Binomial Distribution (n={n}, p={p})')
plt.legend()
plt.show()