Jointplot

import seaborn as sns
sns.set_theme(style="white")

penguins = sns.load_dataset("penguins")
sns.jointplot(data=penguins, x="bill_length_mm", y="bill_depth_mm")

sns.jointplot(data=penguins, x="bill_length_mm", y="bill_depth_mm", hue="species")

sns.jointplot(data=penguins, x="bill_length_mm", y="bill_depth_mm", hue="species", kind="kde")

sns.jointplot(data=penguins, x="bill_length_mm", y="bill_depth_mm", kind="reg")

sns.jointplot(data=penguins, x="bill_length_mm", y="bill_depth_mm", kind="hist")

sns.jointplot(data=penguins, x="bill_length_mm", y="bill_depth_mm", kind="hex")

sns.jointplot(
    data=penguins, x="bill_length_mm", y="bill_depth_mm",
    marker="+", s=100, marginal_kws=dict(bins=25, fill=False),
)

sns.jointplot(data=penguins, x="bill_length_mm", y="bill_depth_mm", height=5, ratio=2, marginal_ticks=True)

g = sns.jointplot(data=penguins, x="bill_length_mm", y="bill_depth_mm")
g.plot_joint(sns.kdeplot, color="r", zorder=0, levels=6)
g.plot_marginals(sns.rugplot, color="r", height=-.15, clip_on=False)