plot.fits.compare {RobustAFT} | R Documentation |
Comparative plots for objects of class "fits.compare".
## S3 method for class 'fits.compare' plot(x, xplots = FALSE, ask = TRUE, which = 1:4, leg.position = c("topleft", "topleft", "topleft"), ...)
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 |
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
.
fits.compare
, plot.default
, plot.TML
## 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)