Я сделал отрицательную биномиальную модель, используя R, используя следующий код:
myModel <- glm.nb(V7 ~ V8 + V3 + V10 + V12 + V11 + V5, data = data)
summary(myModel)
V7 является зависимой переменной. V8, V3, V10, V12, V11 и V5 являются независимыми переменными. V8, V3, V10 и V12 являются постоянными параметрами. V11 и V5 являются случайными параметрами. Исходя из их базовых c характеристик, параметры в статистических моделях могут работать как фиксированные параметры или случайные параметры. Является ли параметр случайным или фиксированным, можно определить по полученному стандартному отклонению для каждого параметра. В том случае, если стандартное отклонение плотности параметра является статистически значимым при доверительном интервале 95%, соответствующий параметр является случайным. С другой стороны, если стандартное отклонение параметра не является статистически значимым, параметр считается фиксированным для всего населения. Код для отрицательной биномиальной модели. Код правильный? Мне нужно получить масштабный параметр для избыточной дисперсии и логарифмического правдоподобия при сходимости модели RPNB и логарифмического правдоподобия только с константой.
Я также хочу рассчитать предельный эффект для отрицательных биномиальных моделей со случайными параметрами, используя указанные выше переменные, такие как V8, V3, V10, V12, V11 и V5. Я попытался использовать следующий код, но я получил сообщение об ошибке.
myModel <- glm.nb(V7 ~ V8 + V3 + V10 + V12 + V11 + V5, data = data)
cplot(myModel, what = "effect")
Какой код мне нужно использовать для расчета эластичностей для моделей с отрицательным биномиальным параметром со случайным параметром
Я получил следующий результат R с использованием приведенного выше кода. Как я могу рассчитать или получить t-статистику на основе значения z и df для случайных параметров отрицательной биномиальной модели?
Coefficients: Estimate Std. Error z value Pr(>|z|)
(Intercept) -4.470e-01 4.838e-01 -0.924 0.356
V8 2.377e-05 1.326e-05 1.793 0.073 . V3 1.963e-02 1.958e-02 1.002 0.316
V10 6.144e-03 4.671e-03 1.315 0.188
V12 -1.219e-03 5.726e-03 -0.213 0.831
V11 -6.452e-01 3.924e-01 -1.644 0.100
V5 2.528e-03 6.681e-03 0.378 0.705
Signif. codes: 0 ‘’ 0.001 ‘’ 0.01 ‘’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for Negative Binomial(0.6279) family taken to be 1)
Null deviance: 160.68 on 159 degrees of freedom
Residual deviance: 150.30 on 153 degrees of freedom (126 observations deleted due to missingness) AIC: 491.81