如何估算包括时变的进口商和出口商固定效应的重力模型？

如何解决如何估算包括时变的进口商和出口商固定效应的重力模型？

我正在尝试使用Poisson拟最大似然（PPML）估计引力模型，但是在增加进出口时变固定影响方面却步履艰难。

我尝试使用以下代码生成虚拟对象：

Gravity_data <- cbind(Gravitydata,dummy(Gravity_data$year,sep = "_"))
Gravity_data <- cbind(Gravitydata,dummy(Gravity_data$iso_o,dummy(Gravity_data$iso_d,sep = "_"))

但是，当我尝试将它们相乘时，例如，

Final <- cbind(Gravity_data,sep = "_")*dummy(Gravity_data$iso_origin,sep = "_")) 我尝试过

Reg_glm <- glm(Tradevolume ~ log(distw) + fta_wto + contig + comlang_off + lgdpo+ lgdpd+ factor(iso_o)*factor(year) + factor(iso_d)*factor(year),family='quasipoisson',data = Final)

或

PPML_model <- ppml( dependent_variable = "Tradevolume",distance = "distw",additional_regressors = c("fta_wto","contig","comlang_off","lgdpo","lgdpd","Both_GATT","year","iso_o*year"," iso_d*year","pairid"),robust = TRUE,data = Final)

或

gravity_pois <- feglm(Tradevolume ~ log(distw)+ fta_wto + contig + comlang_off + lgdpo+ lgdpd+Both_GATT
| factor(iso_o)*factor(year)+ factor(iso_d)*factor(year),family = 'quasipoisson',Final)

所有这些均返回错误。还尝试写iso_d代替factor(iso_d)，但无济于事。

基本上，我正在R中寻找类似于Stata中的内容：

egen exp_time = group(exporter year)
tabulate exp_time,generate(EXPORTER_TIME_FE)
egen imp_time = group(importer year)
tabulate imp_time,generate(IMPORTER_TIME_FE) 

感谢您的帮助。