multiecdf            package:geneplotter            R Documentation

_M_u_l_t_i_p_l_e _e_m_p_i_r_i_c_a_l _c_u_m_u_l_a_t_i_v_e _d_i_s_t_r_i_b_u_t_i_o_n _f_u_n_c_t_i_o_n_s (_e_c_d_f)
_a_n_d _d_e_n_s_i_t_i_e_s

_D_e_s_c_r_i_p_t_i_o_n:

     Plot multiple empirical cumulative distribution functions (ecdf)
     and densities with user interface similar to that of 'boxplot'.

_U_s_a_g_e:

     multiecdf(x, ...)
     ## S3 method for class 'formula':
     multiecdf(formula, data = NULL, ..., na.action = NULL)
     ## Default S3 method:
     multiecdf(x, xlim, col, do.points=FALSE, subsample=TRUE, ...)

     multidensity(x, ...)
     ## S3 method for class 'formula':
     multidensity(formula, data = NULL, main, xlab, ..., na.action = NULL)
     ## Default S3 method:
     multidensity(x, bw="nrd0", xlim, ylim, col, main, xlab, ...)

_A_r_g_u_m_e_n_t_s:

 formula: a formula, such as 'y ~ grp', where 'y' is a numeric vector
          of data values to be split into groups according to the
          grouping variable 'grp' (usually a factor).

    data: a data.frame (or list) from which the variables in 'formula'
          should be taken.

na.action: a function which indicates what should happen when the data
          contain 'NA's.  The default is to ignore missing values in
          either the response or the group.

       x: a list of numeric vectors.

      bw: the smoothing bandwidth, see the manual page for 'density'.
          If 'bw' is a character string specifying a rule to choose the
          bandwidth, this rule is applied to 'x[[1]]' and then the same
          numerical value of 'bw' is used throughout.

    xlim: Range of the x axis. If missing, the data range is used.

    ylim: Range of the y axis. If missing, the range of the density
          estimates is used.

     col: Line colors. If missing, the 'Set1' palette from RColorBrewer
          is used.

    main: Plot title.

    xlab: x-axis label.

do.points: logical; if 'TRUE', also draw points at the knot locations.

subsample: logical; if 'TRUE', subsamples of size 1000 are used to
          compute and plot the ecdf for list items with many
          observations ('>1000')

     ...: Further arguments that get passed on to the 'plot' functions.

_D_e_t_a_i_l_s:

     The usefulness of 'multidensity' can vary, depending on the data
     and because of smoothing artifacts. 'multiecdf' will in many cases
     be preferable.

_V_a_l_u_e:

     For the 'multidensity' functions, a list of 'density' objects.

_A_u_t_h_o_r(_s):

     Wolfgang Huber <URL: http://www.ebi.ac.uk/huber>

_S_e_e _A_l_s_o:

     'boxplot', 'ecdf' 'density'

_E_x_a_m_p_l_e_s:

       f = 1 + (runif(1000)>0.5)
       x = rnorm(length(f), mean=f, sd=f)
       
       multiecdf(x~f)
       multidensity(x~f)

