dyeswapfilter            package:maanova            R Documentation

_G_e_n_e _f_i_l_t_e_r _f_o_r _d_y_e-_s_w_a_p _e_x_p_e_r_i_m_e_n_t

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

     This function is used to flag the questionable spot in any kind of
     dye-swap experiment.

     This function only works for 2-dye arrays.

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

     dyeswapfilter(dataobj, r=4)

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

 dataobj: An object of class 'madata' or 'rawdata'.

       r: A cut-off value for bad spot. The genes with log-ratio
          difference larger than r times standard deviation will be
          flagged.

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

     For each pair of dye-swap, the difference in log ratios (d) are
     computed. Then compute the IQR (interquartile range) of d and
     convert that to Standard Deviation by SD = IQR/1.35. Any gene with
     d larger than r times SD will be flagged.

     Note that I assume in the input data object, the adjecent arrays
     is a dye-swap pair.

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

     An object of class 'rawdata' or 'madata' with the 'flag' field
     created or updated.

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

     Hao Wu hao@jax.org

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

     ## Not run: 
     data(kidney)
     # riplot before filtering
     riplot(kidney.raw, array=1)
     # filter the gene
     rawdata <- dyeswapfilter(kidney.raw)
     # riplot again - some genes are highlighted
     riplot(rawdata, array=1)
     ## End(Not run)

