Generate random data for test purposes

Share:

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

1
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

1
2
3
4
5
6
# Generate random data with 100 observations and 500 markers
rand_data(100,500)

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