rmvtnorm.pedigree function

Simulate residual multivariate Gaussian data from a polygenic model