This function simulates phenotype data by random sampling of markers available on `Glist`. Default parameters for the simulated phenotype reflect the genetic architecture assumed by BayesC prior(Habier et al., 2011). This function is under active development.
gsimC(Glist = NULL, h2 = NULL, m = NULL, prp.cau = NULL, n = NULL)
A list containing:
y
: Vector of simulated phenotypes.
g
: Vector of simulated genetic values.
e
: Vector of simulated residual effects.
b
: Vector of effect sizes of the simulated causal markers.
causal
: Vector of ids for the simulated causal markers.
h2
: Estimated heritability of the simulated phenotype.
A list containing genetic data. If NULL, the function will stop with an error.
Heritability. If NULL, heritability of 0.5 is assumed.
Number of causal markers. The values for either `m` or `prp.cau` should be provided at any given time. If the list of quality controlled markers is not available then list of raw markers is used. If `m` is NULL and `prp.cau` is also NULL, `prp.cau` will default to 0.001.
Proportion of causal markers. The values for either `m` or `prp.cau` should be provided at any given time.
Number of individuals randomly sampled from `Glist`. If NULL, all the individuals on `Glist` is used.
Peter Soerensen
Merina Shrestha