def get_tukey_pvalue(k, df, q):
    """
    return adjusted p-values for Tukey's HSD

    Parameters
    ----------
    k : int in {2, ..., 10}
        number of tests
    df : int
        degrees of freedom of error term
    q : scalar, array_like; q >= 0
        quantile value of Studentized Range

    """
    from statsmodels.stats.libqsturng import psturng
    return psturng(q, k, df)
