Converted to R6 classes.  Still having trouble sampling from the Dirichlet distribution with the R's.