1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
|
> # test.plotmo.args..R: test dot and other argument handling in plotmo
>
> source("test.prolog.R")
> library(earth)
Loading required package: Formula
Loading required package: plotmo
Loading required package: plotrix
> data(ozone1)
> options(warn=1) # print warnings as they occur
>
> options(warn=2) # treat warnings as errors
> lm.mod <- lm(O3~wind, data=ozone1)
>
> expect.err(try(plotmo(lm.mod, se=2, leve=.95)), "plotmo's 'se' argument is deprecated, please use 'level' instead")
Error : plotmo's 'se' argument is deprecated, please use 'level' instead
Got expected error from try(plotmo(lm.mod, se = 2, leve = 0.95))
>
> expect.err(try(plotmo(lm.mod, se=T)), "plotmo's 'se' argument is deprecated, please use 'level=.95' instead")
Error : (converted from warning) plotmo's 'se' argument is deprecated, please use 'level=.95' instead
Got expected error from try(plotmo(lm.mod, se = T))
>
> expect.err(try(plotmo(lm.mod, se=.8)), "plotmo's 'se' argument is deprecated, please use 'level=.95' instead")
Error : plotmo's 'se' argument is deprecated, please use 'level=.95' instead
Got expected error from try(plotmo(lm.mod, se = 0.8))
>
> expect.err(try(plotmo(lm.mod, level=2)), "level=2 is out of range, try level=.95")
Error : level=2 is out of range, try level=.95
Got expected error from try(plotmo(lm.mod, level = 2))
>
> oz2 <- ozone1[1:40,]
> set.seed(2015)
> a <- earth(O3~temp+wind, dat=oz2, deg=2, nk=21, ncr=3, nfo=3, varmod.me="lm")
>
> expect.err(try(plotmo(a, lw=2, trace=1, thresh=.9, SHOWCALL=TRUE)), "predict.earth ignored argument 'lw'")
stats::predict(earth.object, NULL, type="response", lw=2, thresh=0.9)
Error : (converted from warning) predict.earth ignored argument 'lw'
Got expected error from try(plotmo(a, lw = 2, trace = 1, thresh = 0.9, SHOWCALL = TRUE))
>
> options(warn=1)
>
> # test col.response and friends
> plotmo(a, col.response=2, pch.response=c(1, 2, 20), type2="co", SHOWCALL=TRUE) # pch.response tests back compat
plotmo grid: temp wind
53.5 4
> plotmo(a, pt.col=c(1,2,3), pt.pch=c(1, 2, 20), type2="im", SHOWCALL=TRUE)
plotmo grid: temp wind
53.5 4
> plotmo(a, pt.col=c(1,2,3), pt.pch=paste(1:nrow(oz2)), pt.cex=.8, type2="im", do.par=2, SHOWCALL=TRUE)
plotmo grid: temp wind
53.5 4
> plotmo(a, pt.col=c(1,2,3), pt.pch=paste(1:nrow(oz2)), pt.cex=.8, type2="co", degree1=0, do.par=F)
Warning: plotmo: nothing to plot
> par(org.par)
> plotmo(a, col=2, SHOWCALL=TRUE) # will cause red response points
plotmo grid: temp wind
53.5 4
> plotmo(a, pt.col=4, col=3, persp.col="pink", SHOWCALL=TRUE) # col now goes to lines
plotmo grid: temp wind
53.5 4
>
> # test cex and nrug and smooth
> plotmo(a, cex=.8, SHOWCALL=TRUE, nrug=-1, rug.col=2, rug.lwd=1, smooth.col=3,
+ bty="n", col.lab="darkorange", xlab="an x label", cex.lab=1.2) # esoteric, but they work
plotmo grid: temp wind
53.5 4
> plotmo(a, SHOWCALL=TRUE, density.col=2, density.lty=2, smooth.col=3, smooth.f=.3, col="lightblue")
plotmo grid: temp wind
53.5 4
> plotmo(a, cex=1.2, SHOWCALL=TRUE, nrug="density")
plotmo grid: temp wind
53.5 4
>
> # test caption, grid, interval options
> plotmo(a, caption.col=3, caption.font=2, grid.col="pink",
+ level=.8, SHOWCALL=TRUE)
plotmo grid: temp wind
53.5 4
> plotmo(a, caption.col=2, caption.font=2, caption.cex=.8, grid.col=TRUE, bty="n",
+ level=.8, level.shade="lightblue", level.shade2="red",
+ grid.lty=3, grid.lwd=4, grid.nx=NA, SHOWCALL=TRUE)
plotmo grid: temp wind
53.5 4
>
> # test overall plot args handled by par() and graphics args outside do.par
> par(mfrow=c(2,2), mar = c(3,3,3,1), mgp = c(1.5,.5,0), oma=c(0,0,4,0))
> par(col.main="#456789")
> old.mar <- par("mar")
> old.mfcol <- par("mfcol")
> cat("before par: cex=", par("cex"), " col.main=", par("col.main"),
+ " col.axis=", par("col.axis"), " mar=", par("mar"), " mfcol=", par("mfcol"),
+ "\n", sep="")
before par: cex=0.83 col.main=#456789 col.axis=black mar=3331 mfcol=22
> plotmo(a, mfcol=c(2,3), cex.main=1.4, oma=c(5,5,5,5), SHOWCALL=TRUE)
plotmo grid: temp wind
53.5 4
> plotmo(a, caption="no cex")
plotmo grid: temp wind
53.5 4
> plotmo(a, cex=1, caption="cex=1, plot should be identical to previous page")
plotmo grid: temp wind
53.5 4
> plotmo(a, cex=1.2, caption="cex=1.2")
plotmo grid: temp wind
53.5 4
> plotmo(a, do.par=FALSE, degree2=0, degree1=1, main="do.par=FALSE no cex", caption="a test graphics args with do.par=FALSE")
plotmo grid: temp wind
53.5 4
> plotmo(a, do.par=FALSE, degree2=0, degree1=1, cex=1, main="do.par=FALSE cex=1")
plotmo grid: temp wind
53.5 4
> plotmo(a, do.par=FALSE, degree2=0, degree1=1, cex=.8, main="do.par=FALSE cex=.8")
plotmo grid: temp wind
53.5 4
> plotmo(a, do.par=FALSE, degree2=0, degree1=1, cex=1.1, xlab="xlab", col.main=2, col.axis="blue", col.lab=3, font.lab=2,
+ main="do.par=FALSE cex=1.1, col.main=2\ncol.axis=\"blue\", col.lab=3, font.lab=2")
plotmo grid: temp wind
53.5 4
> plotmo(a, do.par=FALSE, degree1=1, degree2=1, persp.ticktype="d",
+ main="do.par=FALSE persp.ticktype=\"d\"")
Warning: 'degree2' specified but no degree2 plots (maybe use all2=TRUE?)
plotmo grid: temp wind
53.5 4
> # all of these should have been restored
> cat("after par: cex=", par("cex"), " col.main=", par("col.main"),
+ " col.axis=", par("col.axis"), " mar=", par("mar"), " mfcol=", par("mfcol"),
+ "\n", sep="")
after par: cex=0.83 col.main=#456789 col.axis=black mar=3331 mfcol=22
> stopifnot(par("col.main") == "#456789")
> stopifnot(par("mar") == old.mar)
> stopifnot(par("mfcol") == old.mfcol)
> par(col.main=1)
>
> # test aliasing of col with other args, and back compat of col.degree1 vs degree1.col
> data(etitanic)
> a20 <- earth(pclass ~ ., data=etitanic, degree=2)
> plotmo(a20, nresponse=1, col=2, col.degree1=3, persp.col="pink", SHOWCALL=1, degree1=1:2, degree2=1:2)
plotmo grid: survived sex age sibsp parch
0 male 28 0 0
> plotmo(a20, nresponse=1, lty=2, persp.lty=1, SHOWCALL=1, degree1=1:2, degree2=1:2)
plotmo grid: survived sex age sibsp parch
0 male 28 0 0
>
> # test "prednames." with a long predictor name
> data(trees)
> trees.with.long.predname <- trees
> trees.with.long.predname$a_quite_long_variable_name <- trees.with.long.predname$Girth
> trees.with.long.predname$Girth <- NULL
> mod <- earth(Volume~.,data=trees.with.long.predname)
> par(mfrow=c(3,2), mar = c(3,3,3,1), mgp = c(1.5,.5,0), oma=c(0,0,4,0))
> plotmo(mod, do.par=FALSE)
plotmo grid: Height a_quite_long_variable_name
76 12.9
> plotmo(mod, do.par=FALSE, prednames.abbreviate=FALSE)
plotmo grid: Height a_quite_long_variable_name
76 12.9
> expect.err(try(plotmo(mod, do.par=FALSE, prednames.abbreviate=c(1,2))), "the prednames.abbreviate argument is not FALSE, TRUE, 0, or 1")
Error : the prednames.abbreviate argument is not FALSE, TRUE, 0, or 1
Got expected error from try(plotmo(mod, do.par = FALSE, prednames.abbreviate = c(1, 2)))
> plotmo(mod, do.par=FALSE, prednames.minlength=3)
plotmo grid: Height a_quite_long_variable_name
76 12.9
>
> source("test.epilog.R")
|
