admix.simulate.binary_pheno#

admix.simulate.binary_pheno(dset: Dataset, hsq: float = 0.5, cor: float | None = None, n_causal: int | None = None, case_prevalence: float = 0.1, beta: ndarray | None = None, cov_cols: List[str] | None = None, cov_effects: List[float] | None = None, n_sim=10, method='probit')[source]#

Simulate under liability threshold. First simulate quantative traits in the liability scale. With the given case_prevalence, convert liability to binary traits.

Parameters:
  • hsq (if method is “probit”, this corresponds to the proportion of variance explained)

  • by the genotype in the liability scale.

  • if method is “logit”, this corresponds to the variance explained of (Xeta),

  • and hsq can be larger than 1, in the case.