Calculates distribution of track expressions' values over the given set of bins.
gdist(
...,
intervals = NULL,
include.lowest = FALSE,
iterator = NULL,
band = NULL
)pairs of 'expr', 'breaks' where 'expr' is a track expression and the breaks determine the bin
genomic scope for which the function is applied
if 'TRUE', the lowest value of the range determined by breaks is included
track expression iterator. If 'NULL' iterator is determined implicitly based on track expressions.
track expression band. If 'NULL' no band is used.
N-dimensional vector where N is the number of 'expr'-'breaks' pairs.
This function calculates the distribution of values of the numeric track expressions over the given set of bins.
The range of bins is determined by 'breaks' argument. For example: 'breaks=c(x1, x2, x3, x4)' represents three different intervals (bins): (x1, x2], (x2, x3], (x3, x4].
If 'include.lowest' is 'TRUE' the the lowest value will be included in the first interval, i.e. in [x1, x2]
'gdist' can work with any number of dimensions. If more than one 'expr'-'breaks' pair is passed, the result is a multidimensional vector, and an individual value can be accessed by [i1,i2,...,iN] notation, where 'i1' is the first track and 'iN' is the last track expression.
# \dontshow{
options(gmax.processes = 2)
# }
gdb.init_examples()
## calculate the distribution of dense_track for bins:
## (0, 0.2], (0.2, 0.5] and (0.5, 1]
gdist("dense_track", c(0, 0.2, 0.5, 1))
#> (0,0.2] (0.2,0.5] (0.5,1]
#> 12463 1518 10
#> attr(,"breaks")
#> attr(,"breaks")[[1]]
#> [1] 0.0 0.2 0.5 1.0
#>
## calculate two-dimensional distribution:
## dense_track vs. sparse_track
gdist("dense_track", seq(0, 1, by = 0.1), "sparse_track",
seq(0, 2, by = 0.2),
iterator = 100
)
#> (0,0.2] (0.2,0.4] (0.4,0.6] (0.6,0.8] (0.8,1] (1,1.2] (1.2,1.4]
#> (0,0.1] 0 58 114 8 3 0 0
#> (0.1,0.2] 0 377 828 38 0 0 0
#> (0.2,0.3] 0 5 294 68 21 0 0
#> (0.3,0.4] 0 2 13 20 4 2 1
#> (0.4,0.5] 0 0 0 0 1 0 0
#> (0.5,0.6] 0 0 0 1 0 0 0
#> (0.6,0.7] 0 0 0 0 0 0 0
#> (0.7,0.8] 0 0 0 0 0 0 0
#> (0.8,0.9] 0 0 0 0 0 0 0
#> (0.9,1] 0 0 0 0 0 0 0
#> (1.4,1.6] (1.6,1.8] (1.8,2]
#> (0,0.1] 0 0 0
#> (0.1,0.2] 0 0 0
#> (0.2,0.3] 0 0 0
#> (0.3,0.4] 0 0 0
#> (0.4,0.5] 0 0 0
#> (0.5,0.6] 0 0 0
#> (0.6,0.7] 0 0 0
#> (0.7,0.8] 0 0 0
#> (0.8,0.9] 0 0 0
#> (0.9,1] 0 0 0
#> attr(,"breaks")
#> attr(,"breaks")[[1]]
#> [1] 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
#>
#> attr(,"breaks")[[2]]
#> [1] 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
#>