hapiIdentifyCV: Indentify crossovers in gamete cells

Description Usage Arguments Value Author(s) Examples

View source: R/crossover.R

Description

Indentify crossovers in gamete cells

Usage

1
hapiIdentifyCV(hap, gmt, hmm = NULL)

Arguments

hap

a dataframe of the two haplotypes

gmt

a dataframe of genotype data of gamete cells

hmm

a list containing probabilities of a HMM. Default is NULL

Value

a dataframe containing crossover information in each gamete cell

Author(s)

Ruidong Li

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
ref <- sample(c('A','T'),500, replace=TRUE)
alt <- sample(c('C','G'),500, replace=TRUE)

hap <- data.frame(hap1=ref, hap2=alt, stringsAsFactors = FALSE)
rownames(hap) <- seq_len(500)

gmt <- data.frame(gmt1=ref, gmt2=alt, gmt3=ref,
    gmt4=ref, gmt5=c(alt[1:250], ref[251:500]),
    stringsAsFactors = FALSE)
    
cvOutput <- hapiIdentifyCV(hap=hap, gmt=gmt)

Hapi documentation built on May 2, 2019, 6:52 a.m.