Searched refs:p_lrt (Results 1 – 4 of 4) sorted by relevance
| /dports/biology/gemma/GEMMA-0.98.3/src/ |
| H A D | lm.cpp | 268 double &beta, double &se, double &p_wald, double &p_lrt, in LmCalcP() argument
279 p_lrt = gsl_cdf_chisq_Q((double)n_size * (log(yPwy) - log(yPxy)), 1); in LmCalcP()
303 double beta = 0, se = 0, p_wald = 0, p_lrt = 0, p_score = 0; in AnalyzeGene() local
358 p_lrt, p_score); in AnalyzeGene()
363 SUMSTAT SNPs = {beta, se, 0.0, 0.0, p_wald, p_lrt, p_score, -0.0 }; in AnalyzeGene()
395 double beta = 0, se = 0, p_wald = 0, p_lrt = 0, p_score = 0; in AnalyzeBimbam() local
472 p_lrt, p_score); in AnalyzeBimbam()
477 SUMSTAT SNPs = {beta, se, 0.0, 0.0, p_wald, p_lrt, p_score, -0.0}; in AnalyzeBimbam()
511 double beta = 0, se = 0, p_wald = 0, p_lrt = 0, p_score = 0; in AnalyzePlink() local
616 p_lrt, p_score); in AnalyzePlink()
[all …]
|
| H A D | lmm.cpp | 160 << st.p_lrt << endl; in WriteFiles()
169 << st.p_lrt << "\t" in WriteFiles()
1371 double p_lrt = 0, p_score = 0; in AnalyzeGene() local
1441 p_lrt = gsl_cdf_chisq_Q(2.0 * (logl_H1 - logl_H0), 1); in AnalyzeGene()
1525 double p_lrt = 0.0, p_score = 0.0; in Analyze() local
1548 p_wald, p_lrt, p_score, logl_H1}; in Analyze()
1714 double p_lrt = 0, p_score = 0; in AnalyzePlink() local
1870 p_wald, p_lrt, p_score, logl_H1}; in AnalyzePlink()
2245 double p_lrt = 0, p_score = 0; in AnalyzeBimbamGXE() local
2348 p_lrt = gsl_cdf_chisq_Q(2.0 * (logl_H1 - logl_H0), 1); in AnalyzeBimbamGXE()
[all …]
|
| H A D | mvlmm.cpp | 195 outfile << sumStat[t].p_lrt << endl; in WriteFiles()
199 outfile << sumStat[t].p_wald << "\t" << sumStat[t].p_lrt << "\t" in WriteFiles()
3310 if (p_lrt < p_nr) { in AnalyzeBimbam()
3321 p_lrt = PCRT(2, d_size, p_lrt, crt_a, crt_b, crt_c); in AnalyzeBimbam()
3787 if (p_lrt < p_nr) { in AnalyzePlink()
3797 p_lrt = PCRT(2, d_size, p_lrt, crt_a, crt_b, crt_c); in AnalyzePlink()
4298 p_lrt = MphCalcP(eval, &X_row2.vector, &X_sub2.matrix, Y, V_g, V_e, in AnalyzeBimbamGXE()
4302 if (p_lrt < p_nr) { in AnalyzeBimbamGXE()
4313 p_lrt = PCRT(2, d_size, p_lrt, crt_a, crt_b, crt_c); in AnalyzeBimbamGXE()
4779 if (p_lrt < p_nr) { in AnalyzePlinkGXE()
[all …]
|
| H A D | param.h | 61 double p_lrt; // p value from a likelihood ratio test. variable
73 double p_lrt; // p value from a likelihood ratio test. variable
|