View source: R/5_get_deg_all.R
| get_deg_all | R Documentation |
Do differential analysis according to an expression set and group information.
get_deg_all(
exp,
group_list,
ids = NULL,
symmetry = TRUE,
my_genes = NULL,
show_rownames = FALSE,
cluster_cols = TRUE,
color_volcano = c("#2874C5", "grey", "#f87669"),
logFC_cutoff = 1,
pvalue_cutoff = 0.05,
adjust = TRUE,
entriz = TRUE,
n_cutoff = 2,
annotation_legend = FALSE,
lab = NA,
change_col = "change",
species = "human"
)
exp |
A numeric matrix |
group_list |
A character or factor vector with one label per sample; character input will be converted to factor and unused levels will be dropped. |
ids |
optional probe-to-symbol annotation table. If 'NULL', 'exp' is treated as a gene-level matrix and no ID conversion is performed. |
symmetry |
a logical value ,would you like to get your plot symmetrical |
my_genes |
genes for pheatmap |
show_rownames |
logical,show rownames or not |
cluster_cols |
boolean values determining if columns should be clustered or |
color_volcano |
color for volcano |
logFC_cutoff |
Cutoff value of logFC,1 by default. |
pvalue_cutoff |
Cutoff value of pvalue,0.05 by default. |
adjust |
a logical value; if TRUE, use |
entriz |
logical, if TRUE, convert symbols to Entrez IDs. |
n_cutoff |
numeric cutoff used after scaling; values outside
|
annotation_legend |
logical,show annotation legend or not |
lab |
label for x axis in volcano plot; when 'NA', pkg = 4 uses 'log2FoldChange', and the other package codes keep their original labels |
change_col |
column name for the up/down/status label. Default uses 'change' when it exists; otherwise the function falls back to the existing automatic classification logic. |
species |
choose human, mouse, or rat; default is human |
a list with a deg data.frame, volcano plot, PCA plot, heatmap, and a list with DEGs.
Xiaojie Sun
get_deg;multi_deg_all
## Not run:
if(requireNamespace("Biobase",quietly = TRUE)&
requireNamespace("AnnoProbe",quietly = TRUE)){
gse = "GSE42872"
geo = geo_download(gse,destdir=tempdir())
group_list = rep(c("A","B"),each = 3)
group_list = factor(group_list)
find_anno(geo$gpl)
ids <- AnnoProbe::idmap(geo$gpl,destdir = tempdir())
dcp = get_deg_all(geo$exp,group_list,ids,entriz = FALSE)
head(dcp$deg)
dcp$plots
}else{
if(!requireNamespace("AnnoProbe",quietly = TRUE)) {
warning("Package 'AnnoProbe' needed for this function to work.
Please install it by install.packages('AnnoProbe')",call. = FALSE)
}
if(!requireNamespace("Biobase",quietly = TRUE)) {
warning("Package 'Biobase' needed for this function to work.
Please install it by BiocManager::install('Biobase')",call. = FALSE)
}
}
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.