[R-sig-Geo] some help

[Edzer Pebesma]

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-request at r-project.org <r-sig-geo-request at r-project.org> wrote:
> 
> From: r-sig-geo-request at r-project.org <r-sig-geo-request at r-project.org>
> Subject: R-sig-Geo Digest, Vol 113, Issue 21
> To: r-sig-geo at r-project.org
> Date: Tuesday, January 22, 2013, 3:00 AM
> 
> Send R-sig-Geo mailing list submissions to
>     r-sig-geo at 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-request at r-project.org
> 
> You can reach the person managing the list at
>     r-sig-geo-owner at 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 <pgalpern at gmail.com>
> To: "r-sig-geo at r-project.org" <r-sig-geo at r-project.org>
> Subject: [R-sig-Geo] Class nb (spdep) from class igraph
> Message-ID: <50FD99F1.3040805 at 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 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