Downsample read counts

Share:

Description

This function downsamples the library sizes of a read counts table to the lowest library size, according to the methdology used in (Soneson and Delorenzi, BMC Bioinformatics, 2013).

Usage

1
    downsample.counts(counts,seed=42)

Arguments

counts

the read counts table which is subjected to downsampling.

seed

random seed for reproducible downsampling.

Value

The downsampled counts matrix.

Author(s)

Panagiotis Moulos

Examples

1
2
3
4
5
6
7
## Dowload locally the file "bottomly_count_table.txt" from
## the ReCount database
#download.file(paste("http://bowtie-bio.sourceforge.net/",
#    "recount/countTables/bottomly_count_table.txt",sep=""),
#    destfile="~/bottomly_count_table.txt")
#M <- as.matrix(read.delim("~/bottomly_count_table.txt",row.names=1))
#D <- downsample.counts(M)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.