paralell_pairwise_discrete: Parallel Pairwise Discrete Estimation

View source: R/pairwise_rf.R

paralell_pairwise_discreteR Documentation

Parallel Pairwise Discrete Estimation

Description

This function performs parallel pairwise estimation of recombination fractions using discrete dosage scoring via a C++ backend.

Usage

paralell_pairwise_discrete(
  mrk.pairs,
  input.seq,
  geno,
  dP,
  dQ,
  count.cache,
  tol = .Machine$double.eps^0.25,
  ll = ll
)

Arguments

mrk.pairs

A matrix of dimensions 2*N, containing N pairs of markers to be analyzed.

input.seq

An object of class mappoly.sequence.

geno

Genotype matrix.

dP

Vector of probabilities for the first allele.

dQ

Vector of probabilities for the second allele.

count.cache

An object of class cache.info containing pre-computed genotype frequencies.

tol

The tolerance level for the estimation accuracy (default is .Machine$double.eps^0.25).

ll

Logical; if TRUE, the function returns log-likelihood values instead of LOD scores. For internal use.

Value

Depending on the ll parameter, returns either log-likelihood values or formatted LOD scores from pairwise recombination fraction estimation.


mappoly documentation built on May 29, 2024, 6:05 a.m.