group {onemap} | 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) ## S3 method for class 'group' print(x, detailed=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. |
x |
an object of class |
detailed |
logical. If |
... |
further arguments, passed to other methods. Currently ignored. |
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 overridden 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
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.
data(example.out) twopts <- rf.2pts(example.out) all.data <- make.seq(twopts,"all") link_gr <- group(all.data) link_gr