sim_meiosis: Simulate meiosis

View source: R/sim_meiosis.R

sim_meiosisR Documentation

Simulate meiosis

Description

Output a random meiotic product from an input individual.

Usage

sim_meiosis(parent, m = 10, p = 0, obligate_chiasma = FALSE, Lstar = NULL)

Arguments

parent

An individual object, as output by create_parent() or cross()

m

interference parameter for chi-square model

p

Proportion of chiasmata coming from no-interference process.

obligate_chiasma

If TRUE, require an obligate chiasma on the 4-strand bundle at meiosis.

Lstar

Adjusted chromosome length, if obligate_chiasma=TRUE. Calculated if not provided.

Details

Simulations are under the Stahl model with the interference parameter being an integer. This is an extension of the chi-square model, but with chiasmata being the superposition of two processes, one following the chi-square model and the other exhibiting no interference.

Value

A list with alleles in chromosome intervals (as integers) and locations of the right endpoints of those intervals.

References

Copenhaver, G. P., Housworth, E. A. and Stahl, F. W. (2002) Crossover interference in arabidopsis. Genetics 160, 1631–1639.

Foss, E., Lande, R., Stahl, F. W. and Steinberg, C. M. (1993) Chiasma interference as a function of genetic distance. Genetics 133, 681–691.

Zhao, H., Speed, T. P. and McPeek, M. S. (1995) Statistical analysis of crossover interference using the chi-square model. Genetics 139, 1045–1056.

See Also

create_parent(), cross(), sim_crossovers(), calc_Lstar()

Examples

ind <- create_parent(100, 1:2)
prod <- sim_meiosis(ind)

kbroman/simcross documentation built on Jan. 13, 2024, 10:31 p.m.