# NOT RUN {
data(trade)
gravity = feols(log(Euros) ~ log(dist_km) | Destination + Origin, trade)
# Extracting the 'working' number of observations used to compute the pvalues
fitstat(gravity, "G")
# Idem, but when coimputing two-way SEs
fitstat(gravity, "G", se = "standard")
# }
Run the code above in your browser using DataLab