Downsampling

zUMIs has powerful downsampling capabilities. Independent of downsampling mode, the full data is always exported as well.

• Adaptive downsampling: According to the recommendation of the Scater package (McCarthy et al., 2017) reads are downsampled to be within 3 times median absolute deviation. This is the default setting.
• downsampling to a fixed depth: Reads are downsampled to a user-specified depth. Any barcodes that do not reach the requested depth are omitted. Example: 10000
• downsampling to a depth range: Barcodes with read depth above the maximum of the range are downsampled to this value. All barcodes within the range are reported without downsampling and barcodes below the minimum specified read depth are ommited. Example: 10000-20000
• downsampling to several depths: Several depths can be requested by comma separation. Combinations of fixed depth and depth ranges may be given. Example: 10000,10000-20000,30000

Output

Downsampled count tables are reported in .dgecounts.rds for each feature type (exons, introns, intron+exon). You will find a list object of all the downsamplings requested in each counting type "readcount" & "umicount".

For an example how to access these tables, see the code below:

AllCounts <- readRDS("zUMIs_output/expression/example.dgecounts.rds")

downsamp <- AllCounts$umicount$exon$downsampling$downsampled_10000