Thu, 29 Jan 2009 12:52:37 +0000, Pauli Virtanen wrote: [clip] > It appears that scipy.stats.gamma doesn't have a scale parameter. Oops, obviously it has a scale parameter: In Scipy: >>> scipy.stats.gamma.pdf(5, 2, 0, 1.0/5) 1.7359929831205026e-09 In R: > dgamma(5,2,5) [1] 1.735993e-09 So, no bugs present, just different order of arguments.