MENU
R'MES general view rmes.gaussien rmes.poisson.composee rmes.poisson
rmes.gfam rmes.format
rmes.histo
rmes.compar rmes.pyramide


Splus Function: rmes.histo
Histograms of Results from R'MES Programs.



Description Usage Required Arguments Optional Arguments Details Examples


Description

The Splus function rmes.histo draws histograms of the word statistics calculated with a R'MES program: rmes.gaussien, rmes.poisson, or rmes.poisson.composee.

Usage

 

rmes.histo (don, dessin, fact, nclass=20, lang="fr")

Required Arguments
don
Character string equal to the pathname of the file containing the results from R'MES. This file should have been created by one of the following commands: rmes.gaussien, rmes.poisson, rmes.poisson.composee.
There will be as many drawn histograms as word lengths considered in this file.

Sommaire

Optional Arguments

dessin
Character string equal to the pathname of the PostScript file where the graphics will be stored. If this file exists, it will be replaced.
When this option is not set, the current graphical device is used. If there is no current active graphical device, X11 is used.
fact
Integer vector to define the scale on the vertical axis on each histogram.
Each element corresponds then to a word length. It is equal to the exponent of the power of 2 by which the number of words (i.e 4^word length) must be divided to obtain the maximum on the y-axis (see example 3).
By default, the greatest exponent that is suitable for the current histogram is calculated.
This argument is specially relevant for comparing histograms of different data sets.
nclass
Number of classes for the histograms.
lang
When equal to "fr", labels and messages are in french. When equal to "en", labels and messages are in english. By default, they are in french.

Sommaire

Details

If lg.max (resp. lg.min) is the maximal (resp. minimal) length of the words analyzed in the file don, the statistics of the words of length lg.min to lg.max are represented by a histogram. There is one histogram per word length.

Sommaire

Examples

1- Histograms of the statistics of the words of length 7 to 8 in the lambda genome under the model M1 using a Gaussian approximation. The two plots are stored in a PostScript file. 

rmes.gaussien -o lambda.gaussien.7-8_1 -seq lambda 
              -hmin 7 -hmax 8 -m 1 


> rmes.histo("lambda.gaussien.7-8_1.0", dessin="histo1.ps", lang="en")
2- Same representation but 40 classes for each histogram. 

> rmes.histo("lambda.gaussien.7-8_1.0", dessin="histo2.ps", nclass=40)
3- Same representation but maxima are set on the y-axis. They are respectively: 
  8192 = 4^7 / 2^2 for the 7-words
 32768 = 4^8 / 2^1 for the 8-words.


> rmes.histo("lambda.gaussien.7-8_1.0", dessin="histo3.ps", 
              fact=c(2,1), nclass=40)

Sommaire