symmetryDetection: Search for symmetries in the loaded model

In dMod: Dynamic Modeling and Parameter Estimation in ODE Models

Description

This function follows the method published in [1].

The function calls a python script via rPython. Usage problems might occur when different python versions are used. The script was written and tested for python 2.7.12, sympy 0.7.6.

Recently, users went into problems with RJSONIO when rPython was used. Unless a sound solution is available, please try to reinstall RJSONIO in these cases.

Usage

symmetryDetection(
  f,
  obsvect = NULL,
  prediction = NULL,
  initial = NULL,
  ansatz = "uni",
  pMax = 2,
  inputs = NULL,
  fixed = NULL,
  cores = 1,
  allTrafos = FALSE
)

Arguments

f	object containing the ODE for which as.eqnvec() is defined
obsvect	vector of observation functions
prediction	vector containing prediction to be tested
initial	vector containing initial values
ansatz	type of infinitesimal ansatz used for the analysis (uni, par, multi)
pMax	maximal degree of infinitesimal ansatz
inputs	specify the input variables
fixed	variables to concider fixed
cores	maximal number of cores used for the analysis
allTrafos	do not remove transformations with a common parameter factor

References

[1] https://journals.aps.org/pre/abstract/10.1103/PhysRevE.92.012920

Examples

## Not run: 
eq <- NULL
eq <- addReaction(eq, "A", "B", "k1*A")
eq <- addReaction(eq, "B", "A", "k2*B")

observables <- eqnvec(Aobs = "alpha * A")

symmetryDetection(eq, observables)


## End(Not run)

dMod documentation built on Jan. 27, 2021, 1:07 a.m.