[R-sig-Geo] colocated cokriging with empty neighbourhood causes crash

Tom Gottfried tom.gottfried at tum.de
Fri Jun 24 17:24:29 CEST 2011

Dear list,

while doing local colocated cokriging it seems that, if some points to be predicted have no value of 
the primary variable in the neighbourhood, this causes R to crash with the messages shown below 
(script with modified code from ASDAR book to generate the crash attached).
Should it be expected that NA is predicted at the respective points (as I understood, this is what 
ordinary kriging and cokriging do when the neighbourhood is empty)?


Output at crash-time:

Intrinsic Correlation found. Good.
[using ordinary cokriging]

  *** caught segfault ***
address 0x50, cause 'memory not mapped'

  1: .Call("gstat_predict", as.integer(nrow(as.matrix(new.X))), as.double(as.vector(raw$locations)), 
     as.vector(new.X), as.integer(block.cols), as.vector(block),     as.vector(bl_weights), 
as.integer(nsim), as.integer(BLUE))
  2: predict.gstat(g.cc, meuse.grid)
  3: predict(g.cc, meuse.grid)

Possible actions:
1: abort (with core dump, if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace

Technische Universität München
Department für Pflanzenwissenschaften
Lehrstuhl für Grünlandlehre
Alte Akademie 12
85350 Freising / Germany
Phone: ++49 (0)8161 715324
Fax:   ++49 (0)8161 713243
email: tom.gottfried at wzw.tum.de
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: test_local_colocated.R
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110624/8669077c/attachment.pl>

More information about the R-sig-Geo mailing list