Simulates omicsarray data based on a given network.

# S4 method for omics_network
gene_expr_simulation(
  omics_network,
  time_label = 1:4,
  subject = 5,
  peak_level = 100,
  act_time_group = 1:4
)

Arguments

omics_network: A omics_network object.
time_label: a vector containing the time labels.
subject: the number of subjects
peak_level: the mean level of peaks.
act_time_group: [NULL] vector ; at which time the groups (defined by sort(unique(group))) are activated ?

Value

A omics_array object.

Author

Bertrand Frederic, Myriam Maumy-Bertrand.

Examples


data(Net)
set.seed(1)

#We simulate gene expressions according to the network Net
Msim<-Patterns::gene_expr_simulation(
  omics_network=Net,
  time_label=rep(1:4,each=25),
  subject=5,
  peak_level=200)
#> Loading required package: VGAM
#> Loading required package: stats4
#> Loading required package: splines
head(Msim)
#> $omicsarray
#>        log(S/US) : P1T1 log(S/US) : P1T2 log(S/US) : P1T3
#> gene 1         86.06709        44.533656       -57.361320
#> gene 2       -146.83138       120.514233       -39.892240
#> gene 3        228.34653        -3.625970       -60.889866
#> gene 4        505.11452        13.929252        -2.786049
#> gene 5        -36.57508        -1.828829        46.784308
#> gene 6       -486.82335       -91.502323      -173.402124
#> 
#> $name
#> [1] "gene 1" "gene 2" "gene 3" "gene 4" "gene 5" "gene 6"
#> 
#> $gene_ID
#> [1] 0
#> 
#> $group
#> [1] 1 1 1 1 1 1
#> 
#> $start_time
#> [1] 0
#> 
#> $time
#> [1] 1 2 3 4
#> 
#> $subject
#> [1] 5
#>