def get_shapes(n_features, n_responses):
    """
    Gets the shapes of the coefficeint and intercept.

    Parameters
    ----------
    n_features: int
        Number of features.

    n_responses: int
        Number of responses

    Output
    ------
    coef_shape, intercept_shape

    coef_shape: tuple
        Shape of the coefficient.

    intercept_shape: tuple
        Shape of the intercept.
    """
    if n_responses == 1:
        coef_shape = (n_features, )
        intercept_shape = ()
    else:
        coef_shape = (n_features, n_responses)
        intercept_shape = (n_responses, )

    return coef_shape, intercept_shape