Chapter 7 Introduction
Circular visualization is popular in Genomics and related omics fields. It is efficient in revealing associations in high dimensional genomic data. In genomic plots, categories are usually chromosomes and data on x axes are genomic positions, but it can also be any kind of general genomic categories.
To make is easy for Genomics analysis, circlize package particularly provides functions which focus on genomic plots. These functions are synonymous to the basic graphic functions but expect special format of input data:
circos.genomicTrack(): create a new track and add graphics.
circos.genomicPoints(): low-level function, add points.
circos.genomicLines(): low-level function, add lines or segments.
circos.genomicRect(): low-level function, add rectangles.
circos.genomicText(): low-level function, add text.
circos.genomicLink(): add links.
The genomic functions are implemented by basic circlize functions (e.g.
circos.points()), thus, the use of genomic functions can be mixed with the basic circlize functions.
7.1 Input data
Genomic data is usually stored as a table where the first three columns define the genomic regions and following columns are values associated with the corresponding regions. Each genomic region is composed by three elements: genomic category (in most case, it is the chromosome), start position on the genomic category and the end position. Such data structure is known as BED format and is broadly used in genomic research.
circlize provides a simple function
generateRandomBed() which generates random genomic data. Positions are uniformly generated from human genome and the number of regions on chromosomes approximately proportional to the length of chromosomes. In the function,
nc control the number of rows and numeric columns that users need. Please note
nr are not exactly the same as the number of rows which are returned by the function.
fun argument is a self-defined function to generate random values.
set.seed(999) bed = generateRandomBed() head(bed)
## chr start end value1 ## 1 chr1 39485 159163 -0.1887635 ## 2 chr1 897195 1041959 -0.2435220 ## 3 chr1 1161957 1177159 0.3749953 ## 4 chr1 1201513 1481406 -0.2600839 ## 5 chr1 1487402 1531773 -0.4633990 ## 6 chr1 1769949 2736215 -0.8159909
bed = generateRandomBed(nr = 200, nc = 4) nrow(bed)
##  205
bed = generateRandomBed(nc = 2, fun = function(k) sample(letters, k, replace = TRUE)) head(bed)
## chr start end value1 value2 ## 1 chr1 98740 566688 e e ## 2 chr1 769960 887938 b q ## 3 chr1 906851 933021 u o ## 4 chr1 1241911 1243537 k f ## 5 chr1 1385344 1410947 v x ## 6 chr1 1498302 1585389 u v
All genomic functions in circlize expect input variable as a data frame which contains genomic data or a list of data frames which contains genomic data in different conditions.