| 12256 |
pd |
1 |
|
| 22999 |
ripley |
2 |
R : Copyright 2003, The R Development Core Team
|
| 26717 |
maechler |
3 |
Version 1.9.0 Under development (unstable) (2003-10-21)
|
| 12256 |
pd |
4 |
|
|
|
5 |
R is free software and comes with ABSOLUTELY NO WARRANTY.
|
|
|
6 |
You are welcome to redistribute it under certain conditions.
|
| 25869 |
maechler |
7 |
Type 'license()' or 'licence()' for distribution details.
|
| 12256 |
pd |
8 |
|
|
|
9 |
R is a collaborative project with many contributors.
|
| 25869 |
maechler |
10 |
Type 'contributors()' for more information.
|
| 12256 |
pd |
11 |
|
| 25869 |
maechler |
12 |
Type 'demo()' for some demos, 'help()' for on-line help, or
|
|
|
13 |
'help.start()' for a HTML browser interface to help.
|
|
|
14 |
Type 'q()' to quit R.
|
| 12256 |
pd |
15 |
|
|
|
16 |
> #### Regression tests for GRAPHICS & PLOTS
|
|
|
17 |
>
|
| 23074 |
ripley |
18 |
> postscript("reg-plot.ps", paper="a4", horizontal=TRUE,
|
|
|
19 |
+ encoding ="ISOLatin1.enc")
|
| 23055 |
ripley |
20 |
> ## since we supply the font metrics, the results depend only on
|
| 23074 |
ripley |
21 |
> ## the encoding used: Windows is different from Unix by default.
|
| 12256 |
pd |
22 |
>
|
|
|
23 |
> options(warn = 1)# print them as they occur
|
|
|
24 |
>
|
|
|
25 |
> plot(0) # this should remain constant
|
|
|
26 |
> str(par(c("usr","xaxp","yaxp")))
|
|
|
27 |
List of 3
|
|
|
28 |
$ usr : num [1:4] 0.568 1.432 -1.080 1.080
|
|
|
29 |
$ xaxp: num [1:3] 0.6 1.4 4.0
|
|
|
30 |
$ yaxp: num [1:3] -1 1 4
|
|
|
31 |
>
|
|
|
32 |
>
|
|
|
33 |
>
|
|
|
34 |
> ### Test for centring of chars. All the chars which are plotted should
|
|
|
35 |
> ### be centred, and there should be no warnings about
|
|
|
36 |
> ### font metrics unknown for character `?'
|
|
|
37 |
>
|
|
|
38 |
> par(pty="s")
|
| 22999 |
ripley |
39 |
> plot(c(-1,16), c(-1,16), type="n", xlab="", ylab="", xaxs="i", yaxs="i")
|
| 12256 |
pd |
40 |
> title("Centred chars in default char set (ISO Latin1)")
|
| 22999 |
ripley |
41 |
> grid(17, 17, lty=1)
|
| 23077 |
ripley |
42 |
> known <- c(32:126, 144:152, 154, 155, 157:255)
|
| 22999 |
ripley |
43 |
>
|
|
|
44 |
> for(i in known) {
|
| 12256 |
pd |
45 |
+ x <- i %% 16
|
|
|
46 |
+ y <- i %/% 16
|
|
|
47 |
+ points(x, y, pch=i)
|
|
|
48 |
+ }
|
|
|
49 |
>
|
|
|
50 |
> par(pty="m")
|
|
|
51 |
>
|
| 13792 |
pd |
52 |
> ## PR 816 (label sizes in dotchart)
|
| 12976 |
pd |
53 |
>
|
|
|
54 |
> ### Prior to 1.2.2, the label sizes were unaffected by cex.
|
|
|
55 |
>
|
|
|
56 |
> data(VADeaths)
|
| 13792 |
pd |
57 |
> dotchart(VADeaths, main = "Death Rates in Virginia - 1940", cex = 0.5)
|
|
|
58 |
> dotchart(VADeaths, main = "Death Rates in Virginia - 1940", cex = 1.5)
|
| 12976 |
pd |
59 |
>
|
| 18374 |
ripley |
60 |
> ## killed by 0 prior to 1.4.0 and in 1.4.1:
|
|
|
61 |
> t1 <- ts(0:100)
|
|
|
62 |
> ## only warnings about values <= 0
|
|
|
63 |
> plot(t1, log = "y")
|
|
|
64 |
Warning in xy.coords(x, NULL, log = log) : 1 y value <= 0 omitted from logarithmic plot
|
|
|
65 |
> plot(cbind(t1, 10*t1, t1 - 4), log="y", plot.type = "single")
|
| 24092 |
ripley |
66 |
Warning in xy.coords(x = matrix(rep.int(tx, k), ncol = k), y = x, log = log) :
|
| 18374 |
ripley |
67 |
7 y values <= 0 omitted from logarithmic plot
|
|
|
68 |
> stopifnot(par("usr")[4] > 3) # log10: ylim[2] = 1000
|
|
|
69 |
>
|
| 19597 |
ripley |
70 |
>
|
|
|
71 |
> ## This one needs to be looked at.
|
|
|
72 |
> ## lty = "blank" killed the fill colour too.
|
|
|
73 |
> plot(1:10, type="n")
|
|
|
74 |
> polygon(c(1, 3, 3, 1), c(1, 1, 3, 3), col="yellow", border="red", lty="blank")
|
|
|
75 |
> rect(6, 6, 10, 10, col="blue", border="red", lty="blank")
|
|
|
76 |
> ## in 1.5.0 all omit the fill colours.
|
|
|
77 |
> data(trees)
|
|
|
78 |
> with(trees, symbols(Height, Volume, circles=Girth/24, inches=FALSE,
|
|
|
79 |
+ lty="blank", bg="blue"))
|
|
|
80 |
> ## in 1.5.0 ignored the lty.
|
|
|
81 |
>
|
| 19891 |
ripley |
82 |
> ## axis() and par(mgp < 0) {keep this example S+ compatible!}:
|
| 19936 |
maechler |
83 |
> lt <- if(is.R()) "31" else 2
|
| 19891 |
ripley |
84 |
> x <- seq(-2,3, len=1001)
|
|
|
85 |
> op <- par(tck= +0.02, mgp = -c(3,2,0))
|
|
|
86 |
> plot(x, x^2 - 1.2, xaxt = "n", xlab="", type ='l', col = 2,
|
|
|
87 |
+ main = "mgp < 0: all ticks and labels inside `frame'")
|
|
|
88 |
> x <- -2:3
|
|
|
89 |
> lines(x, x^2 - 1.2, type ="h", col = 3, lwd=3)
|
| 19936 |
maechler |
90 |
> axis(1, pos = 0, at=-1:1, lty = lt, col=4)## col & lty work only from R 1.6
|
| 19891 |
ripley |
91 |
> par(op)
|
| 19936 |
maechler |
92 |
> axis(1, pos = 0, at=c(-2,2,3), lty = lt, col=4)
|
| 19891 |
ripley |
93 |
> mtext(side=1,"note the x-ticks on the other side of the bars")
|
| 19597 |
ripley |
94 |
>
|
| 23846 |
maechler |
95 |
> ## plot.table(): explicit xlab and ylab for non-1D
|
| 20988 |
ripley |
96 |
> data(UCBAdmissions)
|
|
|
97 |
> plot(UCBAdmissions)# default x- and y-lab
|
|
|
98 |
> plot(UCBAdmissions, xlab = "x label", ylab = "YY")# wrong in 1.5.1
|
| 23846 |
maechler |
99 |
> ## axis suppression
|
|
|
100 |
> plot(tt <- table(c(rep(0,7), rep(1,4), rep(5, 3))), axes = FALSE)
|
|
|
101 |
> plot(tt, xaxt = "n")
|
|
|
102 |
> ## wrong till (incl.) 1.6.x
|
| 20988 |
ripley |
103 |
>
|
| 25869 |
maechler |
104 |
> ## legend with call
|
|
|
105 |
> lo <- legend(2,2, substitute(hat(theta) == that, list(that= pi)))
|
|
|
106 |
> stopifnot(length(lo$text$x) == 1)
|
|
|
107 |
> ## length() was 3 till 1.7.x
|
|
|
108 |
>
|
| 26717 |
maechler |
109 |
> plot(ecdf(c(1:4,8,12)), ylab = "ECDF", main=NULL)
|
|
|
110 |
> ## ylab didn't work till 1.8.0
|
|
|
111 |
>
|