Expand description
PLINK1 --homozyg: runs of homozygosity via the scanning-window algorithm.
A window of window_snp consecutive variants slides one variant at a time.
A window is a hit when it holds at most window_het heterozygous and at most
window_missing missing calls. Each variant is “in a homozygous segment” when
the fraction of windows covering it that are hits is at least window_threshold.
Maximal stretches of in-segment variants on a chromosome become candidate runs,
split wherever consecutive variants lie more than gap kb apart, then kept only
when they clear snp, kb, and density (kb-per-variant ceiling).
Structs§
- Indiv
Row - Params
- Roh
- Sample
Roh - Summary
Row - Per-variant case/control tallies for the
.hom.summaryreport, in.bimorder.
Functions§
- homozyg
- Scan every sample for runs of homozygosity.
- indiv_
rows - summary_
rows - write_
hom - write_
indiv - write_
summary