Next: , Previous: MapView, Up: Visualization


4.2 MMapView

MMapView is designed to draw markers and QTL for one linkage group for several genetic maps. It helps to visualize common markers between maps but also to display marker interval distance heterogeneities.

4.2.1 Command Line Options

Option Usage Type Explanation
-m,--mapdir required string The XML file/directory for the genetic map(s).
-c,--chrom required string The name of the chromosome to display.
-o,--output required string The output file stem.
-r,--refmap optional string The XML file for the reference map.
-q,--withqtl optional integer The threshold on the occurence of the markers.
-t,--tonto optional string The XML file for the trait ontology
--htest optional boolean Test distance homogeneity between the maps and the reference (only with -r).
--hth optional double The threshold of the p-value of the homogeneity test between the maps and the reference (only with -r and -h). The value must be between ]0,1[. A value of 0 or 1 codes for a gradient view.
--mrkt optional integer Threshold on the occurence of the markers.
-p,--parfile optional string The drawing parameter file.
--img optional string The format of the image (use –help for possible values).
Suppose we have a set of genetic maps in XML format in the directory xml and that all the maps have a chromosome called 1. Then we can use MMapView to display all the chromosomes in a single figure as follows:

     
     %java org.metaqtl.main.MMapView \
     >     -c 1 -m xml -o figure
This dumps the image into a file figure_1.jpeg and the parameter file into figure_1.par. We can also represent the same chromosomes but with only the markers which have observed at least in two distinct chromosomes by using the option --mrkt:

     
     %java org.metaqtl.main.MMapView \
     >     -c 1 -m xml -o figure --mrkt 2
Now, suppose we want to draw the same chromosomes but also to represent the eventual interval marker distance heterogenities relatively to a reference chromosome defined in file reference.xml. Then we use the command,

     
     %java org.metaqtl.main.MMapView \
     >     -c 1 -m xml -r reference.xml --htest -o figure

4.2.2 Drawing parameters

MMapView use the same drawing parameters than MapView plus some extra parameters. If you want to display common marker links between adjacent chromosomes set `WITH_COMMON_MARKER' to `true'. The width of the lines which connect the common markers can be modified by using the parameter `COMMON_STROKE_WIDTH'. This is illustrated in the next figure.

fig3.jpg

The link between common markers can also be painted depending on the way of the common sequences in which they are involved. For example, in the above figure the common markers which are involved in common sequence correctly ordered between the two chromosomes have their links painted in blue, otherwise their links are painted in red. Note that one marker is not involved in a common sequence and its link is painted in gray. These colors can be changed by using the parameters `SINGLE_COMMON_COLOR', `POS_COMMON_COLOR' and `NEG_COMMON_COLOR'. For example, to create the previous figure we added the following lines to the parameter file:

     
     #
     # DRAW COMMON MARKER LINKS
     #
     WITH_COMMON_MARKER=true
     COMMON_STROKE_WIDTH=2.0
     SINGLE_COMMON_COLOR=808080
     POS_COMMON_COLOR=0000ff
     NEG_COMMON_COLOR=ff0000
When MMapView is used to visualize the test of homogeneity between a reference chromosome and several other chromosomes, you can parametrize both the number of bin colors to represent the probabilities and the type of gradient. For example, to create the following figure

fig4.jpg

we added the extra parameters:

     
     #
     # PROBABILITY GRADIENT
     #
     PROBA_BIN=20
     PROBA_FROM_COLOR=ffffff
     PROBA_TO_COLOR=ff0000

This means that the interval between 0 and 1 is divided into 20 intervals and the gradient starts in white (ffffff) and ends in red (ff0000). Note that the probability gradient is represented into the legend box.

fig5.jpg

Finally the chromosomes displayed in the figure can be aligned either relatively to the first marker of each chromosome, CHROM_ALIGN_MODE=1, or relatively to the first common marker between the first chromosome (the left one) and the other ones, CHROM_ALIGN_MODE=0.