simulate_data_3gr: Simulate differentially linearly coexpressed genes between...

View source: R/simulate_data_3gr.R

simulate_data_3grR Documentation

Simulate differentially linearly coexpressed genes between three groups.

Description

Simulates sets of genes coexpressed with an eigengene which show differential coexpression among three groups. This tool simulates 9 possible classes of possible gene coexpression changes (positive, negative, null correlation in group A, group B, and groupC) but does not test all 27 possible classes.

Usage

simulate_data_3gr(samples_per_group = 50, total_genes = 2000,
  maxCor_true = 0.95, minCor_true = 0.3, maxCor_null = 0.2,
  verbose = T, corr_slope = 1, corr_intercept = 0,
  group1_offset = 0, group2_offset = 0, group3_offset = 0)

Arguments

samples_per_group

The number of samples per group to simulate. Default=50

total_genes

The number of genes among 9 differential coexpression classes to simulate. Will actually return a number of genes which is a multiple of 9 that is less than this number.

maxCor_true

The absolute value of the maximum correlation rho value for highly correlated pairs.

minCor_true

The absolute value of the minimum correlation rho value for highly correlated pairs.

maxCor_null

The absolute value of the maximum correlation rho value for null correlated pairs. The minimum is 0 correlation.

verbose

Whether or not to print additional information about the simulation when it runs.

corr_slope

The slope of the correlations to simulate. All correlated pairs will have this slope.

corr_intercept

The correlation intercept of the genes in matrix B.

group1_offset

The correlation intercept of group1 (e.g. its mean expression value)

group2_offset

The correlation intercept of group2 (e.g. its mean expression value)

group3_offset

The correlation intercept of group3 (e.g. its mean expression value)

Value

Returns a list of inputs for chowCor, including a matA, matB, design_mat, subgroup names (conditions), and the actual classes of differential coexpression.


ryananeff/superNOVA documentation built on March 29, 2024, 5:31 p.m.