Summarizes the results of the median comparison function
calcMedians
and the pvalue calculation function
calcPvalues
. Should be used with the matrices originating from
these two functions.
1  calcSummary(medians, pvalues)

medians 
a matrix containing median difference values calculated by
the 
pvalues 
a matrix containing pvalues calculated by the

All genes with a median comparison value >= 0.05 and a pvalue of <= 0.01 are summarized into a data.frame. These genes are assumed to contain a batch effect
A data frame with the columns "gene" containing the gene name, "batch" containing the batch number from which the gene was found, "median" and "pvalue" containing the calculated median difference values and the pvalues, respectively.
calcMedians
calcPvalues
correctBatchEffect
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16  ## Shortly running example. For a more realistic example that takes
## some more time, run the same procedure with the full BEclearData
## dataset.
## Whole procedure that has to be done to use this function.
data(BEclearData)
ex.data < ex.data[31:90,7:26]
ex.samples < ex.samples[7:26,]
# Calculates median difference values and pvalues from the example data
med < calcMedians(data=ex.data, samples=ex.samples, parallel=FALSE)
pvals < calcPvalues(data=ex.data, samples=ex.samples, parallel=FALSE,
adjusted=TRUE, method="fdr")
# Summarize pvalues and median differences for batch affected genes
sum < calcSummary(medians=med, pvalues=pvals)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.