[R] Support for compartment models / matrix exponential method

[Thomas Hoffmann]

I try to fit reaction kinetic data to a given hypothetical model
consisting of 3 ODEs. After some reading in Bates/Watts and the literature
cited there (e.g. Jennrich/Bright) I looked in R for functions that support
the use of the matrix exponential method. Until now I used for a similar
(but simpler) problem the nls package and fitted against the the parameterized
solution function. But this works only for simple kinetics.

The algorithms in the nls package (e.g. SSfol) work on "closed forms" of the
solution function, too.

My question is: Are there examples or functions in R that show the usage of
the matrix exponential method or do I have to build the steps from primitive
functions myself, e.g. following the recipe in Jennrich/Bright (or is there 
even a superior method for this kind of tasks that I am not aware of)?


=============================================================================
Thomas Hoffmann, Institut fuer Halbleiter- und Mikrosystemtechnik, TU Dresden 
E-mail: hoffmann at ehmgs2.et.tu-dresden.de
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._