compareorder: Compare two orderings of markers on a chromosome

In qtl: Tools for Analyzing QTL Experiments

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)

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.

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.

Author(s)

Karl W Broman, [e-mail]

See Also

ripple, switch.order, movemarker

Examples

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

qtl documentation built on Sept. 11, 2024, 5:43 p.m.