Мой вопрос является частично статистическим и частично кодирующим.
Скажите, что я установил гамма-линку GLM. do c пример
# Instantiate a gamma family model with the default link function.
gamma_model = sm.GLM(y, X, family=sm.families.Gamma())
gamma_results = gamma_model.fit()
Гамма - это двухпараметрическое распределение формы и масштаба.
Насколько я понимаю, как и в случае регрессии OLS, подобранный GLM предсказывает значения y , которые являются оценками среднего значения и формы гамма-распределения для комбинации X входов. Поэтому я могу использовать GLM для оценки значений параметров формы для данного Xi . Но как мне получить шкалу?
Чего я хотел бы добиться, так это сделать прогнозы на Xi и затем произвести выборку полученного распределения в соответствии с do c пример,
# sample gamma distribution
shape, scale = 2., 2.
s = np.random.gamma(shape, scale, 1000)