def chooseKpos(Lsca, Lrand, stddev=2):
    """Calculate the significant number of eigenvalues.

    Parameters
    ----------
    Lsca : :class:`numpy.array`
        Eigenvector from coupling strengths
    Lrand : :class:`numpy.array`
        Matrix of eigenvalues from randomized coupling strengths
    stddev : int, optional
        Number of standard deviations to use

    Returns
    -------
    Number of significant eigenvalues
    """
    value = Lrand[:, 1].mean() + ((stddev + 1) * Lrand[:, 1].std())
    return Lsca[Lsca > value].shape[0]