JuliaStats/Distributions.jl — pull request #327 by davidanthoff (2014-12-24 01:12:16+00:00)
Title: Gamma distribution should compute mode for shape=1

I think, at least I don't see any reason why it would not.