setwd("C://Users//EamonnM//Google Drive//PhD//Transport//CarSTOCK//Danish CarSTOCK")
list.files()

Data <- read.csv("Elasticities.csv",stringsAsFactors = FALSE)

Reg1 <- lm(log(Data$Activity) ~ log(Data$GDP.PPP)+log(Data$Fuel.Price))
Reg2 <- lm(log(Data$Sales) ~ log(Data$GDP.PPP)+log(Data$Fuel.Price))
Reg3 <- lm(log(Data$Mileage) ~ log(Data$GDP.PPP)+log(Data$Fuel.Price))
Reg4 <- lm(log(Data$Stock) ~ log(Data$GDP.PPP)+log(Data$Fuel.Price))

par(mfrow=c(2,2))
plot(log(Data$GDP.PPP),log(Data$Activity))
plot(log(Data$GDP.PPP),log(Data$Sales))
plot(log(Data$GDP.PPP),log(Data$Mileage))
plot(log(Data$GDP.PPP),log(Data$Stock))




Reg9 <- lm(log(100*Data$Petrol..1400cc) ~ log(Data$GDP.PPP))
Reg10 <- lm(log(100*Data$Petrol.1400...2000) ~ log(Data$GDP.PPP))
Reg11 <- lm(log(100*Data$Petrol..2000) ~ log(Data$GDP.PPP))
Reg12 <- lm(log(100*Data$Diesel..1400cc) ~ log(Data$GDP.PPP))
Reg13 <- lm(log(100*Data$Diesel.1400...2000) ~ log(Data$GDP.PPP))
Reg14 <- lm(log(100*Data$Diesel..2000) ~ log(Data$GDP.PPP))

Reg15 <- lm(log(100*Data$X.1400) ~ log(Data$GDP.PPP))
Reg16 <- lm(log(100*Data$X1400...2000) ~ log(Data$GDP.PPP))
Reg17 <- lm(log(100*Data$X.2000) ~ log(Data$GDP.PPP))
Reg18 <- lm(log(1000000*Data$Electricity) ~ log(Data$GDP.PPP))

summary(Reg1)
summary(Reg2)
summary(Reg3)
summary(Reg4)
summary(Reg5)
summary(Reg6)
summary(Reg7)
summary(Reg8)
summary(Reg9)
summary(Reg10)
summary(Reg11)
summary(Reg12)
summary(Reg13)
summary(Reg14)
summary(Reg15)
summary(Reg16)
summary(Reg17)