compareorder: Compare two orderings of markers on a chromosome

Description

Compare the likelihood of an alternative order for markers on a chromosome to the current order.

Usage

compareorder(cross, chr, order, error.prob=0.0001,
             map.function=c("haldane","kosambi","c-f","morgan"),
             maxit=4000, tol=1e-6, sex.sp=TRUE)

Value

A data frame with two rows: the current order in the input cross object, and the revised order. The first column is the log10 likelihood of the new order relative to the original one (positive values indicate that the new order is better supported). The second column is the estimated genetic length of the chromosome for each order. In the case of sex-specific maps, there are separate columns for the female and male genetic lengths.

Arguments

cross: An object of class cross. See read.cross for details.
chr: The chromosome to investigate. Only one chromosome is allowed. (This should be a character string referring to the chromosomes by name.)
order: The alternate order of markers on the chromosome: a numeric vector that is a permutation of the integers from 1 to the number of markers on the chromosome.
error.prob: Assumed genotyping error rate used in the calculation of the penetrance Pr(observed genotype | true genotype).
map.function: Indicates whether to use the Haldane, Kosambi, Carter-Falconer, or Morgan map function when converting genetic distances into recombination fractions.
maxit: Maximum number of EM iterations to perform.
tol: Tolerance for determining convergence.
sex.sp: Indicates whether to estimate sex-specific maps; this is used only for the 4-way cross.

Author

Karl W Broman, broman@wisc.edu

Examples

Run this code

data(badorder)
compareorder(badorder, chr=1, order=c(1:8,11,10,9,12))