sim_seq: Simulate sequence data for a pedigree

Description Usage Arguments Value

View source: R/sim_seqFunctions.R

Description

Simulate sequence data for a pedigree

Usage

1
2
sim_seq(ped_file, founder_genos, SNV_map, chrom_map, RV_marker,
  burn_in = 1000, gamma_params = c(2.63, 2.63/0.5))

Arguments

ped_file

Data frame. A single pedigree. Must match format of pedigree simulated by sim_RVped

founder_genos

Dataframe. A dataframe with rows corresponding to founders, and columns corresponding to markers. Markers must be listed in same order as SNV_map.

SNV_map

Data frame. A data frame that catalogs the SNVs in haplos. If the read_slim function was used to import SLiM data to R, the data frame Mutations is of the proper format for SNV_map. However, users must add the variable is_CRV to this data frame, see details.

chrom_map

Data.frame. A data.frame consisting of three columns: column 1 contains the chromosome numbers, column 2 start postion of chromosome (in cM), column 3 end position of chromosome (in cM).

RV_marker

character. The marker name of the RV locus.

burn_in

Numeric. The "burn-in" distance in centiMorgan, as defined by Voorrips and Maliepaard (2012), which is required before simulating the location of the first chiasmata with interference. By default, burn_in = 1000. The burn in distance in cM. By default, burn_in = 1000.

gamma_params

Numeric list of length 2. The respective shape and rate parameters of the gamma distribution used to simulate distance between chiasmata. By default, gamma _params = c(2.63, 2*2.63), as discussed in Voorrips and Maliepaard (2012).

Value

offspring_sequences


SimRVSequences documentation built on July 1, 2020, 6:03 p.m.