[R-sig-Geo] cokriging question

[Edzer Pebesma]

Dave,

12000 observations fit, in the c representation, in less than 1 Mb (64 
bytes per observation).

The issue is that you think that passing maxdist to predict.gstat has an 
effect. It doesn't; you need to pass it to function gstat().

The same thing happened in this
https://stat.ethz.ch/pipermail/r-sig-geo/2008-September/004182.html
message, where nmax was passed to predict.gstat, and simulation took 
forever. The other issue in that question was, I suspect, lack of 
standardization of coordinates, used in a trend surface.
--
Edzer

Dave Depew wrote:
> Is there a limit to the # of observations or size of file that can be 
> co-kriged in gstat?
> I have a ~12000 observation data set (2 variables), the variograms, 
> cross variogram and lmc are fit well, and co-kriging starts ok
>
> Linear Model of Coregionalization found. Good.
> [using ordinary cokriging]
>
> then immediately outputs
>
> "memory.c", line 57: can't allocate memory in function m_get()
> Error in predict.gstat(fit.ck, newdata = EcoSAV.grid, maxdist = 100) :
>  m_get
>
> Iv tried different maxdist from 10 to 1000, with exactly the same result.
> I recently upgraded my RAM to 4Gb and flipped the windows XP /3GB switch.
>
>

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Weseler Straße 253, 48151 Münster, Germany.  Phone: +49 251
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de/
http://www.springer.com/978-0-387-78170-9