elliptical.diag.plots: Diagnostic Plots for Elliptical Regression Models

Description

This function generate diagnostic measures plots for the fitted elliptical regression models.

Usage

elliptical.diag.plots(
  object,
  ellipticaldiag = NULL,
  which,
  subset = NULL,
  iden = FALSE,
  labels = NULL,
  ret = FALSE,
  ...
)

Arguments

object

an object with the result of the fitted elliptical regression model.

ellipticaldiag

objects containing the diagnostic measures. If is NULL (by default) is obtained from object.

which

an optional numeric value with the number of plot returned.

subset

an optional numeric vector specifying a subset of observations to be used in the fitting process.

iden

a logical value used to identify observations. If TRUE the observations are identified in the graphic window.

labels

a optinal string vector specifying a labels plots.

ret

a logical value to return the diagnostic measures computing. If FALSE (by default) not return the diagnostic measures.

...

graphics parameters to be passed to the plotting routines.

Value

Return an interactive menu with eleven options to make plots. This menu contains the follows graphics : 1: plot: All. 2: plot: Response residual against fitted values. 3: plot: Response residual against index. 4: plot: Standardized residual against fitted values. 5: plot: Standardized residual against index. 6: plot: QQ-plot of response residuals. 7: plot: QQ-plot of Standardized residuals. 8: plot: Generalized Leverage. 9: plot: Local influence on the response against index. 10: plot: Local influence on the scale against index. 11: plot: Local influence for case-weight against index. If which is provided return an unique graphic selected. If ret is TRUE returns a list of diagnostic arrays (see elliptical.diag for more details).

References

Galea, M., Paula, G. A., and Cysneiros, F. J. A. (2005). On diagnostics in symmetrical nonlinear models. Statistics & Probability Letters, 73(4), 459-467. https://doi.org/10.1016/j.spl.2005.04.033

Examples

Run this code

# NOT RUN {
data(luzdat)
y <- luzdat$y
x1 <- luzdat$x1 ; x1 <- factor(x1) ; x1 <- C(x1,treatment)
x2 <- luzdat$x2
x3 <- (luzdat$x2)^2
luz <- data.frame(y,x1,x2,x3)
elliptical.fitt <- elliptical(y ~ x1+x2+x3, family = Student(df=5),
data=luz)
elliptical.diag.plots(elliptical.fitt, which=3)
# }