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-boun...@stat.math.ethz.ch [mailto:r-sig-geo-boun...@stat.math.ethz.ch] On Behalf Of Edzer Pebesma Sent: Saturday, November 20, 2010 12:39 PM To: r-sig-geo@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@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.pebe...@wwu.de _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
