Message-ID: <Pine.LNX.4.44.0304281057060.1247-100000@gannet.stats>
Date: 2003-04-28T10:01:44Z
From: Brian Ripley
Subject: sum(..., na.rm=TRUE) oddity
In-Reply-To: <3EACD1EE.8050804@agric.uwa.edu.au>
I can't reproduce this: can you give us the random seed used or a
reproducible example in some other way. (Preferably some actual data, as
the random-number generator default has changed since 1.6.1.)
On Mon, 28 Apr 2003, rohan sadler wrote:
> I get two different results when using sum() and the switch na.rm. The
> result is correct when na.rm=FALSE.
> Linux Redhat 7.3, R version 1.6.1.
>
> I've had no luck searching the mail archives, so I was hoping somebody
> could explain/check this one for me. I will need to apply the function
> to missing data, simple as it is.
>
> Code:
>
> x<-matrix(runif(20,0,5)%/%1,4,5) # random matrix data set 4x5,
> values are integers between 0 and 4
>
> coords<-c(2,3,3,4) # choose any sub-matrix, with
> c(ymin,ymax,xmin,xmax) being row and column indices.
>
> occupied.counter<-function(coords){sum(x[coords[3]:coords[4],coords[1]:coords[2]])}
>
> occupied.counter(coords)
> [1] 7
>
> occupied.counter<-function(coords){sum(x[coords[3]:coords[4],coords[1]:coords[2]],
> na.rm=TRUE)}
>
> occupied.counter(coords)
> [1] 8
--
Brian D. Ripley, ripley at stats.ox.ac.uk
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UK Fax: +44 1865 272595