From b9b3c72cf076ce0393aa53a4d723a195b856c99c Mon Sep 17 00:00:00 2001 From: Johannes Ranke 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 ++++++++++----- 1 file changed, 10 insertions(+), 5 deletions(-) (limited to 'R/saem.R') 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) -- cgit v1.2.1