def samplets(features, confounds, x_var=None, y_var=None, group_by=None,
    color_by=None, sort_x_by=None):
    """
    Plot of the feature distribution for each samplet
    along a axis of a given variable,
    grouped and coloured by chosen variables.
    """
    raise NotImplementedError()
