aboutsummaryrefslogtreecommitdiff
path: root/R/saem.R
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2022-10-13 03:48:54 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2022-10-13 03:50:05 +0200
commit37bd36fe8a75163cbf0f97cb7a0e2f7466a53617 (patch)
tree96bbd9b61006731b3d295e517def3b66ecd959ef /R/saem.R
parente7e8105390ebf3d6f034811bc7cce1d9640b7357 (diff)
Cope with failed FIM inversions
Diffstat (limited to 'R/saem.R')
-rw-r--r--R/saem.R3
1 files changed, 2 insertions, 1 deletions
diff --git a/R/saem.R b/R/saem.R
index 05cce682..99712c92 100644
--- a/R/saem.R
+++ b/R/saem.R
@@ -160,7 +160,7 @@ saem.mmkin <- function(object,
if (!fit_failed) {
if (any(is.na(f_saemix@results@se.fixed))) FIM_failed <- c(FIM_failed, "fixed effects")
if (any(is.na(c(f_saemix@results@se.omega, f_saemix@results@se.respar)))) {
- FIM_failed <- c(FIM_failed, "random effects and residual error parameters")
+ FIM_failed <- c(FIM_failed, "random effects and error model parameters")
}
if (!is.null(FIM_failed) & fail_with_errors) {
stop("Could not invert FIM for ", paste(FIM_failed, collapse = " and "))
@@ -208,6 +208,7 @@ saem.mmkin <- function(object,
so = f_saemix,
call = call,
time = fit_time,
+ FIM_failed = FIM_failed,
mean_dp_start = attr(m_saemix, "mean_dp_start"),
bparms.fixed = object[[1]]$bparms.fixed,
data = return_data,

Contact - Imprint