From b9b3c72cf076ce0393aa53a4d723a195b856c99c Mon Sep 17 00:00:00 2001
From: Johannes Ranke <jranke@uni-bremen.de>
Date: Mon, 7 Nov 2022 05:58:50 +0100
Subject: Make print and summary for saem.mmkin more robust

If the likelihood computed by importance sampling, these methods failed.
Now they report "Not available" or NA, respectively.
---
 R/saem.R               | 15 ++++++++++-----
 R/summary.saem.mmkin.R | 11 ++++++++---
 2 files changed, 18 insertions(+), 8 deletions(-)

diff --git a/R/saem.R b/R/saem.R
index 1c54ed93..09a4832f 100644
--- a/R/saem.R
+++ b/R/saem.R
@@ -257,11 +257,16 @@ print.saem.mmkin <- function(x, digits = max(3, getOption("digits") - 3), ...) {
     cat("\nFit did not terminate successfully\n")
   } else {
     cat("\nLikelihood computed by importance sampling\n")
-    print(data.frame(
-        AIC = AIC(x$so, type = "is"),
-        BIC = BIC(x$so, type = "is"),
-        logLik = logLik(x$so, type = "is"),
-        row.names = " "), digits = digits)
+    ll <- try(logLik(x$so, type = "is"), silent = TRUE)
+    if (inherits(ll, "try-error")) {
+      cat("Not available\n")
+    } else {
+      print(data.frame(
+          AIC = AIC(x$so, type = "is"),
+          BIC = BIC(x$so, type = "is"),
+          logLik = logLik(x$so, type = "is"),
+          row.names = " "), digits = digits)
+    }
 
     cat("\nFitted parameters:\n")
     conf.int <- parms(x, ci = TRUE)
diff --git a/R/summary.saem.mmkin.R b/R/summary.saem.mmkin.R
index ea2f02e2..2754e9f0 100644
--- a/R/summary.saem.mmkin.R
+++ b/R/summary.saem.mmkin.R
@@ -176,9 +176,14 @@ summary.saem.mmkin <- function(object, data = FALSE, verbose = FALSE, distimes =
   object$verbose <- verbose
 
   object$fixed <- object$mmkin_orig[[1]]$fixed
-  object$AIC = AIC(object$so)
-  object$BIC = BIC(object$so)
-  object$logLik = logLik(object$so, method = "is")
+  ll <-try(logLik(object$so, method = "is"), silent = TRUE)
+  if (inherits(ll, "try-error")) {
+    object$logLik <- object$AIC <- object $BIC <- NA
+  } else {
+    object$logLik = logLik(object$so, method = "is")
+    object$AIC = AIC(object$so)
+    object$BIC = BIC(object$so)
+  }
 
   ep <- endpoints(object)
   if (length(ep$ff) != 0)
-- 
cgit v1.2.1