generate_kinetics() samples the kinetics of genes in the feature network for which the kinetics have not yet been defined. kinetics_default() is used to configure parameters pertaining this process. kinetics_random_distributions() will do the same, but the distributions are also randomised.

generate_kinetics(model)

kinetics_default()

kinetics_random_distributions()

Arguments

model

A dyngen intermediary model for which the feature network has been generated with generate_feature_network().

Value

A dyngen model.

Details

To write different kinetics settings, you need to write three functions with interface function(feature_info, feature_network, cache_dir, verbose). Described below are the default kinetics samplers.

sampler_tfs() mutates the feature_info data frame by adding the following columns:

  • transcription_rate: the rate at which pre-mRNAs are transcribed, in pre-mRNA / hour. Default distribution: U(1, 2).

  • translation_rate: the rate at which mRNAs are translated into proteins, in protein / mRNA / hour. Default distribution: U(100, 150).

  • mrna_halflife: the half-life of (pre-)mRNA molecules, in hours. Default distribution: U(2.5, 5).

  • protein_halflife: the half-life of proteins, in hours. Default distribution: U(5, 10).

  • splicing_rate: the rate at which pre-mRNAs are spliced into mRNAs, in reactions / hour. Default value: log(2) / (10/60), which corresponds to a half-life of 10 minutes.

  • independence: the degree to which all regulators need to be bound for transcription to occur (0), or whether transcription can occur if only one of the regulators is bound (1).

sampler_nontfs() samples the transcription_rate, translation_rate, mrna_halflife and protein_halflife from a supplementary file of Schwannhäusser et al., 2011, doi.org/10.1038/nature10098. splicing_rate is by default the same as in sampler_tfs(). independence is sampled from U(0, 1).

sampler_interactions() mutates the feature_network data frame by adding the following columns.

  • effect: the effect of the interaction; upregulating = +1, downregulating = -1. By default, sampled from -1, 1 with probabilities .25, .75.

  • strength: the strength of the interaction. Default distribution: 10^U(0, 2).

  • hill: the hill coefficient. Default distribution: N(2, 2) with a minimum of 1 and a maximum of 10.

See also

dyngen on how to run a complete dyngen simulation

Examples

model <- initialise_model( backbone = backbone_bifurcating(), kinetics_params = kinetics_default() ) # \donttest{ data("example_model") model <- example_model %>% generate_kinetics()
#> | | 0 % ~calculating |+ | 2 % ~00s |++ | 4 % ~00s |+++ | 5 % ~00s |++++ | 7 % ~00s |+++++ | 9 % ~00s |++++++ | 11% ~00s |+++++++ | 13% ~00s |++++++++ | 15% ~00s |+++++++++ | 16% ~00s |++++++++++ | 18% ~00s |++++++++++ | 20% ~00s |+++++++++++ | 22% ~00s |++++++++++++ | 24% ~00s |+++++++++++++ | 25% ~00s |++++++++++++++ | 27% ~00s |+++++++++++++++ | 29% ~00s |++++++++++++++++ | 31% ~00s |+++++++++++++++++ | 33% ~00s |++++++++++++++++++ | 35% ~00s |+++++++++++++++++++ | 36% ~00s |++++++++++++++++++++ | 38% ~00s |++++++++++++++++++++ | 40% ~00s |+++++++++++++++++++++ | 42% ~00s |++++++++++++++++++++++ | 44% ~00s |+++++++++++++++++++++++ | 45% ~00s |++++++++++++++++++++++++ | 47% ~00s |+++++++++++++++++++++++++ | 49% ~00s |++++++++++++++++++++++++++ | 51% ~00s |+++++++++++++++++++++++++++ | 53% ~00s |++++++++++++++++++++++++++++ | 55% ~00s |+++++++++++++++++++++++++++++ | 56% ~00s |++++++++++++++++++++++++++++++ | 58% ~00s |++++++++++++++++++++++++++++++ | 60% ~00s |+++++++++++++++++++++++++++++++ | 62% ~00s |++++++++++++++++++++++++++++++++ | 64% ~00s |+++++++++++++++++++++++++++++++++ | 65% ~00s |++++++++++++++++++++++++++++++++++ | 67% ~00s |+++++++++++++++++++++++++++++++++++ | 69% ~00s |++++++++++++++++++++++++++++++++++++ | 71% ~00s |+++++++++++++++++++++++++++++++++++++ | 73% ~00s |++++++++++++++++++++++++++++++++++++++ | 75% ~00s |+++++++++++++++++++++++++++++++++++++++ | 76% ~00s |++++++++++++++++++++++++++++++++++++++++ | 78% ~00s |++++++++++++++++++++++++++++++++++++++++ | 80% ~00s |+++++++++++++++++++++++++++++++++++++++++ | 82% ~00s |++++++++++++++++++++++++++++++++++++++++++ | 84% ~00s |+++++++++++++++++++++++++++++++++++++++++++ | 85% ~00s |++++++++++++++++++++++++++++++++++++++++++++ | 87% ~00s |+++++++++++++++++++++++++++++++++++++++++++++ | 89% ~00s |++++++++++++++++++++++++++++++++++++++++++++++ | 91% ~00s |+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s |++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s |+++++++++++++++++++++++++++++++++++++++++++++++++ | 96% ~00s |++++++++++++++++++++++++++++++++++++++++++++++++++| 98% ~00s |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s
# }