The UCSC data base provides a wealth of annotation information. This
function can be used to access UCSC, to retrieve the data available
there and to return it as an annotation track object ameanable to
plotting with plotTracks
.
clearSessionCache
is can be called to remove all cached items
from the session which are generated when connecting with the UCSC
data base.
UcscTrack(track, table=NULL, trackType=c("AnnotationTrack",
"GeneRegionTrack", "DataTrack", "GenomeAxisTrack"), genome, chromosome,
name=NULL, from, to, ...)
clearSessionCache()
NULL
, in which case the default selection of tables is
used. To find out about available tables for a given track please
consult the online table browser at
http://genome.ucsc.edu/cgi-bin/hgTables?command=start.c("AnnotationTrack",
"GeneRegionTrack", "DataTrack", "GenomeAxisTrack")
. The function
will try to coerce the downloaded data in an object of this
class. See below for details.trackType
argument. See Details section for more
information.trackType
. The data stored at the UCSC data bases can be of different formats:
gene or transcript model data, simple annotation features like CpG
Island locations or SNPs, or numeric data like conservation or
mapability. This function presents a unified API to download all kinds
of data and to map them back to one of the annotation track objects
defined in this package. The type of object to hold the data has to be
given in the trackType
argument, and subsequently the function
passes all data on to the respective object constructor. All
additional named arguments are considered to be relevant for the
constructor of choice, and single character scalars are replaced by
the respective data columns in the dowloaded UCSC tables if
available. For instance, assuming the table for track 'foo' contains
the columns 'id', 'type', 'fromLoc' and 'toLoc', giving the featuer
identifier, type, start end end location. In order to create an
AnnotationTrack
object from that data, we have to
pass the additional named arguments id="id"
,
feature="type"
, start="fromLoc"
and codeend="toLoc" to the
UcscTrack
function. The complete function call could look like
this:
UcscTrack(track="foo", genome="mm9", chromosome=3,
from=1000, to=10000, trackType="AnnotationTrack", id="id",
feature="type", start="from", end="to")
To reduce the bandwidth, some caching of the UCSC connection takes
place. In order to remove these cached session items, call
clearSessionCache
.
The Gviz.ucscUrl
option controls which URL is being used to
connect to UCSC. For instance, one could switch to the European UCSC
mirror by calling
options(Gviz.ucscUrl="http://genome-euro.ucsc.edu/cgi-bin/"
.
AnnotationTrack
DataTrack
GeneRegionTrack
GenomeAxisTrack