WLinfer (version 1.2.0)

plot.WL: Some plots for WL

Description

plot method for a class "WL".

Usage

# S3 method for WL
plot(
  x,
  which = c(1, 2, 3, 4),
  ask = prod(par("mfcol")) < length(which) && dev.interactive(),
  ...
)

Arguments

x

an object of class "WL" made by the function WL.

which

if a subset of the plots is required, specify a subset of 1:4; see 'Details' for a description of the plots.

ask

logical; if TRUE, the user is asked before each plot.

...

other parameters to be passed through to plotting functions.

Details

The first plot (which=1) is histogram of input values and pdf from estimated lambda and phi. This is a good plot to judge how well the estimated values match the actual data.

The second plot (which=2) is the boxplot of the input data. It shows the simple information about median and outlier points. Additionally, you can easily see which the "Outlier Points" are with box on the right side.

The third plot (which=3) is Q-Q plot. This provides the same information as the normal Q-Q plot.

The last plot (which=4) is contour plot of lambda and phi. You can observe the location of estimator, and also compare estimator to Cox&Snell bias corrected MLE and MLEc.

Plots are basically displayed one by one in the order mentioned above, and it is also possible to view only one plot or four at a time. How to do this is shown in the example below.

See Also

See also WL and plot.

Examples

example <- lifetime_alum
result <- WL(example)

plot(result)
plot(result, which=1)
plot(result, which=3)
par(mfrow = c(2, 2)) ; plot(result)
par(mfrow = c(1, 1))