Re: [R-sig-Geo] cokriging question

2008-09-19 Thread Edzer Pebesma
s Taylor; r-sig-geo@stat.math.ethz.ch Subject: Re: [R-sig-Geo] cokriging question Thanks for the quick responses. I've often done global OK or UK that can take ~ 2-3 days to complete. I always assumed that it was b/c the matrices were so large. Looking at task manager indicated that Rgui.exe only

Re: [R-sig-Geo] cokriging question

2008-09-19 Thread Hengl, T.
From: [EMAIL PROTECTED] on behalf of Dave Depew Sent: Fri 9/19/2008 5:00 PM To: Edzer Pebesma Cc: Chris Taylor; r-sig-geo@stat.math.ethz.ch Subject: Re: [R-sig-Geo] cokriging question Thanks for the quick responses. I've often done global OK or UK that can take ~ 2-3 days to complete. I a

Re: [R-sig-Geo] cokriging question

2008-09-19 Thread Dave Depew
Thanks for the quick responses. I've often done global OK or UK that can take ~ 2-3 days to complete. I always assumed that it was b/c the matrices were so large. Looking at task manager indicated that Rgui.exe only consumed ~ 800 Mb of RAM during the processes. I'll try it by passing the maxd

Re: [R-sig-Geo] cokriging question

2008-09-19 Thread Edzer Pebesma
Good morning Dave (late afternoon here), Chris Taylor wrote: Good morning Edzer and Dave, Thanks for bringing up this point. I had a similar issue recently using krige(). Observations at 5800 locations, attempting to krige() predictions at 112,000 locations resulted in the same "memory.c" er

Re: [R-sig-Geo] cokriging question

2008-09-19 Thread Chris Taylor
Good morning Edzer and Dave, Thanks for bringing up this point. I had a similar issue recently using krige(). Observations at 5800 locations, attempting to krige() predictions at 112,000 locations resulted in the same "memory.c" error message. Reducing predicted locations to <<50,000 and red

Re: [R-sig-Geo] cokriging question

2008-09-19 Thread 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/pip

[R-sig-Geo] cokriging question

2008-09-19 Thread Dave Depew
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]