auxillary_functions: Auxillary functions

In jebbd/sim_vcf: A simple package to simulate a vcf file of random variants in HWE

Description

get_pos: get a random postion on a chromosome

mut_ref: mutate a reference base semi randomly, with a Ti/Tv of roughly 2, needs to wrapped in a map/lapply for working with vectors

rev_encode: reverse the encoding of genotypes from numbers to that use in VCF i.e. 0/0, 0/1, 1/1

Usage

get_pos(chrom)

mut_ref(base)

rev_encode(x)

Arguments

chrom	A character vector of chromsomes formatted as e.g. "chr1"
base	an single DNA nucleotide (ATCG) to mutate
x	a vector of numbers from 0 - 2 to encode. ANy other numbers will become no calls ("./.")

jebbd/sim_vcf documentation built on Dec. 20, 2021, 10:08 p.m.