It indicates a singular or near-singular covariance matrix. The two most common reasons for this are: (i) you have one or more duplicate observations (try: zerodist(Dembia.utm)) , or (ii) using a Gaussian variogram model.
Deciphering your "script", it seems (ii) might apply. On 01/22/2013 01:09 PM, Hailu wrote: > Dear all,what is LDLfactor in kriging?I got at error message while I was > using OK. the error statement is like this "Error in predict.gstat(g, newdata > = newdata, block = block, nsim = nsim, : LDLfactor " > I was using the following optional script for kriging to run: 1. > <Dembia.OM.krg <- krige(OM_perc~1, loc=Dembia.utm, newdata=Dembia.grid, > model=fit.OM_perc.gau)>2. <Dembia.OM.krg <- krige(OM_perc ~1, Dembia.utm, > Dembia.grid, fit.OM_perc.gau)>3. <Dembia.OM.krg <- krige(OM_perc~1, > newdata=Dembia.utm, block=Dembia.grid, model=fit.OM_perc.gau)> > all of the above give the same error, and kriging does not work.Can you give > me some hint?Thanks > > Hailu > --- On Tue, 1/22/13, r-sig-geo-requ...@r-project.org > <r-sig-geo-requ...@r-project.org> wrote: > > From: r-sig-geo-requ...@r-project.org <r-sig-geo-requ...@r-project.org> > Subject: R-sig-Geo Digest, Vol 113, Issue 21 > To: r-sig-geo@r-project.org > Date: Tuesday, January 22, 2013, 3:00 AM > > Send R-sig-Geo mailing list submissions to > r-sig-geo@r-project.org > > To subscribe or unsubscribe via the World Wide Web, visit > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > or, via email, send a message with subject or body 'help' to > r-sig-geo-requ...@r-project.org > > You can reach the person managing the list at > r-sig-geo-ow...@r-project.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of R-sig-Geo digest..." > > > Today's Topics: > > 1. Class nb (spdep) from class igraph (Paul Galpern) > 2. Re: How to show the color ramp legend with discrete data in > spplot ? (Oscar Perpi?an) > 3. Kriging problem (Alemu Tadesse) > 4. Re: Kriging problem (Alemu Tadesse) > 5. Kriging (Alemu Tadesse) > 6. Re: Class nb (spdep) from class igraph (Roger Bivand) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 21 Jan 2013 13:41:37 -0600 > From: Paul Galpern <pgalp...@gmail.com> > To: "r-sig-geo@r-project.org" <r-sig-geo@r-project.org> > Subject: [R-sig-Geo] Class nb (spdep) from class igraph > Message-ID: <50fd99f1.3040...@gmail.com> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Is there a simple way to produce an nb object (in spdep) from a weighted > or unweighted undirected graph (in igraph). > > Could the output of igraph::graph.adjacency(x) be converted to class nb? > > Alternatively, igraph::graph.edgelist(x) could be manipulated into class > Graph, and then spdep::graph2nb() used. > > Am I missing an obvious way to do this? > > Thanks, > Paul > > > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@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.pebe...@wwu.de _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo