rand_data: Generate random data for test purposes

Description Usage Arguments Value Examples

View source: R/ctools.R

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)

cpgen documentation built on May 29, 2017, 9:40 p.m.

Search within the cpgen package
Search all R packages, documentation and source code