simulate.diallel: Simulate diallel data from the BayesDiallel model

Description Usage Arguments Examples

View source: R/diallel_sim.R

Description

This function generates data from the BayesDiallel model. It does not currently allow for sex interactions with strain-level effects, as is possible in the proper BayesDiallel package.

Usage

1
2
3
4
5
6
7
## S3 method for class 'diallel'
simulate(M = 8, mu = 10, inbred.mu = -5,
  female.mu = -2, add.size, inbred.size, epi.sym.size, epi.asym.size,
  maternal.size = 0, add.effect = NULL, inbred.effect = NULL,
  epi.sym.effect = NULL, epi.asym.effect = NULL,
  maternal.effect = NULL, n.each = 5, num.sim = 1,
  strains = LETTERS[1:M])

Arguments

M

DEFAULT: 8. The number of strains in the diallel cross. Default of 8 is the number of founders in the Collaborative Cross.

mu

DEFAULT: 10. The grand mean of the phenotype.

inbred.mu

DEFAULT: -5. The overall inbred effect.

female.mu

DEFAULT: -2. The overall female effect (male is reference).

add.size

The proportion of the variance explained by the additive effects.

inbred.size

The proportion of the variance explained by the inbred effects.

epi.sym.size

The proportion of the variance explained by the symmetric epistatic effects.

epi.asym.size

The proportion of the variance explained by the asymmetric epistatic effects.

maternal.size

The proportion of the variance explained by the maternal effects.

n.each

DEFAULT: 5. The number of individuals simulated per cell of the diallel. Total sample size of diallel cross will be n.each * 64

num.sim

DEFAULT: 1. The number of diallel cross simulations to perform.

strains

DEFAULT: c("A", "B", "C", "D", "E", "F", "G", "H"). The strain names used. Defaults to simple letter representation.

Examples

1

gkeele/DIDACT documentation built on Dec. 9, 2018, 12:27 a.m.