Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2014-09-26 Thread Dirk Eddelbuettel
Christian, On 26 September 2014 at 04:41, Christian Gunning wrote: | The Walker Alias method is now (finally) implemented, along with unit | tests. To my knowledge, the only thing not implemented is R's sample2 | - a special version of sample.int for (!replace && is.null(prob) && n | > 1e+07 &&

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2014-09-26 Thread Christian Gunning
Dear all, The Walker Alias method is now (finally) implemented, along with unit tests. To my knowledge, the only thing not implemented is R's sample2 - a special version of sample.int for (!replace && is.null(prob) && n > 1e+07 && size <= n/2). I also added rmultinom, which behaves how an R API

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-04-12 Thread Christian Gunning
Johnathan, Thanks for formatting this, it looks great. I've added the case of without replacement, and a brief discussion of Walker's alias (that depends on the SVN version of RcppArmadillo, r4296), confirmed that it works with a local install of jekyll, and pushed it back to your github repo. Ca

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-29 Thread Dirk Eddelbuettel
On 29 March 2013 at 18:43, Jonathan Olmsted wrote: | Yup. Progress was waiting on finalizing behavior of sample(). | | Current status is visible here: https://github.com/olmjo/rcpp-gallery/blob/ | gh-pages/src/2013-03-19-using-the-Rcpp-based-sample-implementation.Rmd Awesome :) Will check this

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-29 Thread Jonathan Olmsted
Yup. Progress was waiting on finalizing behavior of sample(). Current status is visible here: https://github.com/olmjo/rcpp-gallery/blob/gh-pages/src/2013-03-19-using-the-Rcpp-based-sample-implementation.Rmd -Jonathan - J.P

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-29 Thread Dirk Eddelbuettel
Christian, Thanks -- applied to RcppArmadillo as rev4296. Any news on the Rcpp Gallery writeup on your nice sample() function? [Hint, Hint] Dirk -- Dirk Eddelbuettel | e...@debian.org | http://dirk.eddelbuettel.com ___ Rcpp-devel mailing list Rc

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-28 Thread Christian Gunning
Dirk, See attached for minor edits related to Walker Alias method. Based on the timing tests from the gallery draft and discussion with Johnathan, we decided it was a bad idea to include code that deviates from R's behavior, and is slower than R to boot. Thus, the C++ sample now throws an inform

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-11 Thread Dirk Eddelbuettel
Christian, Thanks again -- this is now committed as rev4278 in the Rcpp repo. Dirk -- Dirk Eddelbuettel | e...@debian.org | http://dirk.eddelbuettel.com ___ Rcpp-devel mailing list Rcpp-devel@lists.r-forge.r-project.org https://lists.r-forge.r-proj

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-11 Thread Dirk Eddelbuettel
On 12 March 2013 at 10:59, Conrad S wrote: | > #include | > which is not expressive enough. I think we should go for something like boost | > where we have to do | > #include | > but that would require a directory | > (inst)/include/RcppArmadillo/RcppArmadillo/extensions

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-11 Thread Conrad S
> #include > which is not expressive enough. I think we should go for something like boost > where we have to do > #include > but that would require a directory > (inst)/include/RcppArmadillo/RcppArmadillo/extensions/sample.h > > which reads silly. Hm. > Better suggestio

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-11 Thread Dirk Eddelbuettel
On 11 March 2013 at 15:55, Dirk Eddelbuettel wrote: | | On 11 March 2013 at 13:39, Christian Gunning wrote: | | > I don't like add-ons/ all that much as it doesn't show (Rcpp)Armadillo. | | > How about any one of these: | | > | | > inst/include/RcppArmadilloXtra/sample.h | | > ins

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-11 Thread Dirk Eddelbuettel
On 11 March 2013 at 13:39, Christian Gunning wrote: | Dirk (and list for reference), | | On Sun, Mar 10, 2013 at 6:16 PM, Dirk Eddelbuettel wrote: | > | > Christian, Jonathan, | > | > On 8 March 2013 at 20:42, Christian Gunning wrote: | > | With some helpful prodding from Jonathan, I sat down wi

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2013-03-11 Thread Christian Gunning
Dirk (and list for reference), On Sun, Mar 10, 2013 at 6:16 PM, Dirk Eddelbuettel wrote: > > Christian, Jonathan, > > On 8 March 2013 at 20:42, Christian Gunning wrote: > | With some helpful prodding from Jonathan, I sat down with the sample > | code today. I moved sample.h to add-ons as per our

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-12-25 Thread Dirk Eddelbuettel
Perfect timing. Meant to email you yesterday regarding sample(), but that fell off the hectic wagon. We should get back to this. On 25 December 2012 at 01:39, Christian Gunning wrote: | I'm working on writing runit tests for the previously-discussed | sample, and I just hit my competence threshol

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-12-25 Thread Christian Gunning
I'm working on writing runit tests for the previously-discussed sample, and I just hit my competence threshold dealing with inst/include. As discussed, the (overloaded) sample() functions are to be called only from user C++ code and are thus not exposed to R. This means that the function definition

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-11-15 Thread Dirk Eddelbuettel
On 15 November 2012 at 04:00, Christian Gunning wrote: | |     Also, you may need to add RNGScope() if you use R's RNG. | | | | | | | | Yes, I should have mentioned that in my post.  Calling RNGScope() is left | to | | the user, since multiple calls to sample() can con

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-11-15 Thread Christian Gunning
> > | Also, you may need to add RNGScope() if you use R's RNG. > | > | > | > | Yes, I should have mentioned that in my post. Calling RNGScope() is > left to > | the user, since multiple calls to sample() can conceivably happen inside > a > | single RNGScope(). > > With that, maybe best to be i

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-11-14 Thread Dirk Eddelbuettel
On 14 November 2012 at 21:25, Christian Gunning wrote: | | | On Wed, Nov 14, 2012 at 6:14 AM, Dirk Eddelbuettel wrote: | | | Thank you, got both emails!  Will take a closer look -- we were busy with | releasing Rcpp 0.10.0 which will rock :) | | | Sorry, I posted to list from wrong

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-11-14 Thread Christian Gunning
On Wed, Nov 14, 2012 at 6:14 AM, Dirk Eddelbuettel wrote: > > Thank you, got both emails! Will take a closer look -- we were busy with > releasing Rcpp 0.10.0 which will rock :) > Sorry, I posted to list from wrong address; 0.10.0 looks great. I found the copious informational links at the bot

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-11-14 Thread Romain Francois
Le 14/11/12 14:05, Christian Gunning a écrit : Dear all, The attached file is for inclusion in RcppArmadillo/src. It's a templated implementation of R's sample that relies on a few Armadillo functions. It should produce results identical to R, except when R uses Walker's alias method (with rep

Re: [Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-11-14 Thread Dirk Eddelbuettel
Christian, On 14 November 2012 at 06:05, Christian Gunning wrote: | Dear all, | | The attached file is for inclusion in RcppArmadillo/src.  It's a templated | implementation of R's sample that relies on a few Armadillo functions.  It | should produce results identical to R, except when R uses Wa

[Rcpp-devel] Sample function(s) for inclusion in RcppArmadillo

2012-11-14 Thread Christian Gunning
Dear all, The attached file is for inclusion in RcppArmadillo/src. It's a templated implementation of R's sample that relies on a few Armadillo functions. It should produce results identical to R, except when R uses Walker's alias method (with replacement, more than 200 nonzero probabilities giv