make Poisson mean field set lmbda to mean
also make resample() initialize the mean field parameters
add Regression class with Gibbs (plus test coverage)
Merge branch 'regression'
tweak Regression .rvs() to take size named argument
Merge branch 'master' into paralleltempering

Conflicts:
	distributions.py