| group | R Documentation |
Identifies linkage groups of markers, using results from two-point (pairwise) analysis and the transitive property of linkage.
group(input.seq, LOD = NULL, max.rf = NULL, verbose = TRUE)
input.seq |
an object of class |
LOD |
a (positive) real number used as minimum LOD score (threshold) to declare linkage. |
max.rf |
a real number (usually smaller than 0.5) used as maximum recombination fraction to declare linkage. |
verbose |
logical. If |
If the arguments specifying thresholds used to group markers, i.e., minimum
LOD Score and maximum recombination fraction, are NULL (default),
the values used are those contained in object input.seq. If not
using NULL, the new values override the ones in object
input.seq.
Returns an object of class group, which is a list
containing the following components:
data.name |
name of
the object of class |
twopt |
name of the object of class |
marnames |
marker names, according to the input file. |
n.mar |
total number of markers. |
LOD |
minimum LOD Score to declare linkage. |
max.rf |
maximum recombination fraction to declare linkage. |
n.groups |
number of linkage groups found. |
groups |
number of the linkage group to which each marker is assigned. |
Gabriel R A Margarido, gramarga@gmail.com and Marcelo Mollinari, mmollina@usp.br
Lincoln, S. E., Daly, M. J. and Lander, E. S. (1993) Constructing genetic linkage maps with MAPMAKER/EXP Version 3.0: a tutorial and reference manual. A Whitehead Institute for Biomedical Research Technical Report.
rf_2pts and
make_seq
data(onemap_example_out)
twopts <- rf_2pts(onemap_example_out)
all.data <- make_seq(twopts,"all")
link_gr <- group(all.data)
link_gr
print(link_gr, details=FALSE) #omit the names of the markers
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.