sim1: Simulate a single multivariate data set consisting of n...

In fboehm/qtl2pleio: Testing Pleiotropy in Multiparental Populations

Description

Simulate a single multivariate data set consisting of n subjects and d phenotypes for each

Usage

sim1(X, B, Sigma)

Arguments

X	design matrix (incorporating genotype probabilities from two loci), dn by df
B	a matrix of allele effects, f rows by d columns
Sigma	dn by dn covariance matrix

Value

a vector of length dn. The first n entries are for trait 1, the second n for trait 2, etc.

Examples

n_mouse <- 20
geno <- rbinom(n = n_mouse, size = 1, prob = 1 / 2)
X <- gemma2::stagger_mats(geno, geno)
B <- matrix(c(1, 2), ncol = 2, nrow = 1)
sim1(X, B, Sigma = diag(2 * n_mouse))

fboehm/qtl2pleio documentation built on July 17, 2021, 11:34 p.m.