MixThres-package {MixThres}R Documentation

Mixture models to define a hybridization Threshold

Description

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.

Details

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

Author(s)

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

References

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

See Also

hybridization,threshold

Examples

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

[Package MixThres version 1.0 Index]