for l in range(len(DNAvals)):
        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(Tvals[not_nan], ln_mu_LxT[l][not_nan])