diff options
Diffstat (limited to 'src/lmm.cpp')
| -rw-r--r-- | src/lmm.cpp | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/lmm.cpp b/src/lmm.cpp index 80372ee..dbf1ad5 100644 --- a/src/lmm.cpp +++ b/src/lmm.cpp @@ -1144,7 +1144,7 @@ void LMM::CalcRLWald(const double &l, const FUNC_PARAM ¶ms, double &beta, beta = P_xy / P_xx; double tau = (double)df / Px_yy; - se = sqrt(1.0 / (tau * P_xx)); + se = safe_sqrt(1.0 / (tau * P_xx)); p_wald = gsl_cdf_fdist_Q((P_yy - Px_yy) * tau, 1.0, df); gsl_matrix_free(Pab); @@ -1186,7 +1186,7 @@ void LMM::CalcRLScore(const double &l, const FUNC_PARAM ¶ms, double &beta, beta = P_xy / P_xx; double tau = (double)df / Px_yy; - se = sqrt(1.0 / (tau * P_xx)); + se = safe_sqrt(1.0 / (tau * P_xx)); p_score = gsl_cdf_fdist_Q((double)ni_test * P_xy * P_xy / (P_yy * P_xx), 1.0, df); @@ -1201,7 +1201,7 @@ void CalcUab(const gsl_matrix *UtW, const gsl_vector *Uty, gsl_matrix *Uab) { size_t index_ab; size_t n_cvt = UtW->size2; - debug_msg("entering"); + // debug_msg("entering"); gsl_vector *u_a = gsl_vector_safe_alloc(Uty->size); @@ -2018,7 +2018,7 @@ void CalcPve(const gsl_vector *eval, const gsl_matrix *UtW, FUNC_PARAM param0 = {true, ni_test, n_cvt, eval, Uab, ab, 0}; - double se = sqrt(-1.0 / LogRL_dev2(lambda, ¶m0)); + double se = safe_sqrt(-1.0 / LogRL_dev2(lambda, ¶m0)); pve = trace_G * lambda / (trace_G * lambda + 1.0); pve_se = trace_G / ((trace_G * lambda + 1.0) * (trace_G * lambda + 1.0)) * se; @@ -2085,7 +2085,7 @@ void CalcLmmVgVeBeta(const gsl_vector *eval, const gsl_matrix *UtW, // Obtain se_beta. for (size_t i = 0; i < Vbeta->size1; i++) { - gsl_vector_set(se_beta, i, sqrt(gsl_matrix_get(Vbeta, i, i))); + gsl_vector_set(se_beta, i, safe_sqrt(gsl_matrix_get(Vbeta, i, i))); } gsl_matrix_safe_free(Uab); |