MixThres-package          package:MixThres          R Documentation

_M_i_x_t_u_r_e _m_o_d_e_l_s _t_o _d_e_f_i_n_e _a _h_y_b_r_i_d_i_z_a_t_i_o_n _T_h_r_e_s_h_o_l_d

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

     A hybridization threshold is defined based on the histogramm of
     the probe intensity values. The procedure is divided in two steps.
     First the intensity distribution is estimated using a particular
     mixture model. Second a hybridization threshold is defined from
     the components of the mixture.

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


       Package:  MixThres
       Type:     Package
       Version:  1.0
       Date:     2008-09-23
       License:  gpl-2

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

     J. Aubert <julie.aubert@agroparistech.fr>

_R_e_f_e_r_e_n_c_e_s:

     Picard et al., Using mixture models to define a hybridization
     threshold in DNA microarray experiments.

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

     'hybridization','threshold'

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

     ## Data simulation
     Signal <- c(rnorm(100,5,1),rnorm(700,8,1),rnorm(200,12,1))
     Data <- data.frame(Gene=paste("Gene",1:1000,sep=""),Signal=Signal)

     ## Mixture models
     out <- hybridization(Data)

     ## Threshold corresponding to the selected mixture model and a particular precision criterion
     out <- threshold(out,crit= 1e-04)

     ## Printing and graphics
     print(out) 
     hist(out)
     plot(out)

