Learn R Programming

hsphase (version 3.0.0)

phf: Half-sib Family Phasing

Description

Phases a half-sib family using the block structure and an imputed sire haplotype matrix.

Usage

phf(GenotypeMatrix, blockMatrix, sirePhasedMatrix)

Value

Returns a matrix containing the phased parental haplotypes of the half-sibs (two rows per individual). Alleles are coded as 0 (A), 1 (B), and 9 (missing/unphased).

Arguments

GenotypeMatrix

matrix. Half-sib genotypes (one half-sib per row; SNPs ordered by mapping position in the columns). Data should be numeric: 0, 1, 2 for AA, AB, BB. Use 9 for missing data.

blockMatrix

matrix. Blocking structure (output of bmh).

sirePhasedMatrix

matrix. Imputed sire haplotypes (output of ssp).

See Also

aio

Examples

Run this code
genotype <- matrix(c(
  2,1,0,
  2,0,0,
  0,0,2
), byrow = TRUE, ncol = 3)

block <- bmh(genotype)
phf(genotype, block, ssp(block, genotype))

Run the code above in your browser using DataLab