[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