[R-sig-Geo] Selective Kriging error

[Edzer Pebesma]

Thanks for forwarding my off-list answer back to the list.

Do a

library(gstat)
demo(examples)

and look at the Stratified ordinary kriging example.

meuse[meuse$part.a == 0,]

is passed to the krige function, meaning a subset of the full meuse
dataset. In your case this would be something like

yourdata[yourdata$measurementerror == 0,]

On 12/15/2010 02:59 PM, kapo coulibaly wrote:
> Great! But i meant which package and how?
> 
> Thanks
> 
> On Wed, Dec 15, 2010 at 1:40 AM, Edzer Pebesma <
> edzer.pebesma at uni-muenster.de> wrote:
> 
>> Yes, you can.
>>
>> On 12/14/2010 07:24 PM, kapo coulibaly wrote:
>>> I'm doing some interpolations and i would like to know if it is possible
>>> with any of R's geostatistical packages to selectively pick a few points
>>> where I know measurements are "perfect" while the rest of them exhibit
>>> measurement errors?
>>>
>>>       [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>> --
>> 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.52north.org/geostatistics      e.pebesma at wwu.de
>>
> 

-- 
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.52north.org/geostatistics      e.pebesma at wwu.de