Automatic starting parameters for saem.mmkin

For the case of mkin transformations. This gives faster convergence,
and appears to avoid problems with numeric ODE solutions

1 files changed, 15 insertions, 1 deletions

diff --git a/man/saem.Rd b/man/saem.Rd
index 3dc8001a..11463351 100644
--- a/man/saem.Rd
+++ b/man/saem.Rd
@@ -20,9 +20,11 @@ saem(object, ...)
   conf.level = 0.6,
   solution_type = "auto",
   covariance.model = "auto",
+  omega.init = "auto",
   covariates = NULL,
   covariate_models = NULL,
   no_random_effect = NULL,
+  error.init = c(3, 0.1),
   nbiter.saemix = c(300, 100),
   control = list(displayProgress = FALSE, print = FALSE, nbiter.saemix = nbiter.saemix,
     save = FALSE, save.graphs = FALSE),
@@ -41,8 +43,10 @@ saemix_model(
   degparms_start = numeric(),
   covariance.model = "auto",
   no_random_effect = NULL,
+  omega.init = "auto",
   covariates = NULL,
   covariate_models = NULL,
+  error.init = numeric(),
   test_log_parms = FALSE,
   conf.level = 0.6,
   verbose = FALSE,
@@ -81,10 +85,18 @@ for parameter that are tested if requested by 'test_log_parms'.}
 \item{solution_type}{Possibility to specify the solution type in case the
 automatic choice is not desired}
 
-\item{covariance.model}{Will be passed to \code{\link[saemix:SaemixModel-class]{saemix::SaemixModel()}}. Per
+\item{covariance.model}{Will be passed to \code{\link[saemix:saemixModel]{saemix::saemixModel()}}. Per
 default, uncorrelated random effects are specified for all degradation
 parameters.}
 
+\item{omega.init}{Will be passed to \code{\link[saemix:saemixModel]{saemix::saemixModel()}}. If using
+mkin transformations and the default covariance model with optionally
+excluded random effects, the variances of the degradation parameters
+are estimated using \link{mean_degparms}, with testing of untransformed
+log parameters for significant difference from zero. If not using
+mkin transformations or a custom covariance model, the default
+initialisation of \link[saemix:saemixModel]{saemix::saemixModel} is used for omega.init.}
+
 \item{covariates}{A data frame with covariate data for use in
 'covariate_models', with dataset names as row names.}
 
@@ -96,6 +108,8 @@ in the 'covariates' data frame.}
 which there should be no variability over the groups. Only used
 if the covariance model is not explicitly specified.}
 
+\item{error.init}{Will be passed to \code{\link[saemix:saemixModel]{saemix::saemixModel()}}.}
+
 \item{nbiter.saemix}{Convenience option to increase the number of
 iterations}