Econometric help

Costtoimport

AnnualGrowth = as.tibble(wb(country="countries_only", indicator="NY.GDP.MKTP.KD.ZG"))

PerCapita = as.tibble(wb(country="countries_only", indicator="NY.GDP.PCAP.PP.CD"))

PerCapita = rename(PerCapita, PerCapita = value)

AnnualGrowth = rename( AnnualGrowth, AnnualGrowth = value)

AnnualGrowth $date=as.numeric(AnnualGrowth $date)

PerCapita $date=as.numeric(PerCapita $date)

D1=

full_join(db,AnnualGrowth, by =c("Country.code"= "iso3c","DBYear"="date"))

D2=

full_join(D1,PerCapita, by =c("Country.code"= "iso3c","DBYear"="date"))

reg1=

lm(AnnualGrowth ~ Costtoimport, data= D2)

summary(reg1)

reg2=

lm(AnnualGrowth ~ Costtoimport+log(PerCapita)+factor(Country.code)+factor(DBYear), data= D2)

summary(reg2)

stargazer(reg1,reg2, type= "html", out = "/Users/sarahaldafiri/Desktop/reg.html", omit = "factor")

stargazer(data=Final%>%select(GDP, Costtoimport),type="text",summary = TRUE)

ggplot(D2, aes(x= Costtoimport, y= AnnualGrowth))+

geom_point(alpha=.3)+

scale_y_continuous(limits = c(-40,40)) +

scale_x_continuous()+

geom_smooth(method = lm)

4.294e-04

4.161e-04

=============================================================================

Statistic N Mean St. Dev. Min Pctl(25) Pctl(75) Max

-----------------------------------------------------------------------------

AnnualGrowth 9,077 3.894 6.391 -64.047 1.414 6.376 149.973

Costtoimport 1,877 2,247.072 2,072.257 368.000 1,042.000 2,547.000 19,881.000

-----------------------------------------------------------------------------