From William.Venables@cmis.CSIRO.AU  Mon Jan 24 22:46:06 2000
From: William.Venables@cmis.CSIRO.AU (Bill Venables)
Date: Tue, 25 Jan 2000 08:46:06 +1000
Subject: tabulate
In-Reply-To: Your message of "23 Jan 2000 13:22:11 +0100."
 <x2ln5hq8p8.fsf@blueberry.kubism.ku.dk>
Message-ID: <200001242246.JAA13516@snowy.nsw.cmis.CSIRO.AU>

> Bill Venables <William.Venables@cmis.CSIRO.AU> writes:
> 
> > table() should be the usual function for everyday use (and using
> > this makes the change to tabulate() unnecessary), but on R the
> > difference in speed between table() and tabulate() can be
> > dramatic.  I would suggest improving the efficiency of table()
> > should be a fairly high priority.
> 
> Would you happen to have a simple example showing the inefficiency? I
> just had a glimpse at the code and it would seem that there are some
> extraneous coercions from integer to double precision and back, which
> could be fairly easily removed.

OK Peter.  This is the first one I cooked up:

Version 0.90.1  (December 15, 1999)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type    "?license" or "?licence" for distribution details.

R is a collaborative project with many contributors.
Type    "?contributors" for a list.

Type    "demo()" for some demos, "help()" for on-line help, or
        "help.start()" for a HTML browser interface to help.
Type    "q()" to quit R.

[Previously saved workspace restored]

> m <- rpois(100000, 1)
> tabulate(m)
[1] 36891 18399  6064  1519   309    50     4     1
> table(m)
m
    0     1     2     3     4     5     6     7     8 
36763 36891 18399  6064  1519   309    50     4     1 
> system.time(tabulate(m))
[1] 0.11 0.00 0.00 0.00 0.00
> system.time(table(m))
[1] 2.90 0.16 4.00 0.00 0.00
> version
         _                   
platform sparc-sun-solaris2.7
arch     sparc               
os       solaris2.7          
system   sparc, solaris2.7   
status                       
major    0                   
minor    90.1                
year     1999                
month    December            
day      15                  
language R                   
> 

-- 
-----------------------------------------------------------------
Bill Venables, Statistician, CMIS Environmetrics Project.

Physical address:                            Postal address:
CSIRO Marine Laboratories,                   PO Box 120,       
233 Middle St, Cleveland, Qld, 4163          Cleveland, Qld, 4163
AUSTRALIA                                    AUSTRALIA

Tel: +61 7 3826 7251           Email: Bill.Venables@cmis.csiro.au     
Fax: +61 7 3826 7304      http://www.cmis.csiro.au/bill.venables/


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._


tabulate
Bill Venables