Description Usage Arguments Details Value Examples
This function performs OVESEG-test to assess significance of molecular markers.
1 2 3 | OVESEGtest(y, group, weights = NULL, alpha = "moderated",
NumPerm = 999, seed = 111, detail.return = TRUE,
BPPARAM = bpparam())
|
y |
a numeric matrix containing log-expression or logCPM (log2-counts per million) values. Data frame or SummarizedExperiment object will be internally coerced into a matrix. Rows correspond to probes and columns to samples. Missing values are not permitted. |
group |
categorical vector or factor giving group membership of columns of y. At least two categories need to be presented. |
weights |
optional numeric matrix containing prior weights for each spot. |
alpha |
parameter specifying within-group variance estimator to be used.
'moderated': empirical Bayes moderated variance estimator as used in
|
NumPerm |
an integer specifying the number of permutation resamplings (default 999). |
seed |
an integer seed for the random number generator. |
detail.return |
a logical indicating whether more details about posterior probability estimation will be returned. |
BPPARAM |
a BiocParallelParam object indicating whether parallelization should be used for permutation resamplings. The default is bpparam(). |
OVESEG-test is a statistically-principled method that can detect tissue/cell-specific marker genes among many subtypes. OVESEG-test statistics are designed to mathematically match the definition of molecular markers, and a novel permutation scheme are employed to estimate the corresponding distribution under null hypotheses where the expression patterns of non-markers can be highly complex.
a list containing the following components:
pv.overall |
a vector of p-values calculated by all permutations regardless of upregulated subtypes. |
pv.oneside |
a vector of subtype-specific p-values calculated specifically for the upregulated subtype of each probe. |
pv.oneside.max |
subtype-specific p-values when observed test statistic equal to zero. |
pv.multiside |
pv.oneside*K (K-time comparison correction) and truncated at 1. |
W |
a matrix of posterior probabilities for each component null hypothesis given an observed probe. Rows correspond to probes and columns to one hypothesis. |
label |
a vector of group labels. |
groupOrder |
a matrix with each row being group indexes ordered based
on decreasing expression levels.
Group indexes are positions in |
F.p.value |
a matrix with each column giving p-values corresponding to F-statistics on certain groups. |
lfdr |
a matrix with each column being local false discovery rates estimated based on one column of weighted F.p.value matrix. |
fit |
a |
F.p.value
, lfdr
and fit
are returned only when
detail.return
is TRUE.
1 2 3 4 5 6 7 8 9 | data(RocheBT)
rtest <- OVESEGtest(RocheBT$y, RocheBT$group, NumPerm=99,
BPPARAM=BiocParallel::SerialParam())
## Not run:
#parallel computing
rtest <- OVESEGtest(RocheBT$y, RocheBT$group, NumPerm=99,
BPPARAM=BiocParallel::SnowParam())
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.