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

np.random.seed(42)
age_data = [
    1,3,3,3,5,6,6,6,7,8,8,9,9,9,9,9,
    10,10,10,10,10,10,11,11,11,11,11,11,
    12,12,12,12,12,12,13,13,13,13,13,
    14,14,14,14,14,14,15,15,15,15,15,15,15,15,
    16,16,16,16,17,17,17,17,17,17,18,18,18,18,
    19,19,19,20,21,22,22,24,24
]

# Create dataframe
ages = pd.DataFrame({'Age': age_data})

# Create dataframe
patients = pd.DataFrame({'Age': ages})

# Set seaborn style
sns.set_theme(style="white")

# Histogram of ages
plt.figure(figsize=(8,6))
sns.histplot(patients['Age'], bins=12, color="#D44803", edgecolor="black")

plt.title('Histogram: Patient Age Distribution with PikBioStat')
plt.xlabel('Age (years)')
plt.ylabel('Number of Patients')

sns.despine(left=True, bottom=True)
plt.show()