example_2way: Simulated example of a 2 way interaction GxE model.

In AlexiaJM/LEGIT: Latent Environmental & Genetic InTeraction (LEGIT) Model

Description

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(μ=0,σ=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

μ = -1 + 2g + 3e + 4ge

y \sim Normal(μ=μ,σ=\code{sigma}) if logit=FALSE

y \sim Binomial(n=1,p=logit(μ)) if logit=TRUE

Usage

example_2way

Arguments

N	Sample size.
sigma	Standard deviation of the gaussian noise (if logit=FALSE).
logit	If TRUE, the outcome is transformed to binary with a logit link.
seed	RNG seed.

Format

An object of class function of length 1.

Value

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

Examples

example_2way(5,1,logit=FALSE)
example_2way(5,0,logit=TRUE)

AlexiaJM/LEGIT documentation built on Feb. 14, 2021, 8:24 p.m.