globalSeg {GeneR}R Documentation

Class manipulating segments

Description

We made a set of tools manipulation segments

Details

The aim of this class is to describe regions on chromosomes that are discontinuous segments on a line like:

            
            1        10             25             40
Region A:   ##########              #######
Region B:   ####  ####              #############
Region C:                                          ######  ####

We made two kind of class

For a better comprehension of other man pages, we introduce this notation:

Author(s)

Epissage group at CGM.

See Also

See also as.segSet, as.globalSeg, unionSeg

Examples


a = list(
    matrix( c( 1, 15, 17,  5, 45, 38),ncol=2),
    matrix( c( 100 , 120),ncol=2),
    matrix( c( 130, 135, 140, 145),ncol=2),
    matrix( c( 142 , 160),ncol=2))


b = list(
    matrix( c(15, 28, 18, 45),ncol=2),
    matrix( c(1, 15, 25, 10, 20, 40),ncol=2),
    matrix( c(17, 35, 23, 38),ncol=2),
    matrix( c(100, 110, 105, 120),ncol=2))

a = as.globalSeg(a)
b = as.globalSeg(b)

par(mar=c(1,0,1,0))

par(mfrow=c(8,1))
plot(a,xlim=c(1,160),main="A")
plot(b,xlim=c(1,160),main="B")

plot(and(b),xlim=c(1,160),main="and(B)")
plot(or(b),xlim=c(1,160),main="or(B)")
plot(Xor(b),xlim=c(1,160),main="Xor(B)")

plot(a&b,xlim=c(1,160),main="A&B")
plot(a|b,xlim=c(1,160),main="A|B")
plot(Xor(a,b),xlim=c(1,160),main="Xor(A,B)")



[Package GeneR version 2.24.0 Index]