Usage
lm.case(fit, lms = summary.lm(fit), lmi = lm.influence(fit))
plot.case(x, fit,
which=c("stu.res","si","h","cook","dffits",
dimnames(x)[[2]][-(1:8)]), ##DFBETAS
between.in=list(y=4, x=9),
oma=c(0,0,0,4), cex.threshold=if.R(r=1.2, s=1.6),
main.in=list(
paste(deparse(fit$call), collapse=""),
cex=main.cex),
sigma.in=summary.lm(fit)$sigma,
p.in=summary.lm(fit)$df[1]-1,
obs.large=".lm.case.large",
obs.large.env=if.R(r=globalenv(), s=0),
main.cex=NULL,
...)
panel.case(x, y, subscripts, rownames, group.names,
nn, pp, ss, cex.threshold,
par.settings, ## R only. S-Plus ignores this argument
obs.large, obs.large.env,
...)
Arguments
fit
"lm"
object computed with x=TRUE
x
In plot.case
, the matrix output from lm.case
containing case diagnostics on each observation in the original
dataset.
In panel.case
, the x variable to be plotted
which
In plot.case
, the names of the columns of x
that are to be graphed.
between.in
between
trellis/lattice argument.
oma
In S-Plus, change par()$oma
to make room for the
threshold values. A warning is printed when par()$oma
is changed as the delayed printing of trellis objects implies we can't
return it to the original value automaticall
cex.threshold
Multiplier for cex
for the threshold values.
main.in
main
title for xyplot
. The default main title
displays the linear model formula from fit
.
sigma.in
standard error for the fit
.
p.in
The number of degrees of freedom associated with the fitted
model.
obs.large
Object name where the names of
all observations for which the standard
thresholds are exceeded will be stored. The default name is
.lm.case.large
.
obs.large.env
Frame in S-Plus (defaults to 0) and environment in
R (defaults to globalenv()
) where obs.large
will be stored.
main.cex
cex
for main title.
...
other arguments to xyplot
y
the y variable to be plotted.
nn
number of rows in original dataset.
pp
The number of degrees of freedom associated with the fitted model.
ss
standard error for the fit
.
subscripts
trellis/lattice argument, position in the reshaped
dataset constructed by plot.case
before calling xyplot
.
rownames
row name in the original data.frame.
group.names
names of the individual statistics.
par.settings
Used in R as part of the call to xyplot
.
Although this argument is not used in the panel function,
it is needed as a formal argument in S-Plus to absorb it out of ...
and thereby prevent it from being forwarded to