MSL: Maximum likelihood from an inferred likelihood surface

Description Usage Arguments Details Value Examples

Description

This computes the maximum of an object of class SLik representing an inferred (summary) likelihood surface

Usage

1
2
MSL(object, CIs = TRUE, level = 0.95, verbose = interactive(),
    eval_RMSEs = TRUE, cluster_args=list(),...)

Arguments

object

an object of class SLik as produced by infer_surface.logLs

CIs

If TRUE, construct one-dimensional confidence intervals for all parameters.

level

Intended coverage probability of the confidence intervals.

verbose

Whether to display some information about progress and results.

eval_RMSEs

Logical: whether to evaluate prediction uncertainty for likelihoods/ likelihood ratios/ parameters.

cluster_args

A list of arguments, passed to makeCluster. May contain a non-null spec element, in which case the distinct nb_cores argument is ignored.

...

Further arguments passed from or to other methods.

Details

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.

Value

The object is returned invisibly, with the following added members, each of which being (as from version 1.5.0) an environment:

MSL

, containing variables MSLE and maxlogL that match the par and value returned by an optim call.

RMSEs

containing, 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_RMSEs

containing, 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.

Examples

1
## see main documentation page for the package

Infusion documentation built on Feb. 22, 2021, 9:08 a.m.

Related to MSL in Infusion...