[R] generating a sequence

Prof Brian Ripley ripley at stats.ox.ac.uk
Tue Dec 21 16:23:49 CET 1999

OK, so we have now seen essentially three solutions:


A <- matrix(1, n-1, n-1)
rA <- row(A)
rA[rA + col(A) <= n]

Jonathan Rougier:

rA <- outer(1:n, 1:n, "-")

Andy Royle:


Timings for n=500 on a 1994 Sun Sparc 20:

3.29, 4.68 and 23.86 seconds.

I deliberately did not use outer, as it calls apply and 
is normally quite slow.  row and col are internal operations and fast.

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 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

r-help 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-help-request at stat.math.ethz.ch

More information about the R-help mailing list