#portable
setwd("C:/Users/harari/Desktop/Dropbox/enseignement/M2 EPOG econometrics/2019-2020/project")

library(tidyverse)



reduced2<-read.csv("./reduced2.csv", header = TRUE)
names(reduced2)

For_first_plm <- reduced2 %>%
select(Company_Name=2, year=3, Kexp=w_ln_capex0,
   Kexp_lag=w_ln_capex0_l,CoreOff = 11, NonCoreNonEnergyOff = 12) %>% 
  filter(!is.na(Kexp_lag)&!is.na(Kexp)) %>% 
  arrange(Company_Name,year)

library(plm)
#small model
model.fe1<-plm(data=For_first_plm,Kexp~ Kexp_lag,
              index=c("Company_Name","year"),model="within")
summary(model.fe1)

For_first_plot<-For_first_plm%>%
  mutate(fitted=model.fe1$model$Kexp-model.fe1$residuals) 

#pooled model plot
ggplot(data=For_first_plot,
       aes(x=Kexp_lag,y=Kexp))+
  geom_point(aes(colour=Company_Name))+
  theme(legend.position = "none")+
  geom_smooth(method="lm")

#fixed effects plot
ggplot(data=For_first_plot)+
  geom_point(aes(x=Kexp_lag, y=Kexp,colour=Company_Name))+
  geom_line(aes(x=Kexp_lag, y=fitted,group=Company_Name,colour=Company_Name))+
  theme(legend.position = "none")

#now a model with an offshoring variable
model.fe2<-plm(data=For_first_plm,Kexp~ Kexp_lag+CoreOff,
               index=c("Company_Name","year"),model="within")
summary(model.fe2)


#more variables
For_plm3 <- reduced2 %>%
  select(Company_Name=2, year=3, Kexp=w_ln_capex0,
         profit_lag=w_ln_profits_l, sales_lag=w_ln_sales_l, longdebt=w_ln_longdebt,
         Kexp_lag=w_ln_capex0_l,CoreOff_lag = 11, NonCoreNonEnergyOff_lag = 12) %>% 
  filter(complete.cases(reduced2)) %>% 
  arrange(Company_Name,year)

model.fe3a<-plm(data=For_plm3,Kexp~ Kexp_lag+ 
                 profit_lag + sales_lag + longdebt ,
               index=c("Company_Name","year"),model="within")
summary(model.fe3a)

model.fe3b<-plm(data=For_plm3,Kexp~ Kexp_lag+ 
                 profit_lag + sales_lag + longdebt +
                 CoreOff_lag + NonCoreNonEnergyOff_lag,
               index=c("Company_Name","year"),model="within")
summary(model.fe3b)