inveRsion (version 1.20.0)

listInv-methods: Constructor of class inversionList

Description

Determines the full inversion sequence from overlapping trial segments of fixed window sized. It computes the limits of the inversion, the population frequency and retrieves the maximum BIC across the trial segments.

Usage

listInv(object, hapCode, geno, ROI, saveRes, thBic, all, saveROI)

Arguments

object

scan. Results of scanning the genome for an inversion with a fixed window size; output of scanInv.

hapCode

HaploCode. Object with the result of coding haplotypes for each candidate brake-point; output of codeHaplo

geno

logical. Whether original data is genotypes of phased haplotypes.

ROI

numeric. 2-vector specification passes the chromosome segment to be analyzed. 4-vector specification passes the region of interest for the left brake-point (ROI[1] and ROI[2]) and the right brake-point (ROI[3] and ROI[4]). ROI should be specified in mega-basis units.

saveRes

logical.Whether results should be saved into file list Inv.RData

thBic

numeric. BIC threshold above which trial segments are selected.

all

logical. Whether recomputing within the ROI should be done for all possible segment sizes or the window size in scan.

saveROI

logical. Whether saving the blocks for the candidate break-points for all the ROIs.

Value

inversionList

Methods

signature(object = "scan")

Details

listInv is both a method for class scan and constructor of inversionList. It re-runs the inversion model within the ROIs found in the previous scan. However, it is also possible to explicitly pass the ROIs defined by the user. The re-run is done with the same window size of the scan, which is convenient if enough significant trial segments where found within the inversion segment. If a more detailed re-run is needed set all=TRUE. This computes the model for all possible trial segments of any length within the ROI. Note that for high SNP density this can be computational intensive.

See Also

HaploCode,scan,inversionList

Examples

data(scanRes)
data(hapCode)
invList<-listInv(scanRes,hapCode=hapCode,geno=FALSE,saveRes=FALSE,all=FALSE,thBic=0, saveROI=FALSE)
invList