Histogram

# === IMPORTS ===
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns

# === USER-EDITABLE PARAMETERS ===
x_col = 'Age'  # Change: Column name for x-axis
bins = 20  # Change: Number of histogram bins
figsize = (10, 6)  # Change: Figure size (width, height) in inches
title_fontsize = 18  # Change: Font size for the plot title
label_fontsize = 16  # Change: Font size for axis labels and ticks
bar_color = '#1f77b4'  # Change: Hex color for histogram bars
kde_color = '#ff7f0e'  # Change: Hex color for KDE line
kde_linewidth = 2  # Change: Line width for KDE
mean_line_color = '#d62728'  # Change: Hex color for mean vertical line
mean_line_style = '--'  # Change: Line style for mean line
mean_line_width = 1.5  # Change: Line width for mean line

# === DATA PREPARATION ===
# Generate example dataset
np.random.seed(0)
ages = np.random.normal(loc=35, scale=12, size=1000)
ages = np.clip(ages, a_min=0, a_max=None)  # Ensure ages are non-negative
df = pd.DataFrame({x_col: ages})

# Compute statistics
mean_age = df[x_col].mean()
std_age = df[x_col].std()
n_points = len(df)

# Print relevant data values
print(f"Number of data points: {n_points}")
print(f"Mean {x_col}: {mean_age:.1f}")
print(f"Standard deviation of {x_col}: {std_age:.1f}")

# === CREATE PLOT ===
sns.set_style('whitegrid')
fig, ax = plt.subplots(figsize=figsize)

# Histogram without KDE
sns.histplot(
    data=df,
    x=x_col,
    bins=bins,
    kde=False,
    color=bar_color,
    ax=ax
)

# KDE overlay
sns.kdeplot(
    data=df,
    x=x_col,
    color=kde_color,
    linewidth=kde_linewidth,
    ax=ax
)

# Add vertical line at mean
ax.axvline(mean_age, color=mean_line_color, linestyle=mean_line_style, linewidth=mean_line_width)

# Title with insight
ax.set_title(
    f"Age distribution: mean {mean_age:.1f} years, SD {std_age:.1f} years",
    fontsize=title_fontsize
)

# Axis labels
ax.set_xlabel(f"{x_col} (years)", fontsize=label_fontsize)
ax.set_ylabel("Count", fontsize=label_fontsize)

# Tick label sizes
ax.tick_params(labelsize=label_fontsize)

# Annotation box with statistics
ax.text(
    0.95,
    0.85,
    f"Mean: {mean_age:.1f}\nSD: {std_age:.1f}",
    transform=ax.transAxes,
    ha='right',
    va='top',
    bbox=dict(facecolor='white', alpha=0.8)
)

# Final layout adjustments
plt.tight_layout()
plt.show()
# END-OF-CODE

Statistics:
Mean: 41.8 years
Median: 41.2 years
Std Dev: 14.9 years
Min: 18.1 years
Max: 79.8 years
Skewness: 0.23 (slight right skew)