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
|
> # test.dots.plotmo.R: test dots functions with the plotmo and earth libraries
>
> source("test.prolog.R")
> library(plotmo)
Loading required package: Formula
Loading required package: plotrix
> library(earth)
> data(ozone1)
> options(warn=1) # print warnings as they occur
>
> a <- earth(O3~., data=ozone1, degree=2)
> expect.err(try(plotmo(a, persp.s=99)), "'s' matches both the 'sub' and 'scale' arguments of persp()")
plotmo grid: vh wind humidity temp ibh dpg ibt vis doy
5760 5 64 62 2112.5 24 167.5 120 205.5
Error : 's' matches both the 'sub' and 'scale' arguments of persp()
Got expected error from try(plotmo(a, persp.s = 99))
>
> # Commented out because we now silently drop partial plot args like cex.l
> # expect.err(try(plotmo(a, cex.l=.8, cex.la=.9)), "arguments 'cex.l' and 'cex.la' both match 'cex.lab' in draw.plot_degree1")
> # expect.err(try(plotmo(a, persp.shad=1, persp.sh=2)), "'persp.shad' and 'persp.sh' both match the 'shade' argument of persp()")
>
> options(warn=2) # treat warnings as errors
>
> # Commented out because we now silently drop partial plot args like cex.l
> # expect.err(try(plotmo(a, cex.l=.8)), "\"cex.l\" is not a graphical parameter")
> # expect.err(try(plotmo(a, cex.lxx=.8)), "\"cex.lxx\" is not a graphical parameter")
> # expect.err(try(plotmo(a, cex.labx=.8)), "\"cex.labx\" is not a graphical parameter")
> # expect.err(try(plotmo(a, cex.l=.8, cex.lab=.9)), "\"cex.l\" is not a graphical parameter")
>
> expect.err(try(plotmo(a, nonesuch=.8)), "predict.earth ignored argument 'nonesuch'")
stats::predict(earth.object, NULL, type="response", nonesuch=0.8)
Error : (converted from warning) predict.earth ignored argument 'nonesuch'
Got expected error from try(plotmo(a, nonesuch = 0.8))
> expect.err(try(plotmo(a, lw=2)), "predict.earth ignored argument 'lw'")
stats::predict(earth.object, NULL, type="response", lw=2)
Error : (converted from warning) predict.earth ignored argument 'lw'
Got expected error from try(plotmo(a, lw = 2))
> options(warn=1)
>
> # test main, xlab, ylab, etc. arguments with recycling
> a <- earth(O3~., data=ozone1, degree=2)
> plotmo(a, caption="test main, xlab, ylab, ticktype arguments",
+ main=c("main1", "main2", "main3", "main4"), xlab=c("x1", "x2"),
+ persp.nticks=2, persp.ticktype="d", ylab=c("y1", "y2", "y3"))
plotmo grid: vh wind humidity temp ibh dpg ibt vis doy
5760 5 64 62 2112.5 24 167.5 120 205.5
>
> par(mfrow=c(2,2), mar = c(3,3,3,1), mgp = c(1.5,.5,0), oma=c(0,0,4,0))
> plotmo(a, trace=1, do.par=FALSE, degree1=1, degree2=1, caption="top: standard\nbottom: lwd=2 thresh=.9") # no errors or warnings
stats::predict(earth.object, NULL, type="response")
stats::fitted(object=earth.object)
got model response from model.frame(O3 ~ vh + wind + humidity + temp + ib...,
data=call$data, na.action="na.fail")
plotmo grid: vh wind humidity temp ibh dpg ibt vis doy
5760 5 64 62 2112.5 24 167.5 120 205.5
> plotmo(a, lwd=2, trace=1, thresh=.9, do.par=FALSE, degree1=1, degree2=1) # no errors or warnings
stats::predict(earth.object, NULL, type="response", thresh=0.9)
stats::fitted(object=earth.object)
got model response from model.frame(O3 ~ vh + wind + humidity + temp + ib...,
data=call$data, na.action="na.fail")
plotmo grid: vh wind humidity temp ibh dpg ibt vis doy
5760 5 64 62 2112.5 24 167.5 120 205.5
>
> source("test.epilog.R")
|
