pairwise.filter: Filter pairwise comparison statistics between two...

In simpleaffy: Very simple high level analysis of Affymetrix data

Description

Given the results of a pairwise.comparison, filter the resulting gene list on expression level, PMA calls (if available), fold change and t-test statistic.

min.exp and min.exp.no allow the data to be filtered on intensity (where min.exp.no specifies the minimum number of arrays that must be above the threshold 'min.exp' to be allowed through the filter).

PMA filtering is done when min.present.number is greater than 0.

min.present.no allows arrays to be filtered by PMA call. A number or 'all' must be specified. If a number, then the at least this many arrays must be called present, if 'all', then all arrays must be called present.

present.by.group specifies whether PMA filtering is to be done on a per-group basis or for all arrays at once. If 'false' then the experiment is treated as a single group (i.e. a probeset passes the filter if it is called present on at least min.present.number arrays in the whole experiment. If 'true' then it must be called present on at least this many arrays in one or more groups. See the vignette for more details.

Usage

  pairwise.filter(object,min.exp=log2(100),min.exp.no=0,min.present.no=0,present.by.group=T,fc=1.0,tt=0.001)

Arguments

object	a 'PairComp' object
min.exp	Filter genes using a minimum expression cut off
min.exp.no	A gene must have an expression intensity greater than 'min.exp' in at least this number of chips
min.present.no	A gene must be called present on at least this number of chips
present.by.group	If true, then the probeset must be called Present on at least min.present.number arrays in any of the replicate sets used to generate the PairComp object being filtered. If false, then must be called present on at least min.present.no of the arrays in the whole experiment
fc	A gene must show a log2 fold change greater than this to be called significant
tt	A gene must be changing with a p-score less than this to be called significant

Value

A 'PairComp' object filtered to contain only the genes that pass the specified filter parameters.

Author(s)

Crispin J Miller

References

http://bioinformatics.picr.man.ac.uk/

Examples

 
  ## Not run: 
    pc <- pairwise.comparison(eset.rma,"group",c("A","P"))
    pf <- pairwise.filter(pc,tt=0.01);
  
## End(Not run)

simpleaffy documentation built on Nov. 8, 2020, 11:09 p.m.