rand_data | R Documentation |
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).
rand_data(n=500,p_marker=10000,h2=0.3,prop_qtl=0.01,seed=NULL)
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 |
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
# Generate random data with 100 observations and 500 markers
rand_data(100,500)
# check that objects have been created
str(M)
str(y)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.