I have worked on a lot of Gene Chips and used several methods to define genes that make up the background (genes that are not expressed):
- filter out genes where mean expression is below mean expression of eset (normalized expression set)
- filter out genes where expression is below mean expression of eset in 50% or 80% of samples,
- filter out genes where expression is below mean expression of eset in 50% or 80% of either cases or controls etc.
- Use the mean expression of the antigenomic probes (which are considered to be background) as a cutoff for main probes (real genes)
But I don't know which is the best way to go about detecting genes that are expressed above background (or to define genes that make up the background, for that matter). I know there is a function
paCalls() in the oligo package for exon level expression but are there are any functions that could do the same thing for gene level expression? Some function that would give you a p-value at gene-level as opposed to at exon-level as in the
paCalls() function. Any suggestions would be appreciated.