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.

1 files changed, 10 insertions, 5 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)