haplotype / kindred

infer realized kinship via latent IBD states
MIT License
1 stars 0 forks source link

Specifics of -j 3 option #1

Open slcapel opened 8 months ago

slcapel commented 8 months ago

Hello, I was hoping you could tell me which Jacquard states are implemented when setting -j 3. I also notice that when I implement -j 3, the last 4 columns of out.kin.gz do not sum to column 4.

haplotype commented 8 months ago

in table 1 of the manuscript (can be found here https://www.biorxiv.org/content/10.1101/2023.09.02.556031v1.full.pdf) -j 3 assume loadings from \Sigma 1-6 are 0 and only considers \Sigma 7,8,9, which are Jacquard states of no inbreeding.

Grant

On Wed, Dec 20, 2023 at 6:21 PM Samantha Capel @.***> wrote:

Hello, I was hoping you could tell me which Jacquard states are implemented when setting -j 3. Thank you.

— Reply to this email directly, view it on GitHub https://github.com/haplotype/kindred/issues/1, or unsubscribe https://github.com/notifications/unsubscribe-auth/ACWHWL2532MDFW4LMMUPID3YKNXILAVCNFSM6AAAAABA5RY6EWVHI2DSMVQWIX3LMV43ASLTON2WKOZSGA2TCNBUGMZDINY . You are receiving this because you are subscribed to this thread.Message ID: @.***>

slcapel commented 8 months ago

Thank you! Is the Supplement that pairs with the manuscript available? I can't seem to find it on bioRxiv.

haplotype commented 8 months ago

oh. I just uploaded both manuscript and supplementary here at github. the supplementary can be found here: https://github.com/haplotype/kindred/blob/main/kindred-supp.pdf

Grant

On Thu, Dec 21, 2023 at 3:42 PM Samantha Capel @.***> wrote:

Thank you! Is the Supplement that pairs with the manuscript available? I can't seem to find it on bioRxiv.

— Reply to this email directly, view it on GitHub https://github.com/haplotype/kindred/issues/1#issuecomment-1866905338, or unsubscribe https://github.com/notifications/unsubscribe-auth/ACWHWLYFNHKDMWMBYLWXR7TYKSNMVAVCNFSM6AAAAABA5RY6EWVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNRWHEYDKMZTHA . You are receiving this because you commented.Message ID: @.***>