Step 1: splitSplit a VCF (or raw genotype/ancestry matrices) by local ancestry. ancestry_split() and the K-population split_*_multi() primitives work for any number of populations, not just two. |
|
|---|---|
Split a phased VCF by local ancestry (Step 1) |
|
Run ancestry splitting pipeline |
|
Run phased ancestry splitting pipeline |
|
Split genotype matrix by ancestry |
|
Split unphased genotype matrix by ancestry for K populations |
|
Split phased haplotypes by ancestry |
|
Split phased haplotypes into K population-specific dosage matrices |
|
Read a PLINK .bed file of BED-encoded local ancestry calls |
|
Step 2: write GDSMaterialize a Step 1 split result to per-population GDS files. |
|
Write all populations of an ancestry split to GDS files (Step 2) |
|
Write an ancestry-specific dosage matrix to a SeqArray GDS file |
|
Step 3: associateRun GMMAT::SMMAT() across ancestry-split GDS files and Cauchy-combine the per-gene p-values. |
|
Run ancestry-stratified SMMAT + Cauchy combination (Step 3) |
|
Cauchy combination test for multiple p-values |
|
Package |
|
lantern: Ancestry-Specific Rare Variant Analysis |
|