for l in range(len(DNAspace_for_mu0)):
            not_nan = _np.logical_not(_np.isnan(ln_mu_LxT[l]))
            if sum(not_nan) > 1:
                # (enough points for linear regression)
                gradient, intercept, r_value, p_value, std_err = _stats.linregress(