ENH: add the Edgeworth expansion based on the normal distribution

This PR adds the Edgeworth expansion (normal distribution + cumulants), as discussed in this ML thread:
https://groups.google.com/forum/#!topic/pystatsmodels/so83mPP9nkE

(time flies) No significant changes since that email exchange last year.