computePairIBD2: Computes pairwise IBD2 for a specific pair of individuals
In sim1000G: Genotype Simulations for Rare or Common Variants Using Haplotypes from 1000 Genomes
Description
Usage
Arguments
Value
Examples
Description
Computes pairwise IBD2 for a specific pair of individuals
Usage
1
computePairIBD2(i, j)
Arguments
i
Index of first individual
j
Index of second individual
Value
Mean IBD2 as computed from shared haplotypes
Examples
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
library("sim1000G")
examples_dir = system.file("examples", package = "sim1000G")
vcf_file = file.path(examples_dir, "region.vcf.gz")
vcf = readVCF( vcf_file, maxNumberOfVariants = 100 ,
min_maf = 0.12 ,max_maf = NA)
# For realistic data use the function downloadGeneticMap
generateUniformGeneticMap()
startSimulation(vcf, totalNumberOfIndividuals = 200)
ped1 = newNuclearFamily(1)
v = computePairIBD2(1, 3)
cat("IBD2 of pair = ", v, "\n");
sim1000G documentation built on June 10, 2019, 1:01 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Value Examples
Description
Computes pairwise IBD2 for a specific pair of individuals
Usage
1 | computePairIBD2(i, j)
|
Arguments
i |
Index of first individual |
j |
Index of second individual |
Value
Mean IBD2 as computed from shared haplotypes
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | library("sim1000G")
examples_dir = system.file("examples", package = "sim1000G")
vcf_file = file.path(examples_dir, "region.vcf.gz")
vcf = readVCF( vcf_file, maxNumberOfVariants = 100 ,
min_maf = 0.12 ,max_maf = NA)
# For realistic data use the function downloadGeneticMap
generateUniformGeneticMap()
startSimulation(vcf, totalNumberOfIndividuals = 200)
ped1 = newNuclearFamily(1)
v = computePairIBD2(1, 3)
cat("IBD2 of pair = ", v, "\n");
|
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.