R Under development (unstable) (2023-10-29 r85433 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.

> library(intamap)
Loading required package: sp
> 
> # set up data:
> data(meuse)
> coordinates(meuse) = ~x+y
> meuse$value = log(meuse$zinc)
> data(meuse.grid)
> gridded(meuse.grid) = ~x+y
> 
> # set up intamap object:
> obj = createIntamapObject(
+ 	observations = meuse,
+ 	predictionLocations = meuse.grid,
+ 	constantBias = 0,
+ 	outputWhat = list(mean = 1, 
+ 		variance = 1, 
+ 		quantile = 0.05, 
+ 		quantile = 0.5, 
+ 		quantile = 0.95,
+ 		excprob = 5.5,
+ 		excprob = 6.6,
+ 		cumdistr = 5.5,
+ 		cumdistr = 6.6,
+ 		cumdistr = 7.9
+ 		)
+ )
[1] "createIntamapObject: formulaString is missing, using:  value~1"
> class(obj) = "linearVariogram"
> 
> # check:
> checkSetup(obj)
Checking object ... OK
> 
> # do interpolation steps:
> obj = preProcess(obj)
> obj = estimateParameters(obj)
> obj = spatialPredict(obj)
> obj = postProcess(obj)
> summary(obj$predictions)
Object of class SpatialPixelsDataFrame
Coordinates:
     min    max
x 178440 181560
y 329600 333760
Is projected: NA 
proj4string : [NA]
Number of points: 3103
Grid attributes:
  cellcentre.offset cellsize cells.dim
x            178460       40        78
y            329620       40       104
Data attributes:
   var1.pred        var1.var      
 Min.   :4.589   Min.   :  2.797  
 1st Qu.:5.186   1st Qu.: 68.236  
 Median :5.536   Median : 95.018  
 Mean   :5.679   Mean   :116.920  
 3rd Qu.:6.184   3rd Qu.:144.255  
 Max.   :7.563   Max.   :460.035  
> output = obj$outputTable
> 
> proc.time()
   user  system elapsed 
   2.50    0.37    2.84