R Under development (unstable) (2023-10-23 r85389 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> 
> pdf("levelplot.pdf")
> 
> library(lattice)
> 
> data(volcano)
> 
> 
> foo <-
+     data.frame(z = as.vector(volcano),
+                x = rep(log1p(1:87), 61),
+                y = rep(1:61, each = 87))
> 
> system.time(print(levelplot(z ~ x * y, foo, contour = TRUE)))
   user  system elapsed 
   0.10    0.00    0.09 
> 
> system.time(print(levelplot(z ~ x * y, foo, contour = TRUE,
+                             panel = panel.levelplot.raster)))
   user  system elapsed 
   0.03    0.02    0.04 
Warning message:
In (function (x, y, z, subscripts, at = pretty(z), ..., col.regions = regions$col,  :
  'x' values are not equispaced; output may be wrong
> 
> system.time(print(levelplot(z ~ x * y, foo, contour = TRUE,
+                             region.type = "contour")))
   user  system elapsed 
   0.26    0.22    0.48 
> 
> levelplot(z ~ x * y, subset(foo, z > 150), contour = T)
> 
> ## region.type = "contour" does not support missing grid points
> suppressWarnings(plot(levelplot(z ~ x * y, subset(foo, z > 150),
+                                 contour = T,
+                                 main = "Expect error",
+                                 region.type = "contour")))
> 
> 
> levelplot(z ~ x * y, foo, subset = z > 150, contour = T)
> contourplot(z ~ x * y, foo, subset = z > 150, cuts = 10)
> 
> ## manual colorkey labels
> levelplot(z ~ x * y, foo, subset = z > 150,
+           colorkey = list(labels = list(at     = c(145, 150),
+                                         labels = c("below min", 150))))
> ## lattice <= 0.21-5 put the "below min" label at z=150
> 
> ## This should be an error because 'labels' are specified without 'at':
> 
> e <- 
+     try(print(levelplot(z ~ x * y, foo, subset = z > 150,
+                         colorkey = list(labels = list(labels = c("A", "B"))))),
+         silent = TRUE)
> stopifnot(inherits(e, "try-error"))
> 
> ## subset group interaction has problems (does it any longer?):
> 
> cloud(Sepal.Length ~ Petal.Length * Petal.Width, 
+       data = iris, cex = 0.8, groups = Species, subpanel = panel.superpose, 
+       main = "With subset", screen = list(z = 20, x = -70, y = 0), subset = 30:60)
> 
> cloud(Sepal.Length ~ Petal.Length * Petal.Width, 
+       data = iris[30:60, ], cex = 0.8, groups = Species, subpanel = panel.superpose, 
+       main = "Without subset", screen = list(z = 20, x = -70, y = 0))
> 
> 
> levelplot(Sepal.Length ~ Petal.Length * Petal.Width, 
+           data = iris, groups = Species, 
+           subset = 30:60)
> 
> cloud(z ~ x * y, foo)
> cloud(z ~ x * y, foo, subset = z > 150)
> 
> ## long format: NA's clipped at the beginning (??)
> 
> splom(~iris[,1:4], iris, groups = Species)
> splom(~iris[,1:4], iris, groups = Species, subset = 70:130, auto.key = T)
> 
> parallelplot(~iris[,1:4] | Species, iris, subset = 30:130)
> 
> ## NA-s in matrix
> 
> volna <- volcano
> volna[20:40, 20:40] <- NA
> levelplot(volna)
> levelplot(volna, region.type = "contour")
> cloud(volna)
> wireframe(volna)
> 
> dev.off()
null device 
          1 
> 
> 
> 
> 
> # grid.newpage()
> # pushViewport(viewport())
> # grid.lines(x = c(1:5, NA, 7:10) / 10,
> #            y = c(1:5, NA, 7:10) / 10)
> 
> # grid.points(x = c(1:5, NA, 7:10) / 10,
> #             y = c(1:5, 6, 7:10) / 10)
> 
> 
> 
> 
> # grid.points(x = .5, y = .1, pch = 16, gp = gpar(cex  = 1))
> # grid.points(x = .5, y = .2, pch = 16, gp = gpar(cex  = 3))
> 
> # grid.points(x = .5, y = .3, pch = '+', gp = gpar(cex  = 1))
> # grid.points(x = .5, y = .4, pch = '+', gp = gpar(cex  = 3))
> 
> 
> 
> 
> # viewport[ROOT] ->
> 
> #     (viewport[GRIDVP1182]->
> 
> #      (viewport[GRIDVP1185],
> #       viewport[GRIDVP1187],
> #       viewport[GRIDVP1189],
> #       viewport[GRIDVP1190],
> 
> #       viewport[panel.1]->
> 
> #       (viewport[GRIDVP1184]->
>        
> #        (viewport[subpanel.1.1],
> #         viewport[subpanel.1.2], viewport[subpanel.1.3],
> #         viewport[subpanel.2.1], viewport[subpanel.2.2],
> #         viewport[subpanel.2.3], viewport[subpanel.3.1],
> #         viewport[subpanel.3.2], viewport[subpanel.3.3])),
> 
> #       viewport[panel.2]->
> #       (viewport[GRIDVP1186]->
> #        (viewport[subpanel.1.1],
> #         viewport[subpanel.1.2], viewport[subpanel.1.3],
> #         viewport[subpanel.2.1], viewport[subpanel.2.2],
> #         viewport[subpanel.2.3], viewport[subpanel.3.1],
> #         viewport[subpanel.3.2], viewport[subpanel.3.3])),
> 
> #       viewport[panel.3]->
> #       (viewport[GRIDVP1188]->
> #        (viewport[subpanel.1.1],
> #         viewport[subpanel.1.2], viewport[subpanel.1.3],
> #         viewport[subpanel.2.1], viewport[subpanel.2.2],
> #         viewport[subpanel.2.3], viewport[subpanel.3.1],
> #         viewport[subpanel.3.2], viewport[subpanel.3.3]))))
> 
> 
> proc.time()
   user  system elapsed 
   2.00    0.53    2.51