Home
last modified time | relevance | path

Searched refs:u31tod (Results 1 – 17 of 17) sorted by relevance

/dports/biology/plink/plink-ng-79b2df8c/2.0/include/
H A Dplink2_stats.cc520 const double aa = u31tod(df) * 0.5; in gamma_incomplete_imp2()
757 const double aa = u31tod(df) * 0.5; in gamma_incomplete_imp2_ln()
1146 const double n_d = u31tod(nn); in binomial_ccdf_ln()
1324 const double nn_d = u31tod(nn); in beta_small_b_large_a_series_ln()
1481 double aa = u31tod(df1) * 0.5; in ibeta_imp2_ln()
1501 double bb = u31tod(df2) * 0.5; in ibeta_imp2_ln()
1571 const double bbar = bb - u31tod(nn); in ibeta_imp2_ln()
1593 return u31tod(do_invert); in TstatToP2()
1607 const double df_d = u31tod(df); in TstatToLnP()
1622 const double df1_d = u31tod(df1); in FstatToLnP()
[all …]
H A Dplink2_string.cc2279 return result + S_CAST(int32_t, (dxx - u31tod(result)) + banker_round[result & 1]);
2287 remainder += S_CAST(int32_t, (dxx - u31tod(remainder)) + banker_round[remainder & 1]);
2295 remainder += S_CAST(int32_t, (dxx - u31tod(remainder)) + banker_round[remainder & 1]);
2303 remainder += S_CAST(int32_t, (dxx - u31tod(remainder)) + banker_round[remainder & 1]);
2311 remainder += S_CAST(int32_t, (dxx - u31tod(remainder)) + banker_round[remainder & 1]);
2319 remainder += S_CAST(int32_t, (dxx - u31tod(remainder)) + banker_round[remainder & 1]);
2327 remainder += S_CAST(int32_t, (dxx - u31tod(remainder)) + banker_round[remainder & 1]);
2335 remainder += S_CAST(int32_t, (dxx - u31tod(remainder)) + banker_round[remainder & 1]);
2845 if (fabs(dxx_10_6 - u31tod(dec_digits)) >= 0.00000005) { in dtoa_f_probp6_spaced()
2858 if (fabs(dxx_10_6 - u31tod(dec_digits)) >= 0.00000005) { in dtoa_f_probp6_clipped()
H A Dpgenlib_misc.cc2760 …ed_variance_times_2n = (observed_variance_times_2n_hi * 4294967296.0) + (u31tod(extra_phased_het_c… in MultiallelicDiploidMinimac3R2()
H A Dplink2_base.h270 HEADER_INLINE double u31tod(uint32_t uii) { in u31tod() function
/dports/biology/plink/plink-ng-79b2df8c/2.0/
H A Dpgenlib_ffi_support.cc202 const double missing_val = u63tod(numer) / u31tod(denom); in Dosage16ToDoublesMeanimpute()
289 const double denom = u31tod(sample_ct - genocounts[3]); in LinearCombinationMeanimpute()
348 resultx += u31tod(cur_dosage) * weights[sample_uidx]; in LinearCombinationMeanimpute()
352 const double denom = 16384 * u31tod(sample_ct - missing_ct); in LinearCombinationMeanimpute()
H A Dplink2_adjust.cc355 const double harmonic_term = valid_allele_ctd / u31tod(aidx); in Multcomp()
372 const double by_ln_pval = sorted_ln_pvals[aidx - 1] - log(u31tod(aidx)) + ln_harmonic_sum; in Multcomp()
474 const double ln_pv_holm_new = ln_pval + log(u31tod(valid_allele_ct - aidx)); in Multcomp()
H A Dplink2_glm.cc1267 const double sample_ct_recip = 1.0 / u31tod(sample_ct); in CheckMaxCorrAndVif()
1268 const double sample_ct_m1_d = u31tod(sample_ct - 1); in CheckMaxCorrAndVif()
1459 const double sample_ct_recip = 1.0 / u31tod(sample_ct); in CheckMaxCorrAndVifF()
1460 const double sample_ct_m1_d = u31tod(sample_ct - 1); in CheckMaxCorrAndVifF()
1548 double max_ssq = 1.5 * u31tod(sample_ct); in GlmFillAndTestCovars()
1549 double min_ssq_minus_sqmean = 0.0198 * u31tod(sample_ct - 1); in GlmFillAndTestCovars()
1576 covar_ssq -= covar_sum * covar_sum / u31tod(sample_ct); in GlmFillAndTestCovars()
1614 *covar_write_iter++ = u31tod(cur_sample_is_in_cat); in GlmFillAndTestCovars()
1616 double covar_ssq = u31tod(cur_cat_obs_ct); in GlmFillAndTestCovars()
1659 xtx_image[0] = u31tod(sample_ct); in InitNmPrecomp()
[all …]
H A Dplink2_misc.cc2684 const double write_val = u31tod(1 + qt_12); in SplitCatPheno()
2838 const double cur_shifted_mean = shifted_pheno_sum / u31tod(cur_sample_ct); in PhenoVarianceStandardize()
4195 variant_ct_recips[0] = 1.0 / u31tod(variant_ct_nony); in WriteMissingnessReports()
4200 variant_ct_recips[1] = 1.0 / u31tod(variant_ct); in WriteMissingnessReports()
4383 nobs_recip = 1.0 / u31tod(cur_nobs); in WriteMissingnessReports()
4451 cswritep = dtoa_g(u31tod(cur_missing_hc_ct) * nobs_recip, cswritep); in WriteMissingnessReports()
4459 cswritep = dtoa_g(u31tod(cur_hethap_ct) * nobs_recip, cswritep); in WriteMissingnessReports()
5115 const double nonmissing_ct_recip = 1.0 / u31tod(nonmissing_ct);
5116 cswritep = dtoa_g(u31tod(het_a1_ct) * nonmissing_ct_recip, cswritep);
5342 const double female_obs_ct_recip = 1.0 / u31tod(female_obs_ct);
[all …]
H A Dplink2_matrix_calc.cc1986 nonmiss_recip = 1.0 / u31tod(nonmiss_ct); in CalcKing()
1993 cswritetp = dtoa_g(nonmiss_recip * u31tod(hethet_ct), cswritetp); in CalcKing()
2001 cswritetp = dtoa_g(nonmiss_recip * u31tod(ibs0_ct), cswritetp); in CalcKing()
3113 nonmiss_recip = 1.0 / u31tod(nonmiss_ct); in CalcKingTableSubset()
3120 cswritep = dtoa_g(nonmiss_recip * u31tod(hethet_ct), cswritep); in CalcKingTableSubset()
3128 cswritep = dtoa_g(nonmiss_recip * u31tod(ibs0_ct), cswritep); in CalcKingTableSubset()
3137 cswritep = dtoa_g(nonmiss_recip * u31tod(het1hom2_ct), cswritep); in CalcKingTableSubset()
3139 cswritep = dtoa_g(nonmiss_recip * u31tod(het2hom1_ct), cswritep); in CalcKingTableSubset()
4157 *grm_iter++ /= u31tod(variant_ct_base); in CalcGrm()
4160 const double variant_ct_recip = 1.0 / u31tod(variant_ct); in CalcGrm()
[all …]
H A Dplink2_ld.cc2991 nmajsums_d[0] = u31tod(nmaj_cts[0]); in LdConsole()
2992 nmajsums_d[1] = u31tod(nmaj_cts[1]); in LdConsole()
2994 unknown_hethet_d = u31tod(unknown_hethet_ct); in LdConsole()
3012 x_male_nmajsums_d[0] = u31tod(x_male_nmaj_cts[0]); in LdConsole()
3013 x_male_nmajsums_d[1] = u31tod(x_male_nmaj_cts[1]); in LdConsole()
3137 double valid_obs_d = u31tod(valid_obs_ct); in LdConsole()
3145 valid_obs_d -= male_decr * u31tod(valid_x_male_ct); in LdConsole()
H A Dplink2_filter.cc3584 max_missing_cts[0] = S_CAST(int32_t, u31tod(variant_ct - variant_ct_y) * mind_thresh); in MindFilter()
3585 max_missing_cts[1] = S_CAST(int32_t, u31tod(variant_ct) * mind_thresh); in MindFilter()
3631 const uint32_t missing_max_ct_nony = S_CAST(int32_t, geno_thresh * u31tod(sample_ct)); in EnforceGenoThresh()
3632 const uint32_t missing_max_ct_y = S_CAST(int32_t, geno_thresh * u31tod(male_ct)); in EnforceGenoThresh()
H A Dplink2_psam.cc335 const double pheno_ctrld = u31tod(1 - affection_01); in LoadPsam()
H A Dplink2.cc583 …od(nony_possible_obs - tot_nony_missing) / u31tod(sample_ct) + u63tod(y_possible_obs - tot_y_missi… in ReportGenotypingRate()
584 genotyping_rate /= u31tod(variant_ct); in ReportGenotypingRate()
5759 … if (u31tod(hard_call_thresh) + import_dosage_certainty * kDosageMid >= u31tod(kDosageMid)) {
H A Dplink2_export.cc3759 const double nonmale_geno_ct_recip = 1.0 / u31tod(variant_ct - y_ct); in ExportOxSample()
3760 const double male_geno_ct_recip = 1.0 / u31tod(variant_ct); in ExportOxSample()
H A Dplink2_data.cc1971 const double dosage_avg = dosage_sumd / u31tod(obs_ct); in LoadAlleleAndGenoCountsThread()
H A Dplink2_import.cc12520 const double denom_d = u31tod(denom); in OxBgenToPgen()
/dports/biology/plink/plink-ng-79b2df8c/2.0/pgenlibr/src/
H A Dpgenlib_ffi_support.cpp202 const double missing_val = u63tod(numer) / u31tod(denom); in Dosage16ToDoublesMeanimpute()
289 const double denom = u31tod(sample_ct - genocounts[3]); in LinearCombinationMeanimpute()
348 resultx += u31tod(cur_dosage) * weights[sample_uidx]; in LinearCombinationMeanimpute()
352 const double denom = 16384 * u31tod(sample_ct - missing_ct); in LinearCombinationMeanimpute()