NPCT

If anyone finds themselves needing an n-dimensional implementation of gammaln in Tom Minka’s Lightspeed, I posted the modified code here.

After wasting an absurd amount of time trying to figure out where the bug in my code was for fitting Dirichlet parameters, it turns out that the bug actually wasn’t in my code, but was in my test data. It appears that Stixbox has a bug in its gamma random variate generator, rgamma that [...]

When I run the rgamma function in Stixbox, I regularly get this error:

??? Error using ==> reshape
To RESHAPE the number of elements must not change.
Error in ==> rgamma at 42

So, I changed the last line of the function from this:

x = reshape(x, n/nn(2), nn(2));

to this:

x = reshape(x(1:n), n/nn(2), nn(2));