From 0af6a61b84cc29cdbfad16a6fc7ee0e6f88c7d0f Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Thu, 22 Oct 2020 13:14:31 +0200 Subject: Really fix check for nlme::varConstProp And add output for nlme fit translating the mkinfit error model "obs" into nlme::varIdent(). --- R/nlme.mmkin.R | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) (limited to 'R/nlme.mmkin.R') 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 -- cgit v1.2.1