PolyHaplotyper (version 1.0.1)

compareHapMrkDosages: compare haplotyping results with observed markers dosages

Description

compare haplotyping results with observed markers dosages

Usage

compareHapMrkDosages(mrkDosage, hapresults)

Arguments

mrkDosage

a data.frame or matrix with the marker dosages, in the format of inferHaplotypes()

hapresults

a list as returned by inferHaplotypes, with one item (itself a list) per haploblock with at least a matrix hapdos and a character vector markers. All markers in all haploblocks must also occur in mrkDosages, and all individuals in the hapdos matrices must also occur in mrkDosages

Value

a 3-D array with dimensions haploblock, individual, and chkresult: mrkNA (ALL markers in the haploblock have missing data for the individual), hapNA (the haplotype dosages do not sum to ploidy and/or are missing), match (TRUE if the non-missing marker dosages match the haplotype dosages, FALSE is there is a conflict, NA if mrkNA and/or hapNA are TRUE) Each element is itself a list with elements:

Examples

# NOT RUN {
data(PolyHaplotyper_small)
chmd <- compareHapMrkDosages(mrkDosage=phdos, hapresults=phresults)
# show results for first haploblock, first 8 individuals:
chmd[1, 1:8,]
# }