R Under development (unstable) (2025-07-15 r88411 ucrt) -- "Unsuffered Consequences" Copyright (C) 2025 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. > options(error = recover) > set.seed(15331) > library(intamap) Loading required package: sp > library(automap) > library(gstat) > library(psgp) Error in library(psgp) : there is no package called 'psgp' recover called non-interactively; frames dumped, use debugger() to view > #loadMeuse() > > sessionInfo() R Under development (unstable) (2025-07-15 r88411 ucrt) Platform: x86_64-w64-mingw32/x64 Running under: Windows Server 2022 x64 (build 20348) Matrix products: default LAPACK version 3.12.1 locale: [1] LC_COLLATE=C LC_CTYPE=German_Germany.utf8 [3] LC_MONETARY=C LC_NUMERIC=C [5] LC_TIME=C time zone: Europe/Berlin tzcode source: internal attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] gstat_2.1-4 automap_1.1-20 intamap_1.5-9 sp_2.2-0 loaded via a namespace (and not attached): [1] gtable_0.3.6 compiler_4.6.0 Rcpp_1.1.0 FNN_1.1.4.1 [5] parallel_4.6.0 scales_1.4.0 lattice_0.22-7 ggplot2_3.5.2 [9] R6_2.6.1 plyr_1.8.9 intervals_0.15.5 classInt_0.4-11 [13] sf_1.0-21 iterators_1.0.14 MASS_7.3-65 tibble_3.3.0 [17] units_0.8-7 stars_0.6-8 evd_2.3-7.1 DBI_1.2.3 [21] pillar_1.11.0 RColorBrewer_1.1-3 rlang_1.1.6 reshape_0.8.10 [25] spacetime_1.3-3 doParallel_1.0.17 MBA_0.1-2 cli_3.6.5 [29] magrittr_2.0.3 class_7.3-23 xts_0.14.1 foreach_1.5.2 [33] grid_4.6.0 mvtnorm_1.3-3 lifecycle_1.0.4 vctrs_0.6.5 [37] KernSmooth_2.23-26 proxy_0.4-27 glue_1.8.0 farver_2.1.2 [41] codetools_0.2-20 zoo_1.8-14 abind_1.4-8 e1071_1.7-16 [45] tools_4.6.0 pkgconfig_2.0.3 > > > crs = CRS("epsg:28992") > data("meuse") > coordinates(meuse) <- ~x+y > proj4string(meuse) <- crs > data("meuse.grid") > coordinates(meuse.grid) <- ~x+y > gridded(meuse.grid) <- TRUE > proj4string(meuse.grid) <- crs > > meuse$value = log(meuse$zinc) > meuse.grid = meuse.grid[sample(1:dim(meuse.grid)[1], 100),] > output = interpolate(meuse, meuse.grid, list(mean=T, variance=T, nsim = 5), methodName = "automap") R 2025-07-16 14:23:13.100668 interpolating 155 observations, 100 prediction locations [1] "rgdal has been retired. \n As a result of this, some of the checks on projections in the \n intamap package have disappeared" Checking object ... OK Warning message: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace > summary(t(output$outputTable), digits = 4) x y mean variance Min. :178500 Min. :329700 Min. :4.812 Min. :0.08335 1st Qu.:179330 1st Qu.:330540 1st Qu.:5.239 1st Qu.:0.12361 Median :180060 Median :331140 Median :5.523 Median :0.14205 Mean :179994 Mean :331338 Mean :5.663 Mean :0.16902 3rd Qu.:180670 3rd Qu.:332150 3rd Qu.:6.203 3rd Qu.:0.19770 Max. :181300 Max. :333380 Max. :6.770 Max. :0.36419 sim1 sim2 sim3 sim4 Min. :4.643 Min. :4.305 Min. :4.361 Min. :4.291 1st Qu.:5.290 1st Qu.:5.004 1st Qu.:5.082 1st Qu.:5.154 Median :5.678 Median :5.615 Median :5.693 Median :5.678 Mean :5.780 Mean :5.637 Mean :5.730 Mean :5.702 3rd Qu.:6.218 3rd Qu.:6.088 3rd Qu.:6.319 3rd Qu.:6.240 Max. :7.204 Max. :7.226 Max. :7.446 Max. :7.358 sim5 Min. :4.182 1st Qu.:5.257 Median :5.688 Mean :5.673 3rd Qu.:6.169 Max. :7.252 > > output = interpolate(meuse, meuse.grid, + optList = list(idpRange = seq(0.1, 2.9, 0.5), nfold = 3), + outputWhat = list(mean=TRUE), methodName = "idw") R 2025-07-16 14:23:14.104517 interpolating 155 observations, 100 prediction locations [1] "rgdal has been retired. \n As a result of this, some of the checks on projections in the \n intamap package have disappeared" Checking object ... OK [1] "best idp value found is 2.6 rmse 0.471806664612396" Warning message: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace > summary(t(output$outputTable), digits = 4) x y mean Min. :178500 Min. :329700 Min. :4.786 1st Qu.:179330 1st Qu.:330540 1st Qu.:5.348 Median :180060 Median :331140 Median :5.601 Mean :179994 Mean :331338 Mean :5.691 3rd Qu.:180670 3rd Qu.:332150 3rd Qu.:6.100 Max. :181300 Max. :333380 Max. :6.746 > > > output = interpolate(meuse, meuse.grid, list(mean=T, variance=T),methodName = "transGaussian") R 2025-07-16 14:23:18.435528 interpolating 155 observations, 100 prediction locations [1] "rgdal has been retired. \n As a result of this, some of the checks on projections in the \n intamap package have disappeared" Checking object ... OK Warning message: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace > summary(t(output$outputTable), digits = 4) x y mean variance Min. :178500 Min. :329700 Min. :4.811 Min. :0.0611 1st Qu.:179330 1st Qu.:330540 1st Qu.:5.262 1st Qu.:0.1076 Median :180060 Median :331140 Median :5.521 Median :0.1297 Mean :179994 Mean :331338 Mean :5.653 Mean :0.1471 3rd Qu.:180670 3rd Qu.:332150 3rd Qu.:6.192 3rd Qu.:0.1763 Max. :181300 Max. :333380 Max. :6.764 Max. :0.3693 > > > set.seed(15331) > data(meuse) > meuse = meuse[sample(dim(meuse)[1],30),] > meuse$value=meuse$zinc > coordinates(meuse) = ~x+y > mgrid = spsample(meuse,10,"regular") > gridded(mgrid) = TRUE > output1 = interpolate(meuse, mgrid, list(mean=T, variance=T, excprob = 1000,quantile = 0.5), + methodName = "copula") R 2025-07-16 14:23:19.64389 interpolating 30 observations, 12 prediction locations Checking object ... OK There were 31 warnings (use warnings() to see them) > > output2 = interpolate(meuse, mgrid, list(mean=T, variance=T, excprob = 1000,quantile = 0.5), + methodName = "copula",optList = list(methodParameters = output1$methodParameters)) R 2025-07-16 14:23:20.216548 interpolating 30 observations, 12 prediction locations Checking object ... OK Warning messages: 1: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace 2: In predLoess(object$y, object$x, newx = if (is.null(newdata)) object$x else if (is.data.frame(newdata)) as.matrix(model.frame(delete.response(terms(object)), : pseudoinverse used at 30 3: In predLoess(object$y, object$x, newx = if (is.null(newdata)) object$x else if (is.data.frame(newdata)) as.matrix(model.frame(delete.response(terms(object)), : neighborhood radius 30 4: In predLoess(object$y, object$x, newx = if (is.null(newdata)) object$x else if (is.data.frame(newdata)) as.matrix(model.frame(delete.response(terms(object)), : reciprocal condition number 0 5: In predLoess(object$y, object$x, newx = if (is.null(newdata)) object$x else if (is.data.frame(newdata)) as.matrix(model.frame(delete.response(terms(object)), : pseudoinverse used at 30 6: In predLoess(object$y, object$x, newx = if (is.null(newdata)) object$x else if (is.data.frame(newdata)) as.matrix(model.frame(delete.response(terms(object)), : neighborhood radius 30 7: In predLoess(object$y, object$x, newx = if (is.null(newdata)) object$x else if (is.data.frame(newdata)) as.matrix(model.frame(delete.response(terms(object)), : reciprocal condition number 0 > > output3 = interpolate(meuse, mgrid, list(mean=T, variance=T, excprob = 1000,quantile = 0.5), + methodName = "automap",optList = list(model = c("Exp", "Sph")), cv = TRUE) R 2025-07-16 14:23:20.458118 interpolating 30 observations, 12 prediction locations Checking object ... OK Warning message: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace > > output4 = interpolate(meuse, mgrid, list(mean=T, variance=T, excprob = 1000,quantile = 0.5), + methodName = "psgp", cv = TRUE) R 2025-07-16 14:23:20.898702 interpolating 30 observations, 12 prediction locations Checking object ... OK Error in estimateParameters.default(krigingObject, ...) : there is no parameter estimation method or default method for objects of class psgp Calls: interpolate -> estimateParameters -> estimateParameters.default In addition: Warning message: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace recover called non-interactively; frames dumped, use debugger() to view > > output5 = interpolate(meuse, mgrid, list(mean=T, variance=T, excprob = 1000,quantile = 0.5), cv = TRUE, methodName = "automap") R 2025-07-16 14:23:21.033713 interpolating 30 observations, 12 prediction locations Checking object ... OK Warning message: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace > output6 = interpolate(meuse, mgrid, list(mean=T, variance=T, excprob = 1000,quantile = 0.5), optList = list(variogramModel = output5$variogramModel), + cv = TRUE) R 2025-07-16 14:23:21.53356 interpolating 30 observations, 12 prediction locations Checking object ... OK Warning messages: 1: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace 2: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : Could not find time model for method spatialPredict_automatic 3: In interpolate(meuse, mgrid, list(mean = T, variance = T, excprob = 1000, : was not able to estimate prediction time for methodName automatic > output6$variogramModel$range[2] = 1000 > output7 = interpolate(meuse, mgrid, list(mean=T, variance=T, excprob = 1000,quantile = 0.5), + cv = TRUE, optList = list(variogramModel = output6$variogramModel)) R 2025-07-16 14:23:21.883916 interpolating 30 observations, 12 prediction locations Checking object ... OK Warning messages: 1: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace 2: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : Could not find time model for method spatialPredict_automatic 3: In interpolate(meuse, mgrid, list(mean = T, variance = T, excprob = 1000, : was not able to estimate prediction time for methodName automatic > output8 = interpolate(meuse, mgrid, list(mean=T, variance=T, excprob = 1000,quantile = 0.5), + cv = TRUE, optList = list(nclus = 4), methodName = "automap") R 2025-07-16 14:23:22.232336 interpolating 30 observations, 12 prediction locations Checking object ... OK Warning message: In predictTime(nObs = dim(observations)[1], nPred = nPred, formulaString = formulaString, : using standard model for estimating time. For better platform spesific predictions, please run timeModels <- generateTimeModels() and save the workspace > > all.equal(output5$predictions, output6$predictions) # Should be the same [1] TRUE > all.equal(output5$predictions, output8$predictions) # Should be the same [1] TRUE > all.equal(output5$predictions, output7$predictions) # Should be different [1] "Attributes: < Component \"data\": Component \"var1.pred\": Mean relative difference: 0.09743593 >" [2] "Attributes: < Component \"data\": Component \"var1.var\": Mean relative difference: 0.261375 >" [3] "Attributes: < Component \"data\": Component \"residual\": Mean relative difference: 0.228149 >" [4] "Attributes: < Component \"data\": Component \"zscore\": Mean relative difference: 0.279825 >" > > > summary(t(output$outputTable), digits = 4) x y mean variance Min. :178500 Min. :329700 Min. :4.811 Min. :0.0611 1st Qu.:179330 1st Qu.:330540 1st Qu.:5.262 1st Qu.:0.1076 Median :180060 Median :331140 Median :5.521 Median :0.1297 Mean :179994 Mean :331338 Mean :5.653 Mean :0.1471 3rd Qu.:180670 3rd Qu.:332150 3rd Qu.:6.192 3rd Qu.:0.1763 Max. :181300 Max. :333380 Max. :6.764 Max. :0.3693 > > output2$outputTable - output1$outputTable [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] x 0 0 0 0 0 0 0 0 0 0 0 0 y 0 0 0 0 0 0 0 0 0 0 0 0 mean 0 0 0 0 0 0 0 0 0 0 0 0 variance 0 0 0 0 0 0 0 0 0 0 0 0 excprob1000 0 0 0 0 0 0 0 0 0 0 0 0 quantile0.5 0 0 0 0 0 0 0 0 0 0 0 0 attr(,"transposed") [1] TRUE > > summary(output3$predictions, digits = 4) Object of class SpatialPointsDataFrame Coordinates: min max x 178912 181352 y 329714 333611 Is projected: NA proj4string : [NA] Number of points: 30 Data attributes: var1.pred var1.var observed residual Min. : 102.9 Min. : 60071 Min. : 113.0 Min. :-554.140 1st Qu.: 254.9 1st Qu.: 72990 1st Qu.: 195.5 1st Qu.:-127.947 Median : 451.1 Median : 88288 Median : 311.0 Median : 1.718 Mean : 501.2 Mean : 99497 Mean : 510.3 Mean : 9.058 3rd Qu.: 665.4 3rd Qu.:105580 3rd Qu.: 692.0 3rd Qu.: 87.583 Max. :1257.1 Max. :326038 Max. :1571.0 Max. : 618.833 zscore Min. :-1.650533 1st Qu.:-0.412824 Median : 0.006141 Mean : 0.009715 3rd Qu.: 0.287125 Max. : 1.585608 > summary(output4$predictions, digits = 4) Error in h(simpleError(msg, call)) : error in evaluating the argument 'object' in selecting a method for function 'summary': object 'output4' not found Calls: summary -> .handleSimpleError -> h recover called non-interactively; frames dumped, use debugger() to view > > > proc.time() user system elapsed 11.51 0.73 12.23