RおよびStataソフトウェアでの重力モデルの計算:係数が同じで標準誤差が異なるのはなぜですか?
RおよびStataソフトウェアで重力モデルの計算を実行しました。 計算にはglmm、R(パラメーター付きfamily = quasipoisson)とppmlStata の標準パッケージを使用しました。 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の結果は次のとおりです。 同じデータについて、ppml手順を使用して、スタタで計算を実行しました。 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 スタタでの計算結果は次のとおりです。 ご覧のとおり、モデル係数(結果表の2列目)は、少なくとも小数点以下4桁目まで同じです。 …