ppsetApply {affy} | R Documentation |
Apply a function over the ProbeSets in an AffyBatch
ppsetApply(abatch, FUN, genenames = NULL, ...) ppset.ttest(ppset, covariate, pmcorrect.fun = pmcorrect.pmonly, ...)
abatch |
An object inheriting from AffyBatch . |
ppset |
An object of class ProbeSet . |
covariate |
the name a covariate in the slot phenoData . |
pmcorrect.fun |
a function to corrrect PM intensities |
FUN |
A function working on a ProbeSet |
genenames |
A list of Affymetrix probesets ids to work with. All
probe set ids used when NULL . |
... |
Optional parameters to the function FUN |
Returns a list
of objects, or values, as returned by the
function FUN
for each ProbeSet
it processes.
Laurent Gautier <laurent@cbs.dtu.dk>
ppset.ttest <- function(ppset, covariate, pmcorrect.fun = pmcorrect.pmonly, ...) { probes <- do.call("pmcorrect.fun", list(ppset)) my.ttest <- function(x) { y <- split(x, get(covariate)) t.test(y[[1]], y[[2]])$p.value } r <- apply(probes, 1, my.ttest) return(r) } ##this takes a long time - and rowttests is a good alternative ## eg: rt = rowttests(exprs(Dilution), Dilution$liver) ## Not run: data(Dilution) all.ttest <- ppsetApply(Dilution, ppset.ttest, covariate="liver") ## End(Not run)