def standardize_snps(G):
    r"""
    Standardize variantes.

    Parameters
    ----------
    G : (`n_inds`, `n_snps`) array
        Genetic data

    Returns
    -------
    G_out : standardized array
    """
    mean = G.mean(0)
    std = G.std(0)
    return (G - mean) / std