rhalfsib: Simulate a dataset from a half-sib design

View source: R/rhalfsib.r

rhalfsibR Documentation

Simulate a dataset from a half-sib design

Description

Sample responses from the following model:

y_{ijk} = μ + α_i + β_{ij} + ε_{ijk},

where μ is a deterministic mean vector and each α_i, β_{ij}, ε_{ijk} are independent mean-zero multivariate Gaussians with covariance matrices A, B, E respectively.

Usage

rhalfsib(mu, A, I, B, J, E, K)

Arguments

mu

Mean vector.

A, B, E

Covariance matrices of the sires dams and individuals.

I, J, K

Number of sires, dams per sire and individuals per dam.

Value

A tibble with column names trait, value, sire, dam, individual. The value column contains the response values whose dimension is indexed by trait.


damian-t-p/halfsibdesign documentation built on March 14, 2023, 4:55 a.m.