Seaborn Joint Hexplot with marginal KDE + hists.

def joint_hex(x, y, **kwargs):
    """Seaborn Joint Hexplot with marginal KDE + hists."""
    return sns.jointplot(
        x, y, kind='hex', stat_func=None, marginal_kws={'kde': True}, **kwargs
    )