Re: [R-sig-Geo] Univeral kriging

2009-10-28 Thread TANKISO THEJANE
.   Kabeli --- On Mon, 10/26/09, Edzer Pebesma wrote: From: Edzer Pebesma Subject: Re: [R-sig-Geo] Univeral kriging To: "TANKISO THEJANE" Cc: R-sig-Geo@stat.math.ethz.ch Date: Monday, October 26, 2009, 1:34 PM Tankiso, Here is some tips: - make your code such that cut and paste works f

Re: [R-sig-Geo] Univeral kriging

2009-10-26 Thread Edzer Pebesma
Tankiso, Here is some tips: - make your code such that cut and paste works for us; for me, the ^2 below is vanished (and interpreted as superscript right away, I don't know why) - you didn't test your code, as the function inv does not exist - your code doesn't compute the kriging variance - the c

Re: [R-sig-Geo] Univeral kriging

2009-10-26 Thread Ashton Shortridge
Hi there, good routines already exist in several R libraries for doing this. They will save you trouble in developing bug-free and efficient code. Check out gstat and maybe fields. If you still want to implement kriging manually, these scripts, which I developed for a class exercise, might be

Re: [R-sig-Geo] Univeral kriging

2009-10-26 Thread TANKISO THEJANE
R geo-helpers   i am new to r and would like you to help me as i am trying to go through basic steps. I need to compute a predicted value at point (3,3) and its variance given la set of locations and value. i want to make the code generic for any size such that i would be able to estimate univer