[Rcpp-devel] Rcpp updates

2013-02-23 Thread Dirk Eddelbuettel
A few things I seem to not have gotten around to announcing here: i) The Rcpp book is indeed coming, we are still hoping for early May. I blogged a few more words at http://dirk.eddelbuettel.com/blog/2013/02/05#rcpp_book_coming_in_may where I also mentioned the dedicated page

Re: [Rcpp-devel] Working with indicators in Armadillo

2013-02-23 Thread Dirk Eddelbuettel
On 23 February 2013 at 21:33, Simon Zehnder wrote: | Dear Rcpp-Devels, | | I try to write some code in C++, where S is a classification vector (arma::uvec) with entries in between 1 and K: Well, that was an incomplete and unreproducible example. At a minimum, try to make them complete and prov

[Rcpp-devel] Working with indicators in Armadillo

2013-02-23 Thread Simon Zehnder
Dear Rcpp-Devels, I try to write some code in C++, where S is a classification vector (arma::uvec) with entries in between 1 and K: 1 arma::mat repS = arma::repmat(S, 1, K); 2 arma::mat compM = arma::ones(S.n_rows, K); 3 arma::rowvec par_post(K); 4 5 for(unsigned int k = 0; k < K; ++k