LDHScpp

Expectation Maximisation (EM) algorithm

Paternal recombination rate and maternal linkage disequilibrium
(LD) are estimated for pairs of biallelic markers such as single nucleotide
polymorphisms (SNPs) from progeny genotypes and sire haplotypes. The
implementation relies on paternal half-sib families. If maternal half-sib
families are used, the roles of sire/dam are swapped. Multiple families can
be considered. For parameter estimation, at least one sire has to be double
heterozygous at the investigated pairs of SNPs.
Based on recombination rates, genetic distances between markers can be
estimated. Markers with unusually large recombination rate to markers in
close proximity (i.e. putatively misplaced markers) shall be discarded in
this derivation.
A workflow description is attached as vignette.
*A pipeline is available at GitHub*
<https://github.com/wittenburg/hsrecombi>
Hampel, Teuscher, Gomez-Raya, Doschoris, Wittenburg (2018) "Estimation of
recombination rate and maternal linkage disequilibrium in half-sibs"
<doi:10.3389/fgene.2018.00186>.
Gomez-Raya (2012) "Maximum likelihood estimation of linkage disequilibrium
in half-sib families" <doi:10.1534/genetics.111.137521>.

D<c3><b6>rte Wittenburg

hsrecombi

Estimation of Recombination Rate and Maternal LD in Half-Sibs

Dörte Wittenburg

LDHScpp function

<dl><dt>XGF1</dt>
<dd>integer matrix of progeny genotypes in genomic family 1</dd>
<dt>XGF2</dt>
<dd>integer matrix of progeny genotypes in genomic family 2</dd>
<dt>fAA</dt>
<dd>frequency of maternal haplotype 1-1</dd>
<dt>fAB</dt>
<dd>frequency of maternal haplotype 1-0</dd>
<dt>fBA</dt>
<dd>frequency of maternal haplotype 0-1</dd>
<dt>theta</dt>
<dd>paternal recombination rate</dd>
<dt>display</dt>
<dd>logical for displaying additional information</dd>
<dt>threshold</dt>
<dd>convergence criterion</dd></dl>

Arguments

Expectation Maximisation (EM) algorithm — LDHScpp

<dl>

<dt>XGF1</dt>
<dd>integer matrix of progeny genotypes in genomic family 1</dd>


<dt>XGF2</dt>
<dd>integer matrix of progeny genotypes in genomic family 2</dd>


<dt>fAA</dt>
<dd>frequency of maternal haplotype 1-1</dd>


<dt>fAB</dt>
<dd>frequency of maternal haplotype 1-0</dd>


<dt>fBA</dt>
<dd>frequency of maternal haplotype 0-1</dd>


<dt>theta</dt>
<dd>paternal recombination rate</dd>


<dt>display</dt>
<dd>logical for displaying additional information</dd>


<dt>threshold</dt>
<dd>convergence criterion</dd>

</dl>

LDHScpp: Expectation Maximisation (EM) algorithm

Description

Usage

Value

Arguments