Extracts regions of possible inversion events from "scan"

Description

lists the regions of interest in a matrix, given by the overlapping of significant trial segments (of fixed window seize) that may be part of an inversion event.

Usage

getROIs(object, thBic)

Arguments

Details

ROIs are defined as overlapping trial segments with BIC greater than thBIC. The output is a matrix for which each row is one ROI. The first two columns give intervals defining the left brake-points and the two subsequent columns are the intervals for the right brake-points. ROIs are given in mega-basis.

Value

matrix

Methods

signature(object = "scan")

list of regions of interest

Author(s)

Alejandro Caceres acaceres@creal.cat

See Also

scan

Examples

data(scanRes)
ROI <- getROIs(scanRes, thBic = 0)
ROI

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(inveRsion)
Loading required package: haplo.stats


Hola!

welcome to inevRsion package. 
 
 

type: manual() for full manual 
      vignette("inveRsion") for a quick start 

> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/inveRsion/getROIs-methods.Rd_%03d_medium.png", width=480, height=480)
> ### Name: getROIs-methods
> ### Title: Extracts regions of possible inversion events from "scan"
> ### Aliases: getROIs-methods getROIs,scan-method getROIs
> ### Keywords: methods
> 
> ### ** Examples
> 
> data(scanRes)
> ROI <- getROIs(scanRes, thBic = 0)
  1 ROI extracted  

> ROI
  LeftBinf LeftBsup RightBinf RightBsup
1 0.742854 0.757891  1.243445  1.258031
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>