rand_data: Generate random data for test purposes

In cheuerde/cpgen: Parallel Genomic Evaluations

Generate random data for test purposes

Description

Generates a random marker-matrix in {-1,0,1} coding and a phenotype vector. Phenotypic variance times h2 (variance explained by markers) is equally spread among all markers (sampled from uniform distribution).

Usage

rand_data(n=500,p_marker=10000,h2=0.3,prop_qtl=0.01,seed=NULL)

Arguments

n	Number of oberservations
p_marker	Number of markers
h2	Heritability of the trait
prop_qtl	Proportion of QTL of total number of markers
seed	Seed for RNG

Value

No return value. Generates two objects globally (M and y) that can be used after the execution of the function. M is the marker matrix and y the phenotype vector

Examples

# Generate random data with 100 observations and 500 markers
rand_data(100,500)

# check that objects have been created
str(M)
str(y)

cheuerde/cpgen documentation built on June 24, 2024, 2:33 p.m.