- using log directory '/data/gannet/tmp/ripley/R/packages/tests-Sun/akima.Rcheck'
- using R version 2.11.0 Under development (unstable) (2010-02-03 r51097)
- using session charset: UTF-8
- checking for file 'akima/DESCRIPTION' ... OK
- this is package 'akima' version '0.5-4'
- checking package dependencies ... OK
- checking if this is a source package ... OK
- checking for executable files ... OK
- checking whether package 'akima' can be installed ... OK
- checking package directory ... OK
- checking for portable file names ... OK
- checking for sufficient/correct file permissions ... OK
- checking DESCRIPTION meta-information ... OK
- checking top-level files ... OK
- checking index information ... OK
- checking package subdirectories ... OK
- checking R files for non-ASCII characters ... OK
- checking R files for syntax errors ... OK
- checking whether the package can be loaded ... OK
- checking whether the package can be loaded with stated dependencies ... OK
- checking whether the package can be unloaded cleanly ... OK
- checking for unstated dependencies in R code ... OK
- checking S3 generic/method consistency ... OK
- checking replacement functions ... OK
- checking foreign function calls ... OK
- checking R code for possible problems ... OK
- checking Rd files ... OK
- checking Rd metadata ... OK
- checking Rd cross-references ... OK
- checking for missing documentation entries ... OK
- checking for code/documentation mismatches ... OK
- checking Rd \usage sections ... OK
- checking Rd contents ... OK
- checking data for non-ASCII characters ... OK
- checking line endings in C/C++/Fortran sources/headers ... OK
- checking line endings in Makefiles ... OK
- checking for portable use of $BLAS_LIBS ... OK
- checking examples ... ERROR
Running examples in 'akima-Ex.R' failed.
The error most likely occurred in:
> ### * interp
>
> flush(stderr()); flush(stdout())
>
> ### Name: interp
> ### Title: Gridded Bivariate Interpolation for Irregular Data
> ### Aliases: interp interp.new interp.old
> ### Keywords: dplot
>
> ### ** Examples
>
> data(akima)
> plot(y ~ x, data = akima, main = "akima example data")
> with(akima, text(x, y, formatC(z,dig=2), adj = -0.1))
>
> ## linear interpolation
> akima.li <- interp(akima$x, akima$y, akima$z)
> image (akima.li, add=TRUE)
> contour(akima.li, add=TRUE)
> points (akima, pch = 3)
>
> ## increase smoothness (using finer grid):
> akima.smooth <-
+ with(akima, interp(x, y, z, xo=seq(0,25, length=100),
+ yo=seq(0,20, length=100)))
> image (akima.smooth, main = "interp(<akima data>, *) on finer grid")
> contour(akima.smooth, add = TRUE, col = "thistle")
> points(akima, pch = 3, cex = 2, col = "blue")
> # use triangulation package to show underlying triangulation:
> if(library(tripack, logical.return=TRUE))
+ plot(tri.mesh(akima), add=TRUE, lty="dashed")
>
> # use only 15 points (interpolation only within convex hull!)
> akima.part <- with(akima, interp(x[1:15], y[1:15], z[1:15]))
> image(akima.part)
> title("interp() on subset of only 15 points")
> contour(akima.part, add=TRUE)
> points(akima$x[1:15],akima$y[1:15], col = "blue")
>
> ## spline interpolation, two variants (AMS 526 "Old", AMS 761 "New")
> ## -----------------------------------------------------------------
> ## "Old": use 5 points to calculate derivatives -> many NAs
> akima.sO <- interp.old(akima$x, akima$y, akima$z,
+ xo=seq(0,25, length=100), yo=seq(0,20, length=100), ncp=5)
> table(is.na(akima.sO$z)) ## 3990 NA's; = 40 %
FALSE
10000
> akima.sO <- with(akima,
+ interp.old(x,y,z, xo=seq(0,25, length=100), yo=seq(0,20, len=100), ncp = 4))
> sum(is.na(akima.sO$z)) ## still 3429
[1] 0
> image (akima.sO, main = "interp.old(*, ncp = 4) [almost useless]")
> contour(akima.sO, add = TRUE)
>
> ## "New:"
> akima.spl <- with(akima, interp.new(x,y,z, xo=seq(0,25, length=100),
+ yo=seq(0,20, length=100)))
> ## equivalent call via setting linear=FALSE in interp():
> akima.spl <- with(akima, interp(x,y,z, xo=seq(0,25, length=100),
+ yo=seq(0,20, length=100),
+ linear=FALSE))
>
> contour(akima.spl, main = "smooth interp(*, linear = FALSE)")
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(x) : no non-missing arguments to min; returning Inf
Warning in max(x) : no non-missing arguments to max; returning -Inf
Error in pretty(zlim, nlevels) :
NA/NaN/Inf in foreign function call (arg 1)
Calls: contour -> contour.default -> pretty -> .C
Execution halted
Time 0:30.08, 17.78 + 3.01