From 0389d523f049945c85eba42d4a006523595043f0 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Fri, 4 Nov 2022 07:48:04 +0100 Subject: Never fail with errors in saem.mmkin I think it is always preferable to be able to investigate the saemix object returned. --- R/saem.R | 6 ------ 1 file changed, 6 deletions(-) (limited to 'R') diff --git a/R/saem.R b/R/saem.R index 440a187d..1c54ed93 100644 --- a/R/saem.R +++ b/R/saem.R @@ -45,8 +45,6 @@ utils::globalVariables(c("predicted", "std")) #' @param covariate_models A list containing linear model formulas with one explanatory #' variable, i.e. of the type 'parameter ~ covariate'. Covariates must be available #' in the 'covariates' data frame. -#' @param fail_with_errors Should a failure to compute standard errors -#' from the inverse of the Fisher Information Matrix be a failure? #' @param quiet Should we suppress the messages saemix prints at the beginning #' and the end of the optimisation process? #' @param nbiter.saemix Convenience option to increase the number of @@ -145,7 +143,6 @@ saem.mmkin <- function(object, control = list(displayProgress = FALSE, print = FALSE, nbiter.saemix = nbiter.saemix, save = FALSE, save.graphs = FALSE), - fail_with_errors = TRUE, verbose = FALSE, quiet = FALSE, ...) { call <- match.call() @@ -177,9 +174,6 @@ saem.mmkin <- function(object, if (any(is.na(c(f_saemix@results@se.omega, f_saemix@results@se.respar)))) { 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 ")) - } transparms_optim <- f_saemix@results@fixed.effects names(transparms_optim) <- f_saemix@results@name.fixed -- cgit v1.2.1