loglikfun

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

loglikfun function

<dl><dt>counts</dt>
<dd>integer vector of observed 2-locus genotype</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>fBB</dt>
<dd>frequency of maternal haplotype 0-0</dd>
<dt>theta</dt>
<dd>paternal recombination rate</dd></dl>

Arguments

Calculate log-likelihood function — loglikfun

<dl>

<dt>counts</dt>
<dd>integer vector of observed 2-locus genotype</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>fBB</dt>
<dd>frequency of maternal haplotype 0-0</dd>


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

</dl>

loglikfun: Calculate log-likelihood function

Description

Value

Arguments