CNE-class: Class '"CNE"'

CNE-classR Documentation

Class "CNE"

Description

CNE class contains all the meta-data of CNEs, including the pair of assemblies, the thresholds, the intermediate and final CNE sets.

Usage

### Constructors:
CNE(assembly1Fn=character(1), assembly2Fn=character(1),
    axt12Fn=character(), axt21Fn=character(),
    window=50L, identity=50L,
    CNE12=GRangePairs(), CNE21=GRangePairs(),
    CNEMerged=GRangePairs(), CNEFinal=GRangePairs(),
    aligner="blat", cutoffs1=4L, cutoffs2=4L)

### Accessor-like methods:
## S4 method for signature 'CNE'
thresholds(x)
## S4 method for signature 'CNE'
CNE12(x)
## S4 method for signature 'CNE'
CNE21(x)
## S4 method for signature 'CNE'
CNEMerged(x)
## S4 method for signature 'CNE'
CNEFinal(x)

## ... and more (see Methods)

Arguments

assembly1Fn,assembly2Fn

Object of class "character": The twoBit filenames of assembly1, assembly2

axt12Fn,axt21Fn

Object of class "character": The Axt filenames of assembly1 to assembly2, assembly2 to assembly1

window

Object of class "integer": The window size for scanning CNEs. By default, it is 50.

identity

Object of class "integer": The identity over the window size for scanning CNEs. By default, it is 50.

CNE12

Object of class "GRangePairs": The preliminary CNEs from axt file with assembly1 as reference.

CNE21

Object of class "GRangePairs": The preliminary CNEs from axt file with assembly2 as reference.

CNEMerged

Object of class "GRangePairs": The CNEs after merging CNE1 and CNE2.

CNEFinal

Object of class "GRangePairs": The CNEs after being realigned back to reference genome, with blat in current implementation.

aligner

Object of class "character": The method to realign CNEs back to the reference genome.

cutoffs1, cutoffs2

Object of class "integer": The CNEs with more than the cutoff hits on the reference genome are removed.

x

Object of class "CNE": A "CNE" object.

Methods

CNE12

signature(x = "CNE"): Get the CNE1 results.

CNE21

signature(x = "CNE"): Get the CNE2 results.

CNEMerged

signature(x = "CNE"): Get the merged CNE results.

CNEFinal

signature(x = "CNE"): Get the final CNE results.

thresholds

signature(x = "CNE"): Get the thresholds used for scanning CNEs.

Author(s)

Ge Tan

Examples

  library(GenomicRanges)
  ## Constructor
  CNE12 <- GRangePairs(first=GRanges(seqnames=c("chr13", "chr4", "chr4"),
                                     ranges=IRanges(start=c(71727138,150679343,
                                                            146653164),
                                                    end=c(71727224, 150679400,
                                                          146653221)),
                                     strand="+"),
                       second=GRanges(seqnames=c("chr1"),
                                      ranges=IRanges(start=c(29854162, 23432387,
                                                             35711077),
                                                     end=c(29854248, 23432444,
                                                           35711134)),
                                      strand="+")
                       )
  CNE21 <- GRangePairs(first=GRanges(seqnames=c("chr1"),
                                    ranges=IRanges(start=c(29854162, 23432387,
                                                           35711077),
                                                   end=c(29854248, 23432444,
                                                         35711134)),
                                    strand="+"),
                       second=GRanges(seqnames=c("chr13", "chr4", "chr4"),
                                      ranges=IRanges(start=c(71727138,150679343,
                                                             146653164),
                                                     end=c(71727224, 150679400,
                                                           146653221)),
                                      strand="+")
                       )
  cne <- CNE(assembly1Fn=file.path(system.file("extdata",
                                    package="BSgenome.Drerio.UCSC.danRer10"),
                        "single_sequences.2bit"),
             assembly2Fn=file.path(system.file("extdata",
                                    package="BSgenome.Hsapiens.UCSC.hg38"),
                        "single_sequences.2bit"),
             window=50L, identity=50L,
             CNE12=CNE12, CNE21=CNE21, CNEMerged=CNE12, CNEFinal=CNE12,
             aligner="blat", cutoffs1=4L, cutoffs2=4L)
             
  ## Accessor
  CNE12(cne)
  CNE21(cne)
  thresholds(cne)
  CNEMerged(cne)
  CNEFinal(cne)

ge11232002/CNEr documentation built on Oct. 26, 2022, 7:08 p.m.