| calcSummary {BEclear} | R Documentation |
Summarizes the results of the median comparison function
calcMedians and the p-value calculation function
calcPvalues. Should be used with the matrices originating from
these two functions.
calcSummary(medians, pvalues)
medians |
a matrix containing median difference values calculated by
the |
pvalues |
a matrix containing p-values calculated by the
|
calcSummary
All genes with a median comparison value >= 0.05 and a p-value of <= 0.01 are summarized into a data.frame. These genes are assumed to contain a batch effect
Null if there are no batch effects detected, else A data frame with the columns "gene" containing the gene name, "batch" containing the batch number from which the gene was found, "median" and "p-value" containing the calculated median difference values and the p-values, respectively.
## 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 p-values from the example data
library(data.table)
samples <- data.table(ex.samples)
data <- data.table(feature=rownames(ex.data), ex.data)
data <- melt(data = data, id.vars = "feature", variable.name = "sample",
value.name = "beta.value")
setkey(data, "feature", "sample")
med <- calcMedians(data=data, samples=samples)
pvals <- calcPvalues(data=data, samples=samples, adjusted=TRUE,
method="fdr")
# Summarize p-values and median differences for batch affected genes
sum <- calcSummary(medians=med, pvalues=pvals)