This computes the maximum of an object of class SLik representing an inferred (summary) likelihood surface
MSL(object, CIs = prod(dim(object$logLs)) < 12000L, level = 0.95,
verbose = interactive(),
RMSE_n=Infusion.getOption("RMSE_nsim"),
eval_RMSEs=(RMSE_n>1L) * prod(dim(object$logLs))The object is returned invisibly, with the following (possibly) added members, each of which being (as from version 1.5.0) an environment:
MSLcontaining variables MSLE and maxlogL that match the par and value returned by an optim call. Also contains the hessian of summary likelihood at its maximum.
CIobjectThe return value of calling allCIs, converted to an environment.
RMSEscontaining, as variable RMSEs, the root mean square errors of the log-likelihood at its inferred maximum and of the log-likelihood ratios at the CI bounds.
par_RMSEscontaining, as variable par_RMSEs, root mean square errors of the CI bounds.
To ensure backward-compatibility of code to possible future changes in the structure of the objects, the extractor function get_from should be used to extract the RMSEs and par_RMSEs variables from their respective environments, and more generally to extract any element from the objects.
an object of class SLik_j as produced by infer_SLik_joint (or, in the primitive workflow, of class SLik as produced by infer_surface.logLs).
If TRUE, construct one-dimensional confidence intervals for all parameters. See confint.SLik_j to obtain bootstrap confidence intervals.
Intended coverage probability of the confidence intervals.
Whether to display some information about progress and results.
Integer: number of simulation replicates for evaluation of prediction uncertainty for likelihoods/ likelihood ratios/ parameters. The default value (10) provides quick but inaccurate estimates.
Logical: whether to evaluate prediction uncertainty for likelihoods/ likelihood ratios/ parameters.
A list of arguments, passed to makeCluster, to control parallel computation of RMSEs. Beware that parallel computation of RMSEs tends to be memory-intensive. The list may contain a non-null spec element, in which case the nb_cores global Infusion option is ignored. Do *not* use a structured list with an RMSE element as is possible for refine (see Details of refine documentation).
Integer: shortcut for specifying cluster_args$spec.
Initial value for the optimiser. Better ignored.
(effective only for up-to-date workflow using gaussian mixture modelling of a joint distribution of parameters and statistics) a function that returns a vector of prior log-likelihood values, which is then added to the likelihood deduced from the summary likelihood analysis. The function's single argument must handle a matrix similar to the newdata argument of predict.SLik_j.
If Kriging has been used to construct the likelihood surface, RMSEs are computed using approximate formulas for prediction (co-)variances in linear mixed midels (see Details in predict). Otherwise, a more computer-intensive bootstrap method is used.
par_RMSEs are computed from RMSEs and from the numerical gradient of profile log-likelihood at each CI bound. Only RMSEs, not par_RMSEs, are compared to precision.
## see main documentation page for the package
Run the code above in your browser using DataLab