gen.simuHaplo_IBD_compare: Compare proband haplotypes for IBD sharing

In GENLIB: Genealogical Data Analysis

Compare proband haplotypes for IBD sharing

Description

Analyze the simulated haplotypes produced by gen.simuHaplo and identify segments of IBD sharing between the diploid haplotypes of a pair of probands

Usage

gen.simuHaplo_IBD_compare(proID_1, proID_2, BP_len, proband_haplotypes_path)

Arguments

proID_1	integer ID of a proband
proID_2	integer ID of the second proband
BP_len	The chromosome BP length used in the simulation
proband_haplotypes_path	String that gives the location of the "Proband_haplotypes.txt" output file produced by gen.simuHaplo

Value

returns a dataframe with the following columns: "simulNo", "n_seg", "pIBD", "mean_seg_len". Each row is a simulation wherein the pair of probands had some IBD sharing. The "simulNo" columns gives the simulation number, "n_seg" is the number of non-contiguous segments shared between the pair (divide this value by 2 for an average value for the pair), "pIBD" is the percent of the simulated chromosome shared IBD between the pair, and "mean_seg_len" is the mean length of IBD segments. Additionally, the function will print the exact location of each IBD segment to the R console. This can be toggled by disabling messages.

See Also

gen.genealogy gen.simuHaplo gen.simuHaplo_traceback

Examples

gen_table <- data.frame(ind=c(1,2,3,4,5,6,7,8),
                        mother=c(0,0,2,0,0,5,3,3),
                        father=c(0,0,1,0,0,4,6,6),
                        sex=c(1,2,2,1,2,1,1,2)
)

#convert the genealogy table into a GENLIB object
sample_gen<-gen.genealogy(gen_table) 

gen.graph(sample_gen)
# Simulate the gene dropping from founders to probamds

gen.simuHaplo(sample_gen, simulNo=5, model = 1, model_params = c(1,1), cM_len=c(100,100),
                BP_len = 100000000, outDir=tempdir())

IBD_compare <- gen.simuHaplo_IBD_compare(7, 8, 100000000, 
                paste(tempdir(), "Proband_Haplotypes.txt", sep="/"))

GENLIB documentation built on Jan. 17, 2023, 5:16 p.m.