getlength: Retrieves Gene length data
In goseq: Gene Ontology analyser for RNA-seq and other length biased data

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/getlength.R

Gets the length of each gene in a vector.

1	getlength(genes, genome, id)

`genes`	A vector or list of the genes for which length information is required.
`genome`	A string identifying the genome that `genes` refer to. For a list of supported organisms run `supportedGenomes`.
`id`	A string identifying the gene identifier used by `genes`. For a list of supported gene IDs run `supportedGeneIDs`.

Length data is obtained from data obtained from the UCSC genome browser for each combination of genome and id. As fetching this data at runtime is time consuming, a local copy of the length information for common genomes and gene ID are included in the geneLenDataBase package. This function uses this package to fetch the required data.

The length of a gene is taken to be the median length of all its mature, mRNA, transcripts. It is always preferable to obtain length information directly for the gene ID used to summarize your count data, rather than converting IDs and then using the supplied databases. Even when two genes have a one-to-one mapping between different identifier conventions (which is often not the case), they frequently refer to slightly different regions of the genome with different lengths. It is therefore recommended that the user perform the full analysis in terms of only one gene ID, or manually obtain their own length data for the identifier used to bin reads by gene.

Returns a vector of the gene lengths, in the same order as genes. If length data is unavailable for a particular gene NA is returned in that position. The returned vector is intended for use with the bias.data option of the nullp function.

Matthew D. Young myoung@wehi.edu.au

supportedGenomes, supportedGeneIDs, nullp, geneLenDataBase