Уважаемые коллеги!
Мы проводил расчёты гравитационной модели в программных средах R и STATA.
Для расчётов используется типовой пакет glmm в среде R (с параметром family=quasipoisson
)
и ppml в STATA.
Вызов процедуры расчёта в среде R такой
summary(glmm<-glm(formula=exports ~ ln_GDPimporter + ln_GDPexporter + ln_GDPimppc + ln_GDPexppc + ln_Distance + ln_Tariff + ln_ExchangeRate + Contig + Comlang + Colony_CIS + EAEU_CIS + EU_European_Union, family=quasipoisson(link="log"),data=data_pua))
Результаты расчётов в среде R следующие:
На тех же данных мы провели расчёты в среде STATA.
Вызов процедуры расчёта в среде STATA следующий
ppml exports ln_gdpimporter ln_gdpexporter ln_gdpimppc ln_gdpexppc ln_distance ln_tariff ln_exchangerate contig comlang colony_cis eaeu_cis eu_european_union
Результаты расчётов в среде STATA следующие:
Как видно оценки коэффициентов модели (вторая колонка в таблицах результатов) совпадают как минимум до 4-го знака (!)
Однако, остальные результаты (колонки в таблице результатов, начиная с третьей) в выводах обоих пакетов не совпадают.
Как можно объяснить полученные результаты?
В частности, почему совпадают оценки (вторые столбцы таблиц результатов), но не совпадают стандартные отклонения (Std. Err., третьи столбцы таблиц результатов) ?