edge (version 2.4.2)

readDGE: Read and Merge a Set of Files Containing Count Data

Description

Reads and merges a set of text files containing gene expression counts.

Usage

readDGE(files, path=NULL, columns=c(1,2), group=NULL, labels=NULL, ...)

Arguments

files

character vector of filenames, or a data.frame of sample information containing a column called files.

path

character string giving the directory containing the files. Defaults to the current working directory.

columns

numeric vector stating which columns of the input files contain the gene names and counts respectively.

group

optional vector or factor indicating the experimental group to which each file belongs.

labels

character vector giving short names to associate with the files. Defaults to the file names.

...

other arguments are passed to read.delim.

Value

• A DGEList object containing a matrix of counts, with a row for each unique tag found in the input files and a column for each input file.

Details

Each file is assumed to contain digital gene expression data for one genomic sample or count library, with gene identifiers in the first column and counts in the second column. Gene identifiers are assumed to be unique and not repeated in any one file. The function creates a combined table of counts with rows for genes and columns for samples. A count of zero will be entered for any gene that was not found in any particular sample. By default, the files are assumed to be tab-delimited and to contain column headings. Other file formats can be handled by adding arguments to be passed to read.delim. For example, use header=FALSE if there are no column headings and use sep="," to read a comma-separated file. Instead of being a vector, the argument files can be a data.frame containing all the necessary sample information. In that case, the filenames and group identifiers can be given as columns files and group respectively, and the labels can be given as the row.names of the data.frame.

See Also

See read.delim for other possible arguments that can be accepted. DGEList-class, DGEList.

Examples

#  Read all .txt files from current working directory

files <- dir(pattern="*\\.txt$")
RG <- readDGE(files)