Чтобы вычислить псевдо R2 моей модели GLM, мне нужно нулевое отклонение.
формула: 1- (deviance / null.deviance)
Отклонение модели приводится в выходных данных model.summary (), но как мне вычислить / получить доступ к нулевому отклонению?
Мой пример кода:
import statsmodels.formula.api as smf
import statsmodels.api as sm
# generate random dataframe
df = pd.DataFrame(np.random.randint(0,100,size=(100, 3)), columns=list('ABC'))
# negative binomial regression
model = smf.glm(formula = "c ~ a + b", data=df, family=sm.families.NegativeBinomial()).fit()
model.summary()