listInv-methods: Constructor of class inversionList

In inveRsion: Inversions in genotype data

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.

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.

Value

inversionList

Methods

signature(object = "scan")

Author(s)

Alejandro Caceres [e-mail]

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

inveRsion documentation built on Nov. 8, 2020, 5:08 p.m.