aboutsummaryrefslogtreecommitdiff
path: root/R/saem.R
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2022-11-04 07:48:04 +0100
committerJohannes Ranke <jranke@uni-bremen.de>2022-11-04 07:48:04 +0100
commit0389d523f049945c85eba42d4a006523595043f0 (patch)
treecd69527d328f8199700113c7f069302dc2586cbd /R/saem.R
parent48471ca62aeb1aded8aee7fd608607343c90955b (diff)
Never fail with errors in saem.mmkin
I think it is always preferable to be able to investigate the saemix object returned.
Diffstat (limited to 'R/saem.R')
-rw-r--r--R/saem.R6
1 files changed, 0 insertions, 6 deletions
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

Contact - Imprint