[R-sig-Geo] Zerodist

[Matevž Pavlič]

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