hapiAutoPhase: Automatic inference of haplotypes

In Jialab-UCR/Hapi: Inference of Chromosome-length Haplotypes using Genomic Data of Single Gamete Cells

Description

Automatic inference of haplotypes

Usage

hapiAutoPhase(gmt, code = "atcg")

Arguments

gmt	a dataframe of genotype data of gamete cells
code	a character indicating the code style of genotype data. One of 'atcg' and '01'. Default is 'atcg'

Value

a dataframe of inferred consensus haplotypes

Author(s)

Ruidong Li

Examples

ref <- sample(c('A','T'),500, replace=TRUE)
alt <- sample(c('C','G'),500, replace=TRUE)

gmt <- data.frame(chr=rep(1,500), pos=seq_len(500),
    ref=ref, alt=alt, gmt1=ref, gmt2=alt, gmt3=ref,
    gmt4=ref, gmt5=c(alt[1:250], ref[251:500]),
    stringsAsFactors = FALSE)
    
hapOutput <- hapiAutoPhase(gmt=gmt, code='atcg')

Jialab-UCR/Hapi documentation built on May 30, 2019, 11:41 a.m.