reconstructPPregion: Reconstruct perfect phylogeny sequencce across a region

Description Usage Arguments Value Examples

View source: R/reconstructPPregion.R

Description

This function reconstructs perfect phylogenies on each possible focal SNV across a genomic region.

Usage

1
reconstructPPregion(hapMat, minWindow, posn.lb = NULL, posn.ub = NULL)

Arguments

hapMat

A data structure of class hapMat. See the arguments in reconstructPP.

minWindow

Minimum number of SNVs around the focal SNV in the window of SNVs used to reconstruct the perfect phylogeny.

posn.lb

Lower bound of the subregion of hapMat (in base pairs) within which to consider SNVs.

posn.ub

Upper bound of the subregion of hapMat (in base pairs) within which to consider SNVs.

Value

An object of class multiPhylo that contains multiple phylo objects.

Examples

1
2
3
4
5
6
7
8
9
data(ex_hapMatSmall_data)   
                   
# Reconstruct partitions across the region of ex_hapMatSmall_data.
rdends <- reconstructPPregion(hapMat = ex_hapMatSmall_data,
                      minWindow = 1)
                      
# Reconstruct partitions between a given range SNV positions.
rdends_range <- reconstructPPregion(hapMat = ex_hapMatSmall_data, minWindow = 1,
                                      posn.lb = 2000, posn.ub = 7000)                      

perfectphyloR documentation built on March 8, 2021, 9:06 a.m.