example_3way_3latent: Simulated example of a 3 way interaction GxExZ model

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

Description

Simulated example of a 3 way interaction GxExZ model (where G, E and Z are latent variables).

g_j \sim Binomial(n=1,p=.30)

j = 1, 2, 3, 4

e_k \sim Normal(μ=0,σ=1.5)

k = 1, 2, 3

z_l \sim Normal(μ=3,σ=1)

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

z = .15z_1 + .60z_2 + .25z_3

μ = -2 + 2g + 3e + z + 5ge - 1.5ez + 2gz + 2gez

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

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

Usage

example_3way_3latent

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), list containing the data.frame of the genetic variants (G), the data.frame of the e environments (E) and the data.frame of the z environments (Z), vector of the true genetic coefficients, vector of the true e environmental coefficients, vector of the true z environmental coefficients, vector of the true main model coefficients

Examples

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

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