I have a rather embarrassing mistake in the last chunk of code. I cut the j>>= 1
line when I was "tightening the code". As Kernighan and Plauger stated, "efficiency often means getting the wrong answer quickly". On Friday, April 4, 2014 10:04:34 PM UTC-5, Stefan Karpinski wrote: > > Yes, this is great. It does make me wonder if we should have a Subset type > that wraps a vector and an integer indicating the values to include in the > subset. Something like this: > > immutable Subset{T} <: AbstractVector{T} > values::Vector{T} > subset::Uint64 > end > > With this type, you could make combinatorial iterators that are very > efficient. > > > On Fri, Apr 4, 2014 at 9:14 PM, John Myles White > <johnmyl...@gmail.com<javascript:> > > wrote: > >> This is really great, Doug. Would you be willing to share the raw Rmd >> source? Fernando Perez wanted to show what your Rpubs document would look >> like in the new native R implementation of IPython. >> >> -- John >> >> On Apr 4, 2014, at 2:17 PM, Douglas Bates <dmb...@gmail.com <javascript:>> >> wrote: >> >> > In describing to statisticians some of the differences in approach >> between R and Julia programming I have used as examples the evaluation of >> the reference distribution for randomization tests comparing means for >> paired designs and for unpaired designs. The documents are >> > >> > http://rpubs.com/dmbates/15250 >> > >> > and >> > >> > http://nbviewer.ipython.org/gist/dmbates/9746197 >> >> >