knitr::read_chunk("../../analysis/chunks.R")
library(ggplot2) library(ggrastr) library(ggpubr) library(dplyr) library(Matrix) library(parallel) library(dropestr) library(dropEstAnalysis) theme_set(theme_base) set.seed(42) kOutputFolder <- '../../output/' kDataPath <- '../../data/' kEstDataPath <- paste0(kDataPath, 'dropest/SRR1784310/runtimes/est_2018_03_15_real/')
holder <- readRDS(paste0(kEstDataPath, 'cell.counts.rds'))
times_coll <- sapply(1:10, function(i) TimeCollisions(holder)) mean(times_coll, trim=0.2)
times_directional <- sapply(1:10, function(i) TimeUmiErrors(holder$reads_per_umi_per_cell, 'directional', mc.cores=10)) print(times_directional) mean(times_directional, trim=0.1)
times_bayesian <- sapply(1:5, function(i) TimeUmiErrors(holder$reads_per_umi_per_cell, mc.cores=10)) print(times_bayesian) mean(times_bayesian, trim=0.1)
times_quality <- sapply(1:10, function(i) TimeQuality(holder, mit.chromosome.name='chrM')) print(times_quality) mean(times_quality, trim=0.1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.