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")

_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.

       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.

_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).

_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) 

