fitNorm2                package:prada                R Documentation

_F_i_t _b_i_v_a_r_i_a_t_e _n_o_r_m_a_l _d_i_s_t_r_i_b_u_t_i_o_n.

_D_e_s_c_r_i_p_t_i_o_n:

     Fits a bivariate normal distribution into a data set of paired
     values and selects data points according to their standard
     deviation from the fitted distribution.

_U_s_a_g_e:

     fitNorm2(x, y=NA, scalefac=1, method="covMcd", noise, gateName = "fitNorm")

_A_r_g_u_m_e_n_t_s:

       x: Numeric vector containing x-value or n by 2 matrix containing
          x and y values or object of class 'cytoFrame'.

       y: Numeric vector containing y-value (optional). The length of
          'x' must be the same as that of 'y'.

scalefac: Numeric vector giving factor of standard deviations used for
          data selection (all points within 'scalefac' standard
          deviations are selected).

  method: One of 'covMcd' or 'cov.rob' defining method used for
          computation of covariance matrix.

   noise: Numeric or logical index vector defining value pairs in x
          that are not used for fitting of distributions. Can be used
          to deal with noisy data.

gateName: Character giving the name of the gate object.

_D_e_t_a_i_l_s:

     Computes the densities of a bivariate normal distribution from the
     covariance matrix of the paired data. Covariance matrices are
     acquired either by function 'covMcd' (considerably faster) or by
     function 'cov.rob'.

_V_a_l_u_e:

     A list containing items 'mu' (midpoint of distribution), 'S'
     (covariance matrix), 'p' (density values for each data pair),
     'sel' (selection of data points), 'scalefac' (factor of standard
     deviations used for data selection), 'data' (x and y values of
     data points) and 'gate', an object of class 'gate' containing the
     selection.

_A_u_t_h_o_r(_s):

     Florian Hahne

_S_e_e _A_l_s_o:

     'cov.rob', 'covMcd', 'plotNorm2'

_E_x_a_m_p_l_e_s:

      sampdat <- readFCS(system.file("extdata",
        "fas-Bcl2-plate323-04-04.A01", package="prada"))
      nfit    <- fitNorm2(exprs(sampdat[,1:2]), scalefac=2)
      plotNorm2(nfit, selection=TRUE, ellipse=TRUE) 

