R/falconer_threshold.R
qtl_threshold_chisq_mssn.RdCalculates the minimum selected case and control sample sizes after deriving their genotype distributions from a single-trait Falconer threshold model.
qtl_threshold_chisq_mssn(
power,
alpha,
qtl_var,
tau,
pd,
x_upper,
x_lower,
k = 1,
verbose = TRUE
)Target power in \((0,1)\).
Significance level in \((0,1)\).
Variance in the standardized trait attributable to the QTL. Must be in \((0,1)\).
Finite dominance-to-additivity ratio \(\delta/a\).
Frequency of the increaser allele. Must be in \((0,1)\).
Percentage selected from the upper population tail in percentile mode.
Percentage selected from the lower population tail in percentile mode.
Positive control-to-case ratio \(N_{control}/N_{case}\).
Logical. If TRUE, prints a concise summary.
An object of class "qtl_threshold_chisq_mssn" containing
selected MSSN values, target NCP, internal S, thresholds,
penetrances, prevalences, conditional frequencies, Falconer parameters,
and separately labelled expected population screening counts.
The primary MSSN is the selected case-control sample size used by the association test. The separately reported screening quantities are expected population counts needed to obtain those selected samples under simple population sampling; they are not statistical MSSN values. Genotype-specific distributions and threshold selection follow Eqs. 6.1 and 6.3–6.5 (pp. 324–326). The downstream genotype-test NCP is Eq. 1.22 (p. 26). The function solves for selected cases and controls; expected screening counts are reported separately and are not rounded MSSNs.
Gordon, D., Finch, S. J., & Kim, W. (2020). Heterogeneity in Statistical Genetics: How to Assess, Address, and Account for Mixtures in Association Studies. Springer, Eqs. 1.22, 6.1, and 6.3–6.5, pp. 26 and 324–326. doi:10.1007/978-3-030-61121-7 .
qtl_threshold_chisq_mssn(
power = 0.8, alpha = 0.0001,
qtl_var = 0.025, tau = 0.5, pd = 0.15,
x_upper = 5, x_lower = 5, verbose = FALSE
)