Running examples in ‘DiceView-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: sectionview3d.function
> ### Title: Plot a contour view of a prediction model or function, including
> ###   design points if available.
> ### Aliases: sectionview3d.function sectionview3d,function,function-method
> ###   sectionview3d.matrix sectionview3d,matrix,matrix-method
> ###   sectionview3d.character sectionview3d,character,character-method
> ###   sectionview3d.km sectionview3d,km,km-method sectionview3d.Kriging
> ###   sectionview3d,Kriging,Kriging-method sectionview3d.NuggetKriging
> ###   sectionview3d,NuggetKriging,NuggetKriging-method
> ###   sectionview3d.NoiseKriging
> ###   sectionview3d,NoiseKriging,NoiseKriging-method sectionview3d.glm
> ###   sectionview3d,glm,glm-method sectionview3d.list
> ###   sectionview3d,list,list-method sectionview3d
> 
> ### ** Examples
> 
> x1 <- rnorm(15)
> x2 <- rnorm(15)
> 
> y <- x1 + x2 + rnorm(15)
> model <- lm(y ~ x1 + x2)
> 
> sectionview3d(function(x) sum(x),
+                        Xlim=cbind(range(x1),range(x2)), col='black')
> DiceView:::.plot3d(x1, x2, y)
$xyz.convert
function (x, y = NULL, z = NULL) 
{
    xyz <- xyz.coords(x, y, z)
    if (angle > 2) {
        temp <- xyz$x
        xyz$x <- xyz$y
        xyz$y <- temp
    }
    y <- (xyz$y - y.add)/y.scal
    return(list(x = xyz$x/x.scal + yx.f * y, y = xyz$z/z.scal + 
        yz.f * y))
}
<bytecode: 0x56053c38f578>
<environment: 0x56053f6138f8>

$points3d
function (x, y = NULL, z = NULL, type = "p", ...) 
{
    xyz <- xyz.coords(x, y, z)
    if (angle > 2) {
        temp <- xyz$x
        xyz$x <- xyz$y
        xyz$y <- temp
    }
    y2 <- (xyz$y - y.add)/y.scal
    x <- xyz$x/x.scal + yx.f * y2
    y <- xyz$z/z.scal + yz.f * y2
    mem.par <- par(mar = mar, usr = usr)
    if (type == "h") {
        y2 <- z.min + yz.f * y2
        segments(x, y, x, y2, ...)
        points(x, y, type = "p", ...)
    }
    else points(x, y, type = type, ...)
}
<bytecode: 0x56053c38c7e8>
<environment: 0x56053f6138f8>

$plane3d
function (Intercept, x.coef = NULL, y.coef = NULL, lty = "dashed", 
    lty.box = NULL, draw_lines = TRUE, draw_polygon = FALSE, 
    polygon_args = list(border = NA, col = rgb(0, 0, 0, 0.2)), 
    ...) 
{
    if (!is.atomic(Intercept) && !is.null(coef(Intercept))) {
        Intercept <- coef(Intercept)
        if (!("(Intercept)" %in% names(Intercept))) 
            Intercept <- c(0, Intercept)
    }
    if (is.null(lty.box)) 
        lty.box <- lty
    if (is.null(x.coef) && length(Intercept) == 3) {
        x.coef <- Intercept[if (angle > 2) 
            3
        else 2]
        y.coef <- Intercept[if (angle > 2) 
            2
        else 3]
        Intercept <- Intercept[1]
    }
    mem.par <- par(mar = mar, usr = usr)
    x <- x.min:x.max
    y <- 0:y.max
    ltya <- c(lty.box, rep(lty, length(x) - 2), lty.box)
    x.coef <- x.coef * x.scal
    z1 <- (Intercept + x * x.coef + y.add * y.coef)/z.scal
    z2 <- (Intercept + x * x.coef + (y.max * y.scal + y.add) * 
        y.coef)/z.scal
    if (draw_polygon) 
        do.call("polygon", c(list(c(x.min, x.min + y.max * yx.f, 
            x.max + y.max * yx.f, x.max), c(z1[1], z2[1] + yz.f * 
            y.max, z2[length(z2)] + yz.f * y.max, z1[length(z1)])), 
            polygon_args))
    if (draw_lines) 
        segments(x, z1, x + y.max * yx.f, z2 + yz.f * y.max, 
            lty = ltya, ...)
    ltya <- c(lty.box, rep(lty, length(y) - 2), lty.box)
    y.coef <- (y * y.scal + y.add) * y.coef
    z1 <- (Intercept + x.min * x.coef + y.coef)/z.scal
    z2 <- (Intercept + x.max * x.coef + y.coef)/z.scal
    if (draw_lines) 
        segments(x.min + y * yx.f, z1 + y * yz.f, x.max + y * 
            yx.f, z2 + y * yz.f, lty = ltya, ...)
}
<bytecode: 0x56053c3852f8>
<environment: 0x56053f6138f8>

$box3d
function (...) 
{
    mem.par <- par(mar = mar, usr = usr)
    lines(c(x.min, x.max), c(z.max, z.max), ...)
    lines(c(0, y.max * yx.f) + x.max, c(0, y.max * yz.f) + z.max, 
        ...)
    lines(c(0, y.max * yx.f) + x.min, c(0, y.max * yz.f) + z.max, 
        ...)
    lines(c(x.max, x.max), c(z.min, z.max), ...)
    lines(c(x.min, x.min), c(z.min, z.max), ...)
    lines(c(x.min, x.max), c(z.min, z.min), ...)
}
<bytecode: 0x56053c2d37f8>
<environment: 0x56053f6138f8>

$contour3d
function (f, x.count = 10, y.count = 10, type = "l", lty = "24", 
    x.resolution = 50, y.resolution = 50, ...) 
{
    if (inherits(f, "lm")) {
        vars <- all.vars(formula(f))
    }
    else vars <- c("z", "x", "y")
    for (x1 in seq(x.range.fix[1], x.range.fix[2], length = x.count)) {
        d <- data.frame(x1, seq(y.range.fix[1], y.range.fix[2], 
            length = y.resolution))
        names(d) <- vars[-1]
        if (inherits(f, "lm")) {
            d[vars[1]] <- predict(f, newdata = d)
        }
        else d[vars[1]] <- f(d[[1]], d[[2]])
        xyz <- xyz.coords(d)
        if (angle > 2) {
            temp <- xyz$x
            xyz$x <- xyz$y
            xyz$y <- temp
        }
        y2 <- (xyz$y - y.add)/y.scal
        x <- xyz$x/x.scal + yx.f * y2
        y <- xyz$z/z.scal + yz.f * y2
        mem.par <- par(mar = mar, usr = usr)
        if (type == "h") {
            y2 <- z.min + yz.f * y2
            segments(x, y, x, y2, ...)
            points(x, y, type = "p", ...)
        }
        else points(x, y, type = type, lty = lty, ...)
    }
    for (x2 in seq(y.range.fix[1], y.range.fix[2], length = y.count)) {
        d <- data.frame(seq(x.range.fix[1], x.range.fix[2], length = x.resolution), 
            x2)
        names(d) <- vars[-1]
        if (inherits(f, "lm")) {
            d[vars[1]] <- predict(f, newdata = d)
        }
        else d[vars[1]] <- f(d[[1]], d[[2]])
        xyz <- xyz.coords(d)
        if (angle > 2) {
            temp <- xyz$x
            xyz$x <- xyz$y
            xyz$y <- temp
        }
        y2 <- (xyz$y - y.add)/y.scal
        x <- xyz$x/x.scal + yx.f * y2
        y <- xyz$z/z.scal + yz.f * y2
        mem.par <- par(mar = mar, usr = usr)
        if (type == "h") {
            y2 <- z.min + yz.f * y2
            segments(x, y, x, y2, ...)
            points(x, y, type = "p", ...)
        }
        else points(x, y, type = type, lty = lty, ...)
    }
}
<bytecode: 0x56053c2cfec0>
<environment: 0x56053f6138f8>

$par.mar
$par.mar$mar
[1] 5.1 3.1 4.1 3.1


> 
> sectionview3d(function(x) {
+                       x = as.data.frame(x)
+                       colnames(x) <- all.vars(model$call)[-1]
+                       p = predict.lm(model, newdata=x, se.fit=TRUE)
+                       list(mean=p$fit, se=p$se.fit)
+                     }, vectorized=TRUE,
+               add=TRUE)
> 
> X = matrix(runif(15*2),ncol=2)
> y = apply(X,1,branin)
> 
> sectionview3d(X, y)
Warning in min(x) : no non-missing arguments to min; returning Inf
Warning in max(x) : no non-missing arguments to max; returning -Inf
Warning in min(y.prty) : no non-missing arguments to min; returning Inf
Warning in max(y.prty) :
  no non-missing arguments to max; returning -Inf
Error in plot.window(c(x1, x2), c(z.min, z.max + yz.f * y.max), asp = asp) : 
  need finite 'xlim' values
Calls: sectionview3d ... sectionview3d.matrix -> .plot3d -> <Anonymous> -> plot.window
Execution halted