plot.fits.compare {RobustAFT} | R Documentation |
Plot Method for "fits.compare" objects
Description
Comparative plots for objects of class "fits.compare".
Usage
## S3 method for class 'fits.compare'
plot(x, xplots = FALSE, ask = TRUE, which = 1:4,
leg.position = c("topleft", "topleft", "topleft"), ...)
Arguments
x |
An object of class "fits.compare", usually, a result of a call to |
xplots |
If xplots=TRUE, plots of the independent variables versus the residuals are produced. |
ask |
If ask=TRUE, plot.fits.compare() operates in interactive mode. |
which |
If a subset of the plots is required, specify a subset of the numbers |
leg.position |
A vector of character string specifying the legend position of the second, third and fourth plots. |
... |
Optional arguments for |
Details
For clarity reasons, at most three models should be compared.
Four default plots (selectable by which) are produced: histograms of the residuals of each model,
a residual Q-Q plot, response against fitted values and residuals against fitted values.
Additional plots are produced if xplots=TRUE
.
See Also
fits.compare
, plot.default
, plot.TML
Examples
## Not run:
data(D243)
Cost <- D243$Cost # Cost (Swiss francs)
LOS <- D243$LOS # Length of stay (days)
Adm <- D243$Typadm; Adm <- (Adm==" Urg")*1 # Type of admission
# (0=on notification, 1=Emergency)
lwrob <- TML.noncensored(log(Cost)~log(LOS)+Adm, errors="logWeibull")
reg <- lm(log(Cost)~log(LOS)+Adm)
comp <- fits.compare(least.squares=reg, TML.logWeibull=lwrob)
plot(comp, leg.position=c("topleft", "topleft", "bottomleft"), xplots=TRUE)
## End(Not run)