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,
plotmodata(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