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, xlab, na.action = NULL, ...)
     ## S3 method for class 'matrix':
     multiecdf(x, xlab, ...) 
     ## S3 method for class 'list':
     multiecdf(x,
               xlim,
               col = brewer.pal(9, "Set1"),
               main = "ecdf",
               xlab,
               do.points = FALSE,
               subsample = 1000L,
               legend = list(
                 x = "right",
                 legend = if(is.null(names(x))) paste(seq(along=x)) else names(x),
                 fill = col),
               ...)

     multidensity(x, ...)
     ## S3 method for class 'formula':
     multidensity(formula, data = NULL, xlab, na.action = NULL, ...)
     ## S3 method for class 'matrix':
     multidensity(x, xlab, ...) 
     ## S3 method for class 'list':
     multidensity(x,
                  bw = "nrd0",
                  xlim,
                  ylim,
                  col  = brewer.pal(9, "Set1"),
                  main = if(length(x)==1) "density" else "densities",
                  xlab,
                  lty  = 1L,
                  legend = list(
                    x = "topright",
                    legend = if(is.null(names(x))) paste(seq(along=x)) else names(x),
                    fill = col),
                  ...)

_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, lty: Line colors and line type.

    main: Plot title.

    xlab: x-axis label.

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

subsample: numeric or logical of length 1. If numeric, and larger than
          0, subsamples of that size are used to compute and plot the
          ecdf for those elements of 'x' with more than that number of
          observations. If logical and 'TRUE', a value of 1000 is used
          for the subsample size.

  legend: a list of arguments that is passed to the function 'legend'.

     ...: 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:

       words = strsplit(packageDescription("geneplotter")$Description, " ")[[1]]
       factr = factor(sample(words, 2000, replace = TRUE))
       x = rnorm(length(factr), mean=as.integer(factr))
       
       multiecdf(x ~ factr)
       multidensity(x ~ factr)

