def var_data(f, N):
    """
    Variance of MLE of beta for quantitative trait, assuming var(y) = 1

    Args:
      f: minor allele freq
      N: sample size

    Returns:
      variance of MLE beta
    """
    return 1 / (2 * N * f * (1 - f))