考虑在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的形状？