earth
models by plotting them.plot.earth.models(x = stop("no 'x' arg"), which = c(1:2),
caption = "", rlim = c(0,1), jitter = 0,
col.grsq = discrete.plot.cols(length(x)),
col.rsq = 0, col.npreds = 0, col.cum = NULL,
col.vline = "grey", lty.vline = 3, col.legend = 1,
legend.pos = NULL, legend.text = NULL, do.par = TRUE,
main = "Model Comparison", ...)
1:2
, meaning both."string"
string
""
(default) no caption
NULL
generate a caption from the $call
component of the earth
objects.c(min,max)
specifying min and max
values on the y axis of the RSq/GRSq plot.
Default is c(0,1)
. The special value min=-1
means the minimum y axis value
is the smallest GRSq or RSq,
0
, meaning no jitter.
A typical useful value is 0.01.discrete.plot.cols(length(x))
which is vector
of distinguishable colors, the first three of which are also distinguishable on
a monochrome printer.
You can examine the colors 0
, meaning no RSq plot. Default is 0
, meaning no "number of predictors" plot.
The special value NULL
means borrow col.grsq
(or col.rsq
NULL
(default) means borrow col.grsq
(or col.rsq
if col.grsq
is NULL
).col.vline
.
Default is "grey"
.3
.1
, meaning draw a legend.
Use 0
for no legend.
The legend is drawn in the cumulative distribution graph,
if that graph is plotted.
Else the legend is drawn in the model comparison chart.NULL
(default) means position the legend automatically.
Else specify c(x,y)
in user coordinates.NULL
(default) means generate the legend text automatically
from call$formula
.
The following settings are related to par()
and are included so yopar()
for global settings as appropriate.
Default is TRUE
,
which sets mfrow, mar=c(4,4,2,3), mgp=c(1.6,0.6,0), cex=0.7
.
Set to FALSE
if you want to append figures to an existing plot.NULL
, meaning generate figure headings automatically.earth
,
plot.earth
,
plot.earth.models
,
plotd
,
plotmo
data(ozone1)
a1 <- earth(O3 ~ ., data = ozone1, degree = 2)
a2 <- earth(O3 ~ .-wind, data = ozone1, degree = 2, nk = 31)
a3 <- earth(O3 ~ .-humidity, data = ozone1, degree = 2, nk = 31)
plot.earth.models(list(a1,a2,a3), rlim=c(.6,.8))
Run the code above in your browser using DataLab