readCNData: Read Copy Number Data

In facopy: Feature-based association and gene-set enrichment for copy number alteration analysis in cancer

Description

The function reads the copy number output generated by one of the supported copy number analysis tools

Usage

readCNData(folder, method = NULL, sex = c("X"), FUN = NULL, version = NULL, 
           window = NULL, rankThr = NULL, pfbFilename, lengthThr = 10, ...)

Arguments

folder	A character indicating the path where the files with copy number information are located. In order to avoid possible conflicts, try not to place any additional files within the folder.
method	Either NULL or a character. If a character, it specifies the method used generate the copy number output. Currently supported: "seqcna", "cnanorm", "patchwork", "freec", "oncosnp", "oncosnp-seq", "gap" and "exomecnv". If NULL, a generic input will be assumed with the following columns: chromosome, segment start, segment end, segment length and segment type. The segment types should be coded as 1 (deletion), 2 (normal), 3 (amplification), 11 (LOH deletion), 12 (LOH normal) or 13 (LOH amplification).
sex	Either NULL or a character vector. If a character vector, it should contain the names of the sexual chromosomes that should be read.
FUN	A pattern replacement function. Some copy number analysis tools add suffixes and file extension characters to the names of their output files. This function allows keeping just the sample names. Leave to NULL for default replacement based on the selected method.
version	A character with the software version of the copy number analysis tool. The output may vary between versions of the same method. Leave to NULL if unknown and check your software version if reading the copy number data is not successful.
window	An integer. Required for output files from CNAnorm, which do not contain information on the window length used during the copy number analysis process.
rankThr	An integer. Threshold on the Rank parameter provided by OncoSNP-SEQ and the newer versions of OncoSNP (>1.3).
pfbFilename	For GAP, name of a file with population B allele frequencies (PFB), the same one used during the copy number analysis process.
lengthThr	An integer. For GAP, filter those alterations that span less SNPs than this threshold.
...	If a method is not selected (NULL), further arguments passed to the read.delim function, in charge of reading the input files from the selected folder.

Value

A facopyInfo object with copy number data.

Author(s)

David Mosen-Ansorena

Examples

# myCalls = readCNData("~/myFolder/", "seqcna")
# myCalls = readCNData("~/myFolder/", "gap", pfbFilename="~/myPfb.pfb")
# myCalls = readCNData("~/myFolder/", "cnanorm", window=50000)

facopy documentation built on May 2, 2018, 2:30 a.m.