BBB_NMP_ad_subset

A dataset containing real allele depth information for 1,000 SNPs. The data
is available for ten individuals and the first 1,000 SNPs on chromosome 1.
The data corresponds to a low-fold sequencing experiment. There are two
columns per individuals (read counts for allele 1 and for allele 2).

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_ad_subset function

A data frame with 1,000 rows and 25 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_ad1</dt>
<dd>The read count for the first individual at the first markrer</dd> <dt>id1_ad2</dt>
<dd>The read count for the first individual at the second markrer</dd> <dt>id2_ad1</dt>
<dd>The read count for the second individual at the first markrer</dd> <dt>id2_ad2</dt>
<dd>The read count for the second individual at the second markrer</dd> <dt>id3_ad1, id3_ad2, id4_ad1, id4_ad2, id5_ad1, id5_ad2, id6_ad1, id6_ad2, id7_ad1, id7_ad2, id8_ad1, id8_ad2, id9_ad1, id9_ad2, id10_ad1, id10_ad2</dt>
<dd>The read counts for the other individuals</dd>
 
</dl>

Format

Example for "ad" format specification — BBB_NMP_ad_subset

A data frame with 1,000 rows and 25 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_ad1</dt>
<dd>The read count for the first individual at the first markrer</dd>

 <dt>id1_ad2</dt>
<dd>The read count for the first individual at the second markrer</dd>

 <dt>id2_ad1</dt>
<dd>The read count for the second individual at the first markrer</dd>

 <dt>id2_ad2</dt>
<dd>The read count for the second individual at the second markrer</dd>

 <dt>id3_ad1, id3_ad2, id4_ad1, id4_ad2, id5_ad1, id5_ad2, id6_ad1, id6_ad2, id7_ad1, id7_ad2, id8_ad1, id8_ad2, id9_ad1, id9_ad2, id10_ad1, id10_ad2</dt>
<dd>The read counts for the other individuals</dd>
 
</dl>

BBB_NMP_ad_subset: Example for "ad" format specification

Description

Usage

Arguments

Format