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@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: [EMAIL PROTECTED]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
