MixThres-package {MixThres} | R Documentation |
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.
Package: | MixThres |
Type: | Package |
Version: | 1.0 |
Date: | 2008-09-23 |
License: | gpl-2 |
J. Aubert <julie.aubert@agroparistech.fr>
Picard et al., Using mixture models to define a hybridization threshold in DNA microarray experiments.
## 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)