Python statsmodels中的GLM伽玛回归

考虑在Python包statsmodel中适合的GLM gamma函数。

这是代码：

import numpy
import statsmodels.api as sm

model = sm.GLM(ytrain, xtrain, family=sm.families.Gamma(link = sm.genmod.families.links.identity)).fit()

print model.summary()

这为我提供了通过伽玛回归获得的拟合模型参数的摘要。我感兴趣的是上述模型的确切pdf $ P（y | X）$。到目前为止，我可以收集的是model.params * x给出了作为训练数据函数的伽马平均值。如何从摘要中推断pdf的形状？