def format_pvalue(p_value, alpha=0.05, include_equal=True):
    """
    If p-value is lower than 0.05, change it to "<0.05", otherwise, round it to two decimals
    :param p_val: input p-value as a float
    :param alpha: significance level
    :param include_equal: include equal sign ('=') to pvalue (e.g., '=0.06') or not (e.g., '0.06')
    :return: p_val: processed p-value (replaced by "<0.05" or rounded to two decimals) as a str
    """
    if p_value < alpha:
        p_value = "<" + str(alpha)
    else:
        if include_equal:
            p_value = '=' + str(round(p_value, 3))
        else:
            p_value = str(round(p_value, 3))

    return p_value