[R-sig-Geo] readAsciiGrid (maptools) 2

[Roger Bivand]

On Tue, 13 Nov 2007, Ingo Holz wrote:

> Hi again,
>
> OK, it is not a problem of readAsciiGrid (maptools)!
>
> However:
>
> x <- 3539000.5 : 3540000.5   ## OK, readAsciiGrid adds .5 to the coordinates
> length(x)
> x[1:10]
> min(x)                       ## OK, this is the rounded value
> summary(x)                   ## I still do not understand this
> summary(1:1001)
>
> I get the following results:
>
>> x <- 3539000.5 : 3540000.5   ## OK, readAsciiGrid adds .5 to the coordinates
>> length(x)
> [1] 1001
>> x[1:10]
> [1] 3539001 3539002 3539003 3539004 3539005 3539006 3539007 3539008 3539009
> 3539010
>> min(x)                       ## OK, this is the rounded value
> [1] 3539001
>> summary(x)                   ## I still do not understand this
>   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
> 3539000 3539000 3540000 3540000 3540000 3540000

> x <- seq(3539000.5, 3540000.5, 1)
length(x)
[1] 1001
> options(digits=16)
> summary(x)
      Min.   1st Qu.    Median      Mean   3rd Qu.      Max.
3539000.5 3539250.5 3539500.5 3539500.5 3539750.5 3540000.5

OK?

Roger


>> summary(1:1001)
>   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
>      1     251     501     501     751    1001
>>
>
> Ingo
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

-- 
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no