plot.MLEce {MLEce}

Providing some plots for effective closed-form estimators

Description

plot method for a class "MLEce".

Usage

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

Arguments

x

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

which

if a subset of the plots is required, specify a subset of 1:4.

ask

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

...

not used, but exists because of the compatibility.

Details

The boxplot for given data is presented first with which=1. For which=2, a contour line is drawn by the probability density function of the estimated parameter based on effective closed-form estimators. In the counter plot, the x-axis is the first column of data and the y-axis is the second column of data. For which=3, a marginally fitted probability density plot is given for the first column of input data. And a fitted line is added for the efficient closed-form estimator. For which=4, is a marginally fitted probability density plot is given like the former one for the second column of input data. Note that, marginally fitted probability density plots in which=3 and which=4 present comparisons between efficient closed form estimators (MLEces) and correlation based method estimators (CMEs) for the bivariate Weibull distribution. Note that this plot commend is limited at the bivariate distributions.

Examples

data(flood)
est_BiGam <- MLEce(flood, "BiGam")
plot(est_BiGam, c(3))
air_data <- airquality[ ,3:4]
air_data[ ,2] <- air_data[ ,2]*0.1
est_BiWei <- MLEce(air_data, "BiWei")
plot(est_BiWei)
data(fossil_pollen) 
fossil_data <- cbind(fossil_pollen[,1]/100,rowSums(fossil_pollen[,-1]/100))
est_fossil <- MLEce(fossil_data, "Dirichlet")
plot(est_fossil,c(2))