[R-sig-Geo] the "image" question nobody answered some time ago
Luis Ridao Cruz
Luisr at frs.fo
Wed Nov 1 17:04:56 CET 2006
GeoR help,
Below is the code I used to get the kriging
for the object Su[[1]] :
> str(Su[[1]])
List of 2
$ x , y : num [1:1209, 1:2] -4.24 -5.31 -7.57 -7.42
-6.91 ...
..- attr(*, "dimnames")=List of 2
.. ..$ : chr [1:1209] "2235" "2236" "2237" "2238" ...
.. ..$ : chr [1:2] "x" "y"
$ data : num [1:1209] 3.91 3 7.87 3.78 3.43 1.61 1.39 4.22 3.83
4.17 ...
- attr(*, "class")= chr "geodata"
, using the maximum likelihood estimation object "ml" :
> ml
likfit: estimated model parameters:
beta tausq sigmasq phi
"2.901" "1.033" "1.140" "0.429"
likfit: maximised log-likelihood = -1897
########################
luis <- krige.conv(Su[[1]]
, locations = lociLuis1
, krige = krige.control(obj.m = ml))
#######################
But when I try to plot with "image" :
> image(luis, loc = test1, col = grey(seq(1, 0.1, l = 30)))
##################
test <- expand.grid(seq(-6, -5, l = 10), seq(61.5, 62, l = 10))
test1 <- as.matrix(test)
colnames(test1) <- c('x','y')
##################
I get the following:
Warning messages:
1: no non-missing arguments to min; returning Inf
2: no non-missing arguments to max; returning -Inf
What is wrorg?
Thanks in advance
P.S: this message (similar one ) was posted on 31/05/2006 to this list
but no answerwas produced
Windows XP
> version
_
platform i386-pc-mingw32
arch i386
os mingw32
system i386, mingw32
status
major 2
minor 3.0
year 2006
month 04
day 24
svn rev 37909
language R
version.string Version 2.3.0 (2006-04-24)
More information about the R-sig-Geo
mailing list