ellipse(center, shape, radius, log="", center.pch=19, center.cex=1.5,
segments=51, draw=TRUE, add=draw, xlab="", ylab="",
col=palette()[2], lwd=2, fill=FALSE, fill.alpha=0.3, grid=TRUE, ...)
dataEllipse(x, y, weights, log = "", levels = c(0.5, 0.95), center.pch = 19,
center.cex = 1.5, draw = TRUE, plot.points = draw, add = !plot.points,
segments = 51, robust = FALSE, xlab = deparse(substitute(x)),
ylab = deparse(substitute(y)), col = palette()[1:2], lwd = 2,
fill = FALSE, fill.alpha = 0.3, grid = TRUE, labels, id.method = "mahal",
id.n = if (id.method[1] == "identify") Inf else 0, id.cex = 1,
id.col = palette()[1], ...)
confidenceEllipse(model, ...)
## S3 method for class 'lm':
confidenceEllipse(model, which.coef, L, levels=0.95, Scheffe=FALSE, dfn,
center.pch=19, center.cex=1.5, segments=51, xlab, ylab,
col=palette()[2], lwd=2, fill=FALSE, fill.alpha=0.3, draw=TRUE, add=!draw, ...)
## S3 method for class 'glm':
confidenceEllipse(model, chisq, ...)
## S3 method for class 'default':
confidenceEllipse(model, which.coef, L, levels=0.95, Scheffe=FALSE, dfn,
center.pch=19, center.cex=1.5, segments=51, xlab, ylab,
col=palette()[2], lwd=2, fill=FALSE, fill.alpha=0.3, draw=TRUE, add=!draw, ...)
"x"
if the x-axis is logged, "y"
if the y-axis is
logged, and "xy"
FALSE
or NULL
the center point isn't plotted.TRUE
produce graphical output; if FALSE
, only invisibly return coordinates of ellipse(s).TRUE
add ellipse(s) to current plot.y
is missing) a 2-column numeric matrix.x
.FALSE
data ellipses are drawn,
but points are not plotted.TRUE
use the cov.trob
function in the lm
or glm
.L
matrix is given, it takes precedence over the which.coef
argTRUE
scale the ellipse so that its projections onto the
axes give Scheffe confidence intervals for the coefficients.Scheffe
is TRUE
, or to
TRUE
, the confidence ellipse for the coefficients in a generalized linear model are
based on the chisquare statistic, if FALSE
on the $F$-statistic. This corresponds to using the default
and linear-model methods resp2
(see par
).col
(default, FALSE
)?0.3
).plot
and
line
.id.n=0
for labeling no points. See
showLabels
for details of these arguments.ellipse
returns invisibly the (x, y) coordinates of the calculated ellipse.
dataEllipse
and confidenceEllipse
return invisibly the coordinates of one or more ellipses, in the latter instance a list named by
levels
.dataEllipse
superimposes the normal-probability contours over a scatterplot
of the data.cov.trob
, cov.wt
, linearHypothesis
.dataEllipse(Duncan$income, Duncan$education, levels=0.1*1:9, lty=2,
fill=TRUE, fill.alpha=0.1)
confidenceEllipse(lm(prestige~income+education, data=Duncan), Scheffe=TRUE)
confidenceEllipse(lm(prestige~income+education, data=Duncan),
L=c("income + education", "income - education"))
wts <- rep(1, nrow(Duncan))
wts[c(6, 16)] <- 0 # delete Minister, Conductor
with(Duncan, {
dataEllipse(income, prestige, levels=0.68)
dataEllipse(income, prestige, levels=0.68, robust=TRUE, plot.points=FALSE, col="green3")
dataEllipse(income, prestige, weights=wts, levels=0.68, plot.points=FALSE, col="brown")
dataEllipse(income, prestige, weights=wts, robust=TRUE, levels=0.68,
plot.points=FALSE, col="blue")
})
Run the code above in your browser using DataCamp Workspace