library(readxl)
gravity <- read_excel("gravity.xlsx")
gravity$limports<-log(gravity$imports)
gravity$lgdp_exporter<-log(gravity$gdp_exporter)
gravity$lgdp_importer<-log(gravity$gdp_importer)
gravity$ldist<-log(gravity$dist)
gravity2<-subset(gravity,gravity$imports>0)
library(estimatr)
ols<-lm_robust(limports~lgdp_exporter+lgdp_importer+ldist+comlang_off,data = gravity2)
library(modelsummary)
modelsummary(ols)
|
|
(1)
|
|
(Intercept)
|
−33.751
|
|
|
(0.353)
|
|
lgdp_exporter
|
1.226
|
|
|
(0.008)
|
|
lgdp_importer
|
0.951
|
|
|
(0.008)
|
|
ldist
|
−1.374
|
|
|
(0.020)
|
|
comlang_off
|
1.293
|
|
|
(0.050)
|
|
Num.Obs.
|
19978
|
|
R2
|
0.642
|
|
R2 Adj.
|
0.642
|
|
AIC
|
93181.8
|
|
BIC
|
93229.2
|
|
RMSE
|
2.49
|