simulate_GRN: Generate gene regulatory network (GRN) data from structural...

Description Usage Arguments Value Examples

View source: R/simulate_GRN.R

Description

Generate gene regulatory network (GRN) data from structural equation model.

Usage

1
2
3
4
5
6
7
8
9
simulate_GRN(
  ngen,
  nexp,
  nobs,
  prob_edge,
  seed = NULL,
  lwr_edge_weight = -1,
  upr_edge_weight = 1
)

Arguments

ngen

Integer number of genetic markers.

nexp

Integer number of expression traits.

nobs

Integer number of observations in the generated data set.

prob_edge

Numeric probability of having a causal interaction (edge) between two expression traits. Passed on to randomDAG.

seed

Integer seed for reproducible random number generation.

lwr_edge_weight

Lower limit of edge weight (causal strength T_i -> T_j)

upr_edge_weight

Upper limit of edge weight (causal strength T_i -> T_j)

Value

A data set consisting of nobs samples from a random generated gene regulatory network. This includes ngen binomially distributed variables representing the genetic markers associated with the genetic makeup of each individual observation. It also includes nexp columns representing the gene expression levels of nexp traits, which have an underlying causal graph describing their interactions.

Examples

1
simulate_GRN(10, 10, 100, 0.1, 1634)

igbucur/BFCS documentation built on Oct. 20, 2020, 7:39 a.m.