[R-sig-Geo] the "image" question nobody answered some time ago

[Paulo Justiniano Ribeiro Jr]

Dear Luis

Coud you please send an script and data/workspace such that I can try to
reproduce the probblem here?

best
P.J.


Paulo Justiniano Ribeiro Jr
LEG (Laboratório de Estatística e Geoinformação)
Departamento de Estatística
Universidade Federal do Paraná
Caixa Postal 19.081
CEP 81.531-990
Curitiba, PR  -  Brasil
Tel: (+55) 41 3361 3573
Fax: (+55) 41 3361 3141
e-mail: paulojus AT  ufpr  br
http://leg.ufpr.br/~paulojus

On Wed, 1 Nov 2006, Luis Ridao Cruz wrote:

> 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)
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>