Data

Data is imported form World Bank database

library("WDI")

indicators <- c(
 'EN.ATM.CO2E.PC', ## CO2 emission in tons per/capita
 "EG.ELC.RNWX.ZS", ## Electricity from renovable sources
 'NY.GDP.PCAP.CD', ## GDP per capita
 'EG.ELC.COAL.ZS', ## Electricity produced from coal
 'AG.PRD.LVSK.XD') ## Livstock

d0 <- WDI(#country=c("PL", "RU"), 
          indicator = indicators,
          start=2015, end=2015)
#str(d0)

Cleaning the data

Removing aggregates

##levels(as.factor(d0$iso2c))
aggregates <- c("1A", "1W", "4E", "7E", "8S", 'B8', 
                'S1', 'S2', 'S3', 'S4',
                "V1", "V2", "V3", "V4", "Z4", "Z7")

dc  <- d0 %>% filter ( ! iso2c %in% aggregates ) %>%
select (EN.ATM.CO2E.PC, NY.GDP.PCAP.CD, EG.ELC.RNWX.ZS, 
         EG.ELC.COAL.ZS, AG.PRD.LVSK.XD)
summary(dc)
##  EN.ATM.CO2E.PC    NY.GDP.PCAP.CD     EG.ELC.RNWX.ZS    EG.ELC.COAL.ZS  
##  Min.   : 0.0379   Min.   :   293.5   Min.   : 0.0000   Min.   : 0.000  
##  1st Qu.: 0.7813   1st Qu.:  1972.1   1st Qu.: 0.1181   1st Qu.: 0.000  
##  Median : 2.5960   Median :  5589.0   Median : 2.3893   Median : 3.425  
##  Mean   : 4.0864   Mean   : 15669.6   Mean   : 6.6102   Mean   :18.799  
##  3rd Qu.: 5.7909   3rd Qu.: 18076.6   3rd Qu.: 8.0707   3rd Qu.:34.233  
##  Max.   :32.4706   Max.   :167313.3   Max.   :65.4437   Max.   :97.467  
##  NA's   :27        NA's   :9          NA's   :77        NA's   :77      
##  AG.PRD.LVSK.XD  
##  Min.   : 77.41  
##  1st Qu.: 99.10  
##  Median :100.12  
##  Mean   : 99.85  
##  3rd Qu.:101.24  
##  Max.   :116.29  
##  NA's   :48

Plotting the data

plot(dc)

CO2 vs …

fit <- lm(EN.ATM.CO2E.PC ~ NY.GDP.PCAP.CD + EG.ELC.RNWX.ZS +
     EG.ELC.COAL.ZS + AG.PRD.LVSK.XD,  data=dc)

summary(fit)
## 
## Call:
## lm(formula = EN.ATM.CO2E.PC ~ NY.GDP.PCAP.CD + EG.ELC.RNWX.ZS + 
##     EG.ELC.COAL.ZS + AG.PRD.LVSK.XD, data = dc)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -13.6894  -2.4363  -0.6568   0.6929  17.7696 
## 
## Coefficients:
##                   Estimate  Std. Error t value             Pr(>|t|)    
## (Intercept)     3.28248279 11.19217919   0.293             0.769773    
## NY.GDP.PCAP.CD  0.00019126  0.00001977   9.673 < 0.0000000000000002 ***
## EG.ELC.RNWX.ZS -0.12614905  0.03462508  -3.643             0.000388 ***
## EG.ELC.COAL.ZS  0.02432969  0.01512948   1.608             0.110240    
## AG.PRD.LVSK.XD -0.00549063  0.11255599  -0.049             0.961168    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.168 on 130 degrees of freedom
##   (115 observations deleted due to missingness)
## Multiple R-squared:  0.4242, Adjusted R-squared:  0.4065 
## F-statistic: 23.95 on 4 and 130 DF,  p-value: 0.000000000000007442