Skip to content
Back to formatted view

Raw Message

Message-ID: <AD5CA6183570B54F92AA45CE2619F9B9EB59A8@gi-zrmk.si>
Date: 2010-11-20T11:55:34Z
From: Matevz Pavlic
Subject: Zerodist

Hi Edzer, 

I have this line that i would like to use for kriging :

OK.Res1<-krige(Res1~1, loc=DF, newdata=grd, model=vmf1)

Where should i input d[-zerodist(d)[,1],]?

Thanks, m
-----Original Message-----
From: r-sig-geo-bounces at stat.math.ethz.ch [mailto:r-sig-geo-bounces at stat.math.ethz.ch] On Behalf Of Edzer Pebesma
Sent: Saturday, November 20, 2010 12:39 PM
To: r-sig-geo at stat.math.ethz.ch
Subject: Re: [R-sig-Geo] Zerodist

if your data are in d, pass as data argument to the kriging function.

This arbitrarily removes duplicate points, so they are ignored in the kriging. If you want to do something more sophisticated, zerodist can be used to identify the duplicate pairs.

On 11/20/2010 12:29 PM, Matev? Pavli? wrote:
> Hi all,
> 
>  
> 
> I got the
> 
>> "chfactor.c", line 130: singular matrix in function LDLfactor()
> 
> Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  : 
> 
> LDLfactor<
> 
>  
> 
> Error, probably because there are some point pairs with zero distance in the matrix. My question is, how can i delete these duplicate nodes from the data set?
> 
>  
> 
> Thanks, m
> 
>  
> 
>  
> 
>  
> 
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> 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

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo