read_htseq: Read multiple HTSeq-count data into a matrix

In jzsh2000/jzshr: jzsh2000's R Code

Description

HTSeq-count output a tab-delimed count file, while the first column records gene_id, and the second column records raw read count. read_htseq merges multiple HTSeq-count output into one matrix. The count files could be selected using Sys.glob.

Usage

read_htseq(files, labels = NULL, clean = TRUE)

Arguments

files	HTSeq-count output files
labels	sample labels used as colnames in count matrix
clean	clean meta-line (start with '__') or not

Value

a count matrix, with rows as genes and columns as samples

jzsh2000/jzshr documentation built on May 20, 2019, 6:29 a.m.