Infusion (version 2.2.0)

plot.SLik: Plot SLik or SLikp objects

Description

Primarily conceived for exposition purposes, for the two-parameters case. The black-filled points are those for which the observed summary statistic was outside of the convex hull of the simulated empirical distribution. The crosses mark the estimated ML point and the confidence intervals points, that is, the outmost points on the contour defined by the profile likelihood threshold for the profile confidence intervals. There is a pair of CI points for each interval. The smaller black dots mark points added in the latest iteration, if refine was used.

Usage

# S3 method for SLik
plot(x, y, filled = FALSE, decorations = NULL,
                    color.palette = NULL, plot.axes = NULL, 
                    plot.title = NULL, plot.slices=TRUE, ...)
# S3 method for SLik_j
plot(x, y, filled = nrow(x$logLs)>5000L, decorations = NULL, 
                      color.palette = NULL, plot.axes = NULL, 
                      plot.title = NULL, from_refine=FALSE, plot.slices=TRUE, 
                      show_latest=FALSE, ...)

Value

plot.SLik_j returns invisibly a list including coordinates of the plot(s) (at least if the latest version of the spaMM package is installed). The exact format will depend on the nature of the plot but the names of elements should be self-explanatory. plot.SLik returns the plotted object invisibly.

Arguments

x

An object of class SLik or SLikp

y

Not used, but included for consistency with the plot generic.

filled

whether to plot a mapMM or a filled.mapMM.

decorations

Graphic directives added to the default decorations value in calls of mapMM or a filled.mapMM (see the source code of plot.SLik for the latter default values).

color.palette

Either NULL or a function that can replace the default color function used by plot.SLik. The function must have a single argument, giving the number of color levels.

plot.title

statements which replace the default titles to the main plot (see Details).

plot.axes

statements which replace the default axes on the main plot (see Details).

from_refine

For programming purposes, not documented.

plot.slices

boolean: whether to plot “slices” of the summary-likelihood surface for pairs of parameters (p1,p2), when more than two parameters are fitted. In such plots the additional parameters p3, p4... are fixed to their estimates [in contrast to profile plots where p3, p4... take distinct values for each (p1,p2), maximizing the function for each (p1,p2)].

show_latest

Logical: whether to show distinctly the points added in the latest iteration.

...

further arguments passed to or from other methods (currently can be used to pass a few arguments such as map.asp in all cases, or variances to filled.mapMM).

Details

Different graphic functions are called depending on the number of estimated parameters. For two parameters, mapMM or filled.mapMM are called. For more than two parameters, spaMM.filled.contour is called. See the documentation of these functions for the appropriate format of the plot.title and plot.axes arguments.

Examples

if (FALSE) {
## Using 'slik_j' object from the example in help("example_reftable") 
plot(slik_j,filled=TRUE,
     plot.title=quote(title("Summary-likelihood-ratio surface",
                             xlab=expression(mu),
                             ylab=expression(sigma^2))))
}