def false_discovery(alpha,beta,rho):
    """The false discovery rate.
    
    The false discovery rate is the probability that an observed edge is
    incorrectly identified, namely that is doesn't exist in the 'true' network.
    This is one measure of how reliable the results are.
    
    Parameters
    ----------
    alpha : float
        The estimate of the true-positive rate.
    beta : float
        The estimate of the false-positive rate.
    rho : float
        The estimate of network density.
    
    Returns
    -------
    float
        The false discovery rate (probability).
        
    References
    ----------
    .. [1] Newman, M.E.J. 2018. “Network structure from rich but noisy data.” 
        Nature Physics 14 6 (June 1): 542–545. doi:10.1038/s41567-018-0076-1.
    """
    return (1-rho)*beta/(rho*alpha + (1-rho)*beta)