[STATSGRASS] sites import problem

Roger Bivand Roger.Bivand at nhh.no
Mon Mar 17 08:28:35 EST 2003


On Mon, 17 Mar 2003, javier garcia - CEBAS wrote:

> Hi;
> I've imported a sites file into R, but get the following warning messages:
> ------------------------------------------------------------------------------------------------------
> > Cotas <- sites.get(G, slist="cotas_5000", debug=FALSE)
> > ls()
> [1] "Cotas"  "G"      "system"
> > summary(Cotas)
>       east            north              var1
>  Min.   :632206   Min.   :4163710   Min.   :   0.00
>  1st Qu.:657634   1st Qu.:4174570   1st Qu.:  65.13
>  Median :666414   Median :4179580   Median : 140.00
>  Mean   :    NA   Mean   :     NA   Mean   : 161.71
>  3rd Qu.:676889   3rd Qu.:4184840   3rd Qu.: 240.00
>  Max.   :688899   Max.   :4198600   Max.   :1116.18
> Warning messages:
> 1: Integer overflow in sum(.); use sum(as.numeric(.))
> 2: Integer overflow in sum(.); use sum(as.numeric(.))
> ---------------------------------------------------------------------------------------------
> I can't understand the consecuences of this warning nor how to avoid it
> Any help?

Which version of R, and the interface package are you using - it always 
helps to know?

If there are many points, then the sums of the coordinates may exceed 
INT_MAX=2147483647 in calculating the means. The warning means that if you 
need the mean of for example Cotas$east, which is represented as an 
integer, you could do summary(as.double(Cotas$east)), which will run.

For integer northings you run out of integers (at 4000km N of the equator)  
after 540 points, for your eastings you run out at about 3600 points.  
Another way to do this is to divide the coordinates by 1000 to represent
values as km.

I'll look at why sites.get() is returning integers, when it should 
probably be returning doubles.

Roger 

> 
> Javier
> 

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




More information about the grass-stats mailing list