DCP_sim_data: Simulating Data for Rhythmicity Comparison

In DiffCircaPipeline/Rpackage: An R package for DiffCircaPipeline: A framework for multifaceted characterization of differential rhythmicity

Simulating Data for Rhythmicity Comparison

Description

DCP_sim_data simulates a simple two-group study. The simulated data is just used for illustration purposes and is not intended to mimic real data.

Usage

DCP_sim_data(
  ngene = 1000,
  nsample = 30,
  A1 = c(1, 3),
  A2 = c(1, 3),
  phase1 = c(0, 2 * pi),
  phase2 = c(0, 2 * pi),
  M1 = c(4, 6),
  M2 = c(4, 6),
  sigma1 = 1,
  sigma2 = 1
)

Arguments

ngene	integer. Number of features. Shared by both groups
nsample	integer. Number of samples Shared by both groups
A1	vector of length 2. c(min, max) amplitude for group I
A2	vector of length 2. c(min, max) amplitude for group II
phase1	vector of length 2. c(min, max) phase for group I (values not in the 0 to 2π interval will be converted)
phase2	vector of length 2. c(min, max) phase for group II
M1	vector of length 2. c(min, max) MESOR for group I
M2	vector of length 2. c(min, max) MESOR for group I
sigma1	numeric. Noise level for group I.
sigma2	numeric. Noise level for group II.

Value

A list of length two, each is data for one group.

Examples

#simulate data with shifted phase
x = DCP_sim_data(ngene=1000, nsample=30, A1=c(1, 3), A2=c(1, 3),
phase1=c(0, pi/4), phase2=c(pi/4, pi/2),
M1=c(4, 6), M2=c(4, 6), sigma1=1, sigma2=1)

DiffCircaPipeline/Rpackage documentation built on March 17, 2023, 7:32 a.m.