def refactor_confints_todict(means, st_devs, feat_names):
    """Refactor means and confidence intervals into a dictionary

    Args:
        means: list of the means of the WRR coefficients
        st_devs: list of the standard deviations of the WRR coefficients
        feat_names: list of feature names associated with the coefficients

    Returns:
        conf_int: dictionary,
            key = the feature name
            value = confidence interval for the feature (upper, lower bound)
    """

    conf_int = {}
    for name, mean, stdev in zip(feat_names, means, st_devs):
        conf_int[name] = [mean - 1.96 * stdev, mean + 1.96 * stdev]
    return conf_int