onemap (version 2.1.1)

ripple_seq: Compares and displays plausible alternative orders for a given linkage group

Description

For a given sequence of ordered markers, computes the multipoint likelihood of alternative orders, by shuffling subsets (windows) of markers within the sequence. For each position of the window, all possible \((ws)!\) orders are compared.

Usage

ripple_seq(input.seq, ws = 4, ext.w = NULL, LOD = 3, tol = 0.1)

Arguments

input.seq

an object of class sequence with a predefined order.

ws

an integer specifying the length of the window size (defaults to 4).

ext.w

an integer specifying how many markers should be considered in the vicinity of the permuted window. If ext.w=NULL all markers in the sequence are considered. In this veriosn, it is used only in backcross, \(F_2\) or RIL crosses.

LOD

threshold for the LOD-Score, so that alternative orders with LOD less then or equal to this threshold will be displayed.

tol

tolerance for the C routine, i.e., the value used to evaluate convergence.

Value

This function does not return any value; it just produces text output to suggest alternative orders.

Details

Large values for the window size make computations very slow, specially if there are many partially informative markers.

References

Broman, K. W., Wu, H., Churchill, G., Sen, S., Yandell, B. (2008) qtl: Tools for analyzing QTL experiments R package version 1.09-43

Jiang, C. and Zeng, Z.-B. (1997). Mapping quantitative trait loci with dominant and missing markers in various crosses from two inbred lines. Genetica 101: 47-58.

Lander, E. S., Green, P., Abrahamson, J., Barlow, A., Daly, M. J., Lincoln, S. E. and Newburg, L. (1987) MAPMAKER: An interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics 1: 174-181.

Mollinari, M., Margarido, G. R. A., Vencovsky, R. and Garcia, A. A. F. (2009) Evaluation of algorithms used to order markers on genetics maps. Heredity 103: 494-502.

Wu, R., Ma, C.-X., Painter, I. and Zeng, Z.-B. (2002a) Simultaneous maximum likelihood estimation of linkage and linkage phases in outcrossing species. Theoretical Population Biology 61: 349-363.

Wu, R., Ma, C.-X., Wu, S. S. and Zeng, Z.-B. (2002b). Linkage mapping of sex-specific differences. Genetical Research 79: 85-96

See Also

make_seq, compare, try_seq and order_seq.

Examples

# NOT RUN {
# }
# NOT RUN {
 #Outcross example
  data(example_out)
  twopt <- rf_2pts(example_out)
  markers <- make_seq(twopt,c(27,16,20,4,19,21,23,9,24,29))
  markers.map <- map(markers)
  ripple_seq(markers.map)

#F2 example
 data(onemap_example_f2)
 twopt <- rf_2pts(onemap_example_f2)
 all_mark <- make_seq(twopt,"all")
 groups <- group(all_mark)
 LG3 <- make_seq(groups,3)
 LG3.ord <- order_seq(LG3, subset.search = "twopt", twopt.alg = "rcd", touchdown=TRUE)
 LG3.ord
 make_seq(LG3.ord) # get safe sequence
 ord.1<-make_seq(LG3.ord,"force") # get forced sequence
 ripple_seq(ord.1, ws=5)
# }
# NOT RUN {

# }