| plot.gvlmaDel {gvlma} | R Documentation |
Various Plots for a Gvlmadel Object
Description
Plots to display the behavior of the deletion statistics stored in a gvlmaDel object.
Usage
## S3 method for class 'gvlmaDel'
plot(x, which = 1:2, TukeyStyle = TRUE, ask
= prod(par("mfcol")) < max(c(10, 5)[which]) && dev.interactive(),
pointlabels, ...)
Arguments
x |
A |
which |
Vector indicating which, or both, of two types of plots to show. |
TukeyStyle |
If TRUE, determine unusual observations in a robust way based on inter-quartile ranges, else based on standard deviations. |
ask |
If TRUE, prompt the user before beginning a new page of graphs. |
pointlabels |
A vector of length the number of observations in
the linear model fit in the |
... |
Additional arguments that are ignored. |
Details
If which = 1, each of the 5 deletion statistics (deletion
global statistic and each of the 4 directional statistics) is plotted
against the time sequence used for the 4th directional statistic
(assessing heteroscedasticity).
If which = 2, the function display.delstats is called
for each of the 5 deletion statistics. The argument TukeyStyle
is passed directly to display.delstats. See the help for
display.delstats for details.
If which = c(1,2), the default, then all 10 plots are
generated.
The deletion statistics in the gvlmaDel object are the percent
relative change when each observation, in turn, is omitted from the
model fitting.
Value
No value is returned.
Author(s)
Slate, EH slate@stat.fsu.edu and Pena, EA pena@stat.sc.edu.
References
Pena, EA and Slate, EH (2006). “Global validation of linear model assumptions,” J.\ Amer.\ Statist.\ Assoc., 101(473):341-354.
See Also
Examples
data(CarMileageData)
CarModelAssess <- gvlma(NumGallons ~ MilesLastFill + NumDaysBetw,
data = CarMileageData)
CarModelDel <- deletion.gvlma(CarModelAssess)
par(mfrow=c(1,1))
plot(CarModelDel)
par(mfrow=c(2,2))
plot(CarModelDel)
plot(CarModelDel, TukeyStyle = FALSE)
plot(CarModelDel, which = 2)