example_2way | R Documentation |
Simulated example of a 2 way interaction GxE model (where G and E are latent variables).
g_j \sim Binomial(n=1,p=.30)
j = 1, 2, 3, 4
e_l \sim Normal(\mu=0,\sigma=1.5)
l = 1, 2, 3
g = .2g_1 + .15g_2 - .3g_3 + .1g_4 + .05g_1g_3 + .2g_2g_3
e = -.45e_1 + .35e_2 + .2e_3
\mu = -1 + 2g + 3e + 4ge
y \sim Normal(\mu=\mu,\sigma=\code{sigma}) if logit =FALSE |
y \sim Binomial(n=1,p=logit(\mu)) if logit =TRUE
|
example_2way(N, sigma = 1, logit = FALSE, seed = NULL)
N |
Sample size. |
sigma |
Standard deviation of the gaussian noise (if |
logit |
If TRUE, the outcome is transformed to binary with a logit link. |
seed |
RNG seed. |
Returns a list containing, in the following order: data.frame with the observed outcome (with noise) and the true outcome (without noise), data.frame of the genetic variants (G), data.frame of the environments (E), vector of the true genetic coefficients, vector of the true environmental coefficients, vector of the true main model coefficients
example_2way(5,1,logit=FALSE)
example_2way(5,0,logit=TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.