quality               package:PROcess               R Documentation

_Q_u_a_l_i_t_y _C_h_e_c_k _o_n _a _S_e_t _o_f _S_p_e_c_t_r_a

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

     Compute three quality parameters for a set of spectra.

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

     quality(Ma, peakinfofile, cutoff)

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

      Ma: a Matrix where the baseline-subtracted spectra are stored
          column wise.

peakinfofile: a `.csv' file in the same format as  Ciphergen's peak
          info file, with 5 columns data. See Details of 'getPeaks'.

  cutoff: The point in m/z below which spectra are cutoff.

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

     The quality parameters are computed a la fashion de Mani. 1.
     Estimate noise by moving average with a 5 point window. 2.
     Estimate the noise envelop by 3 times the standard deviation of
     noise in a 251 point moving window. 3. Compute the area under the
     baseline-subtracted curve, area0. 4. Compute the area after
     subtracting the noise envelop from the baseline-subtracted curve,
     area1. 5. Parameter 'Quality' is defined as area1/area0. 6.
     Parameter 'Retain' is defined as the number of points with height
     above 5 times the noise envelop over total number of points in the
     spectrum. 7. Detect peaks in each spectrum by 'getPeaks' or
     Ciphergen software. 8. Parameter 'peak' is defined as the number
     of peaks in a spectrum divided by the mean number of peaks across
     spectra.

     A spectrum is considered to be of poor quality if Quality<0.4,
     Retain<0.1 and peak<0.5 simultaneously.

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

     A matrix with three named columns, 'Quality','Retain' and 'peak',
     with spectrum file names as row names.

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

     Xiaochun Li

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

     example(getPeaks)
     qualRes <- quality(testM, peakfile, cutoff=1500)

