Calculating Scale/dispersion Of Gamma Glm Using Statsmodels

I'm having trouble obtaining the dispersion parameter of simulated data using statsmodels' GLM function. import statsmodels.api as sm import matplotlib.pyplot as plt import scipy.

Solution 1:

As Josef noted in the comments, statsmodels uses a different kind of parameterization.