sampling betas using a sparse design-matrix

This function is equivalent in its behaviour to "sample_betas.c".
Instead of a dense vector of doubles it takes an S4-class, namely "dgCMatrix",
and extracts the pointers to indices and data.
Dotproducts and Vector-Vector additions are replaced by iterations over the nonzero
elements for every column in the design matrix