Learn R Programming

qtl (version 1.39-5)

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)

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.

See Also

ripple, switch.order, movemarker

Examples

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

Run the code above in your browser using DataLab