Description Usage Arguments Details Value Examples
Calcuate QC metrics from gene count matrix
1 | calculate_QC_metrics(sce)
|
sce |
a |
get QC metrics using gene count matrix. The QC statistics added are
number_of_genes number of genes detected.
total_count_per_cell sum of read number after UMI deduplication.
non_mt_percent 1 - percentage of mitochondrial gene counts. Mitochondrial genes are retrived by GO term GO:0005739
non_ERCC_percent ratio of exon counts to ERCC counts
non_ribo_percent 1 - percentage of ribosomal gene counts ribosomal genes are retrived by GO term GO:0005840.
an SingleCellExperiment
with updated QC metrics
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | data("sc_sample_data")
data("sc_sample_qc")
sce = SingleCellExperiment(assays = list(counts = as.matrix(sc_sample_data)))
organism(sce) = "mmusculus_gene_ensembl"
gene_id_type(sce) = "ensembl_gene_id"
QC_metrics(sce) = sc_sample_qc
demultiplex_info(sce) = cell_barcode_matching
UMI_dup_info(sce) = UMI_duplication
# The sample qc data already run through function `calculate_QC_metrics`.
# So we delete these columns and run `calculate_QC_metrics` to get them again:
colnames(colnames(QC_metrics(sce)))
QC_metrics(sce) = QC_metrics(sce)[,c("unaligned","aligned_unmapped","mapped_to_exon")]
sce = calculate_QC_metrics(sce)
colnames(QC_metrics(sce))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.