BBB_NMP_GP_subset

A dataset containing real genotype probabilities for 1,000 SNPs. The data is
available for ten individuals and 1,000 SNPs on chromosome 10. The data
corresponds to a low-fold sequencing experiment. There are three columns per
individuals (for genotypes 00, 01 and 11).

datasets

Functions to identify Homozygous-by-Descent (HBD) segments associated with runs of homozygosity (ROH) and to
estimate individual autozygosity (or inbreeding coefficient). HBD segments and autozygosity are assigned to multiple HBD classes
with a model-based approach relying on a mixture of exponential distributions. The rate of the exponential distribution is distinct
for each HBD class and defines the expected length of the HBD segments. These HBD classes are therefore related to the age of the
segments (longer segments and smaller rates for recent autozygosity / recent common ancestor). The functions allow to estimate the
parameters of the model (rates of the exponential distributions, mixing proportions), to estimate global and local autozygosity
probabilities and to identify HBD segments with the Viterbi decoding. The method is fully described in Druet and Gautier (2017)
<doi:10.1111/mec.14324>.

Tom Druet

RZooRoH

Partitioning of Individual Autozygosity into Multiple
Homozygous-by-Descent Classes

BBB_NMP_GP_subset function

A data frame with 1,000 rows and 35 variables:<dl>
 <dt>chr</dt>
<dd>The chromosome number</dd> <dt>marker_name</dt>
<dd>The marker id</dd> <dt>pos</dt>
<dd>The position of the marker</dd> <dt>allele1</dt>
<dd>The name of the first marker allele</dd> <dt>allele2</dt>
<dd>The name of the second marker allele</dd> <dt>id1_gp1</dt>
<dd>The AA genotype probability for the first individual</dd> <dt>id1_gp2</dt>
<dd>The AB genotype probability for the first individual</dd> <dt>id1_gp3</dt>
<dd>The BB genotype probability for the first individual</dd> <dt>id2_gp1</dt>
<dd>The AA genotype probability for the second individual</dd> <dt>id2_gp2</dt>
<dd>The AB genotype probability for the second individual</dd> <dt>id2_gp3</dt>
<dd>The BB genotype probability for the second individual</dd> <dt>id3_gp1, id3_gp2, id3_gp3, id4_gp1, id4_gp2, id4_gp3, id5_gp1, id5_gp2, id5_gp3, id6_gp1, id6_gp2, id6_gp3, id7_gp1, id7_gp2, id7_gp3, id8_gp1, id8_gp2, id8_gp3, id9_gp1, id9_gp2, id9_gp3, id10_gp1, id10_gp2, id10_gp3</dt>
<dd>The genotype probabilities for the other individuals</dd>
 
</dl>

Format

Example for "gp" format specification — BBB_NMP_GP_subset

A data frame with 1,000 rows and 35 variables:<dl>
 <dt>chr</dt>
<dd>The chromosome number</dd>

 <dt>marker_name</dt>
<dd>The marker id</dd>

 <dt>pos</dt>
<dd>The position of the marker</dd>

 <dt>allele1</dt>
<dd>The name of the first marker allele</dd>

 <dt>allele2</dt>
<dd>The name of the second marker allele</dd>

 <dt>id1_gp1</dt>
<dd>The AA genotype probability for the first individual</dd>

 <dt>id1_gp2</dt>
<dd>The AB genotype probability for the first individual</dd>

 <dt>id1_gp3</dt>
<dd>The BB genotype probability for the first individual</dd>

 <dt>id2_gp1</dt>
<dd>The AA genotype probability for the second individual</dd>

 <dt>id2_gp2</dt>
<dd>The AB genotype probability for the second individual</dd>

 <dt>id2_gp3</dt>
<dd>The BB genotype probability for the second individual</dd>

 <dt>id3_gp1, id3_gp2, id3_gp3, id4_gp1, id4_gp2, id4_gp3, id5_gp1, id5_gp2, id5_gp3, id6_gp1, id6_gp2, id6_gp3, id7_gp1, id7_gp2, id7_gp3, id8_gp1, id8_gp2, id8_gp3, id9_gp1, id9_gp2, id9_gp3, id10_gp1, id10_gp2, id10_gp3</dt>
<dd>The genotype probabilities for the other individuals</dd>
 
</dl>

BBB_NMP_GP_subset: Example for "gp" format specification

Description

Usage

Arguments

Format