diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2020-10-22 13:14:31 +0200 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2020-10-22 13:17:40 +0200 |
| commit | 0af6a61b84cc29cdbfad16a6fc7ee0e6f88c7d0f (patch) | |
| tree | 7cb2715570688cac8c94b58e4ef78ba8afa29cc3 /R/nlme.mmkin.R | |
| parent | cb3adccaa4cdebc36c506fcb76c0d7c1d7e4e455 (diff) | |
Really fix check for nlme::varConstProp
And add output for nlme fit translating the mkinfit error model "obs"
into nlme::varIdent().
Diffstat (limited to 'R/nlme.mmkin.R')
| -rw-r--r-- | R/nlme.mmkin.R | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/R/nlme.mmkin.R b/R/nlme.mmkin.R index 7f7e34e9..c8a99d59 100644 --- a/R/nlme.mmkin.R +++ b/R/nlme.mmkin.R @@ -101,7 +101,7 @@ get_deg_func <- function() { #' endpoints(f_nlme_sfo_sfo) #' endpoints(f_nlme_dfop_sfo) #' -#' if (findFunction("varConstProp")) { # tc error model for nlme available +#' if (length(findFunction("varConstProp")) > 0) { # tc error model for nlme available #' # Attempts to fit metabolite kinetics with the tc error model #' #f_2_tc <- mmkin(list("SFO-SFO" = m_sfo_sfo, #' # "SFO-SFO-ff" = m_sfo_sfo_ff, @@ -126,6 +126,7 @@ get_deg_func <- function() { #' "DFOP-SFO" = m_dfop_sfo), #' ds_2, quiet = TRUE, error_model = "obs") #' f_nlme_sfo_sfo_obs <- nlme(f_2_obs["SFO-SFO", ]) +#' print(f_nlme_sfo_sfo_obs) #' # The same with DFOP-SFO does not converge, apparently the variances of #' # parent and A1 are too similar in this case, so that the model is #' # overparameterised |