<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi,<br>
I am trying to do kriging with in strata. I use the example in book
"Applied Spatial Data Analysis with R" page:217. It works fine with
"meuse" data. But when I am trying to use my data it shows some
error massage.<br>
Do I missing something here?<br>
Thanks in advance<br>
Zia<br>
<br>
Note: I am using following codes:<br>
<br>
<br>
<font face="Courier New, Courier, monospace">coordinates(tala) <-
~ x + y<br>
gridded(tala.grid) = ~x+y<br>
</font><font face="Courier New, Courier, monospace">>
names(tala.grid)<br>
<small>[1] "lt.a" "lt.b" "elev" "FDID"
"wfe_was" "was" "was.lt" "elev.sqrt"<br>
</small></font><font face="Courier New, Courier, monospace"><small>#
lt.a= flooding land type 1 and lt.b=</small></font><small><font
face="Courier New, Courier, monospace">flooding </font></small><font
face="Courier New, Courier, monospace"><small>land type 2<br>
<br>
</small></font><font face="Courier New, Courier, monospace">tala$lt.a
<- tala.grid$lt.a[overlay(tala.grid, tala)]<br>
> names(tala)<br>
<small><br>
[1] "YID" "FDID" "elev" "was" "sas" "lt.a" "lt.b" </small></font><font
face="Courier New, Courier, monospace"><small></small><br>
</font><br>
<br>
<font face="Courier New, Courier, monospace">x1 <- krige(sas ~ 1,
tala[tala$lt.a == 0, ],tala.grid[tala.grid$lt.a == 0, ], model =
vgm(35,"Sph", 2000, 15), nmin = 8, nmax = 40, maxdist = 1000)<br>
</font><br>
<font color="#ff0000"><small><font face="Courier New, Courier,
monospace">Error in tala[tala$lt.a == 0, ] : NAs not permitted
in row index Error in krige(sas ~ 1, tala[tala$lt.a == 0, ],
tala.grid[tala.grid$lt.a == : <br>
error in evaluating the argument 'locations' in selecting a
method for function 'krige'</font><br>
</small></font><br>
<font face="Courier New, Courier, monospace">x2<- krige(sas~1,
tala[tala$lt.a == 1, ],tala.grid[tala.grid$lt.a == 1, ], model =
vgm(30,"Sph", 2000), nmin = 8, nmax = 40, maxdist = 1000)</font><br>
<font color="#ff6666"><br>
</font><font color="#ff6666" face="Courier New, Courier, monospace"><small>Error
in tala[tala$lt.a == 1, ] : NAs not permitted in row index Error
in krige(sas ~ 1, tala[tala$lt.a == 1, ],
tala.grid[tala.grid$lt.a == : <br>
error in evaluating the argument 'locations' in selecting a
method for function 'krige'</small></font><font color="#ff6666"><br>
</font><br>
<br>
On 1/22/2011 4:51 PM, Edzer Pebesma wrote:
<blockquote cite="mid:4D3B517F.2050302@uni-muenster.de" type="cite">
<pre wrap="">
On 01/22/2011 11:18 AM, Patrick Giraudoux wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi,
Does anybody knows if it is possible to fix user defined maximum and
minimum values in the key of a xyplot (or spplot) graphics ? I have a
set of seperate drawings and I would like to keep the same color scheme
between fixed values for all of them (which looks like being not
possible, autokey normally adapting the range of colors to the
particular range of values in each plot.
</pre>
</blockquote>
<pre wrap="">
If your object is of class SpatialPointsDataFrame, it may suffice to
pass the argument "cuts" with the vector of cuts. If it is
grid/polygons/lines, use the "at" component of the "colorkey" argument
to specify breaks; see ?levelplot in package lattice.
</pre>
<blockquote type="cite">
<pre wrap="">
Patrick
_______________________________________________
R-sig-Geo mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-sig-Geo@r-project.org">R-sig-Geo@r-project.org</a>
<a class="moz-txt-link-freetext" href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a>
</pre>
</blockquote>
<pre wrap="">
</pre>
</blockquote>
</body>
</html>