1 Introduction

1.1 Background

In an increasing energy demanding society, concerns about climate change and energy security are shifting the global energy matrix towards a low-carbon economy. There are some scientists indicate the pace of global warming and climate change is directly proportional to the concentration of carbon dioxide emissions. The motivation to research this topic is inspired by global climate change and CO2 emissions. As the temperature and sea levels remain risen over years in Malaysia, with a surface mean temperature increase of 0.14°C-0.25°C per decade and the annual rise of sea level in coastal Malaysia is approximately 3mm per year. The fast growth in population and economies at a faster pace than in previous generations contributes to more CO2emissions, hastening sea level and temperature rise.

1.2 Research Question

Q1: Do the CO2 emission lead to climate change?

Q2: Does it result from the growth of population and economic activities?

1.3 Non-Technical Explanation

There are 4 variables that might have a correlation with the CO2 emissions. These variables were measured at the categorical levels such as population and GDP; Affected by CO2 emissions is the rising sea level in Malaysia. This data set with interval and ratio variables is fit to use linear regression, by assessing the interaction terms.

There are two regression models that will be developed by including five explanatory variables to understand the emission value of CO2 in Malaysia. One simple model that involved two variables(set Annual CO2 emissions as an independent variable while the surface mean of temperature used to be a dependent variable) is used to explain whether CO2 emissions affect temperature rise, and another multiple-regression model that consists of three variables(Annual CO2 emissions, the number of populations, Annual GDP, in Malaysia) is used to explain whether economic activities and the number of populations affect CO2 emissions in Malaysia.

Linear regression helps us to understand the linear relationship between the dependent and independent variables by creating scatter plots for good visualization to analyze. A simple inear regression is used to show a relationship between two variables. While a multiple linear regression is used to describe data and to explain the relationship between one dependent variable and two or more independent variables. At the center of the multiple linear regression analysis lies the task of fitting a single line through the scatter plot, as can be seen in Model Visualization. If there are any significant outliers, the regression analysis will increase the spread of the estimated coefficient and this reduces the fit of the regression equation. Linear regression eliminates independent variables that do not affect the outcome of interest in the analysis, referring to the value of standard deviation and outliers.

1.4 The below visualized datasets used for our analysis.

1.4.1 Figure 1 Sea Levels of coastal Malaysia 1890-2020

data1<-read.csv("Malaysia_air_pollution.csv")

# Line Chart
ggplot(data1, aes(x = Year, y = Sea_Level, fill = Annual_CO2_emissions)) + 
geom_line(size = 1)

1.4.2 Figure 2 Annual GDP in Malaysia 1890-2020

# Bar Chart
ggplot(data1, aes(x = Year, y = GDP, fill = GDP)) + 
  geom_bar(stat = "identity", position = "dodge")

1.4.3 Figure 3 Annual CO2 emissions in Malaysia 1890-2020

# Point Chart
ggplot(data1, aes(x = Year, y = Annual_CO2_emissions, fill = Annual_CO2_emissions)) +
  geom_point(shape = 21, size = 2)

1.4.4 Figure 4 Population in Malaysia 1890-2020

# Line Chart
ggplot(data1, aes(x = Year, y = Population, fill = Annual_CO2_emissions)) + 
geom_line(size = 1)

1.4.5 Figure 5 Annual Mean of Tamperatures in Malaysia 1890-2020

# Line Chart
ggplot(data1, aes(x = Year, y = Annual_Mean_Tamperature, fill = Annual_CO2_emissions)) + 
geom_line(size = 1)

2 Model Development and Linear regression

2.1 Model 1 Simple linear Regression

Below shows the regression coefficients for the independent variable including:

  • Annual CO2 emissions 1890-2020

2.1.1 Performing Linear Regression

model1 <- data1 %>% 
  select(Annual_CO2_emissions)
glimpse(model1)
## Rows: 130
## Columns: 1
## $ Annual_CO2_emissions <int> 7328, 18320, 47632, 10992, 47632, 109920, 124576,~

2.1.2 Summary of Coefficients

data1<-read.csv("Malaysia_air_pollution.csv")
colnames(data1)[4]<-"Annual_Mean_Tamperature"
# model 1 (CO2 to temperature)
model1<-lm(Annual_Mean_Tamperature~Annual_CO2_emissions,data=data1)

Using lm function.

equatiomatic::extract_eq(model1, use_coefs = TRUE)

\[ \operatorname{\widehat{Annual\_Mean\_Tamperature}} = 25.2 + 0(\operatorname{Annual\_CO2\_emissions}) \]

The simple linear regression equation as shown above.

summary(model1)
## 
## Call:
## lm(formula = Annual_Mean_Tamperature ~ Annual_CO2_emissions, 
##     data = data1)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.47409 -0.11497  0.01166  0.10385  0.76584 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          2.520e+01  2.048e-02 1230.26   <2e-16 ***
## Annual_CO2_emissions 4.543e-09  2.265e-10   20.05   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.2 on 128 degrees of freedom
## Multiple R-squared:  0.7585, Adjusted R-squared:  0.7567 
## F-statistic: 402.1 on 1 and 128 DF,  p-value: < 2.2e-16

Table 1

options(scipen = 999)
tab.model1 <- tidy(model1, conf.int = TRUE)
kable(tab.model1) %>%
  kable_styling()
term estimate std.error statistic p.value conf.low conf.high
(Intercept) 25.19941 0.0204831 1230.25494 0 25.15888 25.23994
Annual_CO2_emissions 0.00000 0.0000000 20.05281 0 0.00000 0.00000
modelsummary(model1)
Model 1
(Intercept) 25.199
(0.020)
Annual_CO2_emissions 0.000
(0.000)
Num.Obs. 130
R2 0.759
R2 Adj. 0.757
AIC 45.6
BIC 37.0
Log.Lik. 25.793
F 402.115
RMSE 0.20
model1 %>%
  tbl_regression()
Characteristic Beta 95% CI1 p-value
Annual_CO2_emissions 0.00 0.00, 0.00 <0.001
1 CI = Confidence Interval

Based on table 1, CO2 emissions have an impact on temperatures with a high estimated value and coefficient. As shown in the linear regression, our conjecture is confirmed by the data that the temperature keeps rising with the increasing CO2 emissions with an accuracy of 76% of the equation, However, The p-value is the evidence against a null hypothesis. In this case, A p-value less than 0.001 will under normal circumstances mean that there is substantial evidence against the null hypothesis.

2.2 Model 2 Multiple linear Regression

Below shows the regression coefficients for the independent variable including:

  • GDP = Annual GDP in Malaysia 1890-2020
  • Population = Malaysia population 1890-2020

2.2.1 Performing Linear Regression

model2 <- data1 %>% 
  select(GDP, Population)
glimpse(model2)
## Rows: 130
## Columns: 2
## $ GDP        <dbl> 844000000, 806488420, 770196441, 745396116, 728624703, 6958~
## $ Population <int> 1533320, 1585304, 1639062, 1694657, 1752135, 1811562, 18730~

2.2.2 Summary of Coefficients

# model2 (GDP, population to CO2)
model2<-lm(Annual_CO2_emissions~GDP+Population,data=data1)
equatiomatic::extract_eq(model2, use_coefs = TRUE)

\[ \operatorname{\widehat{Annual\_CO2\_emissions}} = -6278387.86 + 0(\operatorname{GDP}) + 1.05(\operatorname{Population}) \]

The multiple linear regression equation as shown above.

summary(model2)
## 
## Call:
## lm(formula = Annual_CO2_emissions ~ GDP + Population, data = data1)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -29816986  -4118924    359663   3152055  24121958 
## 
## Coefficients:
##                       Estimate         Std. Error t value             Pr(>|t|)
## (Intercept) -6278387.857179145  1376505.948151349  -4.561          0.000011973
## GDP                0.000344264        0.000009479  36.320 < 0.0000000000000002
## Population         1.050876145        0.201436588   5.217          0.000000733
##                
## (Intercept) ***
## GDP         ***
## Population  ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7509000 on 125 degrees of freedom
##   (因为不存在,2个观察量被删除了)
## Multiple R-squared:  0.9895, Adjusted R-squared:  0.9894 
## F-statistic:  5903 on 2 and 125 DF,  p-value: < 0.00000000000000022

Table 2

tidy(model2)
termestimatestd.errorstatisticp.value
(Intercept)-6.28e+061.38e+06-4.561.2e-05 
GDP0.0003449.48e-0636.3 2.84e-68
Population1.05    0.201   5.227.33e-07
options(scipen = 999)
tab.model2 <- tidy(model2, conf.int = TRUE)
kable(tab.model2) %>%
  kable_styling()
term estimate std.error statistic p.value conf.low conf.high
(Intercept) -6278387.8571791 1376505.9481513 -4.561105 0.0000120 -9002663.9659307 -3554111.748428
GDP 0.0003443 0.0000095 36.320330 0.0000000 0.0003255 0.000363
Population 1.0508761 0.2014366 5.216908 0.0000007 0.6522081 1.449544
modelsummary(model2)
Model 1
(Intercept) 6278387.857
(1376505.948)
GDP 0.000
(0.000)
Population 1.051
(0.201)
Num.Obs. 128
R2 0.990
R2 Adj. 0.989
AIC 4421.1
BIC 4432.5
Log.Lik. 2206.550
F 5902.821
RMSE 7508825.22
model2 %>%
  tbl_regression()
Characteristic Beta 95% CI1 p-value
GDP 0.00 0.00, 0.00 <0.001
Population 1.1 0.65, 1.4 <0.001
1 CI = Confidence Interval

Table 2 indicates that GDP and population have high estimates and coefficients for the effect of CO2, as shown in the linear regression. The generated equation between CO2 emissions and GDP, population with an accuracy of 99%; however, the p-values are still less than 0.001, implying that there is substantial evidence against the original hypothesis.

2.3 Model Assessment

2.3.1 Model 1

fit_A <- augment(model1)
fit_A
Annual_Mean_TamperatureAnnual_CO2_emissions.fitted.resid.hat.sigma.cooksd.std.resid
25.1732825.2-0.0694 0.0105 0.2010.000646-0.349 
25.11832025.2-0.0795 0.0105 0.2010.000847-0.4   
25.14763225.2-0.13   0.0105 0.2  0.00225 -0.652 
25  1099225.2-0.189  0.0105 0.2  0.00481 -0.952 
25.14763225.2-0.0996 0.0105 0.2010.00133 -0.501 
25.110992025.2-0.0799 0.0105 0.2010.000855-0.402 
25  12457625.2-0.16   0.0105 0.2  0.00342 -0.804 
25.19160025.2-0.14   0.0105 0.2  0.00262 -0.703 
25.112824025.2-0.13   0.0105 0.2  0.00226 -0.654 
25.112091225.2-0.14   0.0105 0.2  0.00262 -0.704 
25.35862425.20.14   0.0105 0.2  0.00264 0.705 
25.44763225.20.21   0.0105 0.2  0.00593 1.06  
25.35129625.20.14   0.0105 0.2  0.00264 0.706 
25.55129625.20.29   0.0105 0.1990.0113  1.46  
25.53297625.20.26   0.0105 0.1990.00908 1.31  
25.34030425.20.0604 0.0105 0.2010.0004890.304 
25.13664025.2-0.0796 0.0105 0.2010.000848-0.4   
25  8793625.2-0.18   0.0105 0.2  0.00433 -0.904 
25.49160025.20.24   0.0105 0.2  0.00772 1.21  
25.416854425.20.16   0.0105 0.2  0.00342 0.803 
25.223449625.20.0195 0.0105 0.2015.09e-050.0981
25.31832025.20.0705 0.0105 0.2010.0006660.354 
25.26228825.20.0103 0.0105 0.2011.42e-050.0518
25.214656025.20.0199 0.0105 0.2015.31e-050.1   
25.220518425.2-0.0203 0.0105 0.2015.53e-05-0.102 
25.454960025.20.148  0.0104 0.2  0.00292 0.744 
25.365219225.20.0576 0.0104 0.2010.0004420.29  
25.367051225.20.0675 0.0104 0.2010.0006060.34  
25.176944025.2-0.0729 0.0104 0.2010.000706-0.367 
25.4111019225.20.146  0.0104 0.2  0.0028  0.732 
25.5141430425.20.254  0.0103 0.1990.00852 1.28  
25.2200420825.2-0.0485 0.0103 0.2010.000308-0.244 
25.2258312025.2-0.0111 0.0102 0.2011.62e-05-0.056 
25.1250617625.2-0.131  0.0102 0.2  0.00223 -0.657 
25.3296417625.20.107  0.0101 0.2010.00149 0.538 
25.3340019225.20.0851 0.0101 0.2010.0009340.428 
25.3337820825.20.125  0.0101 0.2  0.00202 0.629 
25.2380689625.20.0333 0.0101 0.2010.0001420.167 
25.2411467225.2-0.008110.01   0.2018.4e-06 -0.0407
25.4367132825.20.184  0.0101 0.2  0.00435 0.924 
25.9270769625.20.678  0.0102 0.1910.0598  3.41  
25.1179536025.2-0.0776 0.0103 0.2010.00079 -0.39  
25.1159017625.2-0.107  0.0103 0.2010.0015  -0.536 
25.1173307225.2-0.0773 0.0103 0.2010.000784-0.388 
25.1180635225.2-0.0776 0.0103 0.2010.000791-0.39  
25.2204817625.2-0.0487 0.0103 0.2010.000311-0.245 
25.2234496025.20.0199 0.0102 0.2015.18e-050.1   
25.3191260825.20.132  0.0103 0.2  0.00228 0.663 
25.3170376025.20.0928 0.0103 0.2010.00113 0.467 
25.3256846425.20.129  0.0102 0.2  0.00216 0.648 
25.3218374425.20.0507 0.0102 0.2010.0003350.255 
25.265952025.20.0276 0.0104 0.2010.0001010.139 
25.3131904025.20.0946 0.0103 0.2010.00118 0.476 
25.1110286425.2-0.0744 0.0104 0.2010.000733-0.374 
25  61188825.2-0.172  0.0104 0.2  0.00394 -0.866 
24.960456025.2-0.332  0.0104 0.1990.0147  -1.67  
25.368516825.20.0675 0.0104 0.2010.0006050.339 
25.4119026225.20.225  0.0104 0.2  0.0067  1.13  
25.2122011225.2-0.015  0.0103 0.2012.96e-05-0.0752
25.2365534825.20.024  0.0101 0.2017.39e-050.121 
25  522036725.2-0.263  0.0099 0.1990.00874 -1.32  
24.9637088625.2-0.378  0.009780.1980.0179  -1.9   
24.9649889925.2-0.319  0.009770.1990.0127  -1.6   
24.9700416425.2-0.301  0.009720.1990.0112  -1.51  
24.8742972625.2-0.423  0.009670.1970.0221  -2.13  
25.1756252925.2-0.134  0.009660.2  0.0022  -0.672 
24.8351333125.2-0.375  0.0101 0.1980.0181  -1.89  
24.9383209925.2-0.337  0.0101 0.1980.0145  -1.69  
25.1318324425.2-0.0939 0.0101 0.2010.00114 -0.472 
25  420145025.2-0.188  0.01   0.2  0.00454 -0.947 
24.9468491925.2-0.341  0.009960.1980.0147  -1.71  
25.2477654925.20.0289 0.009950.2010.0001060.145 
25.3582796825.20.104  0.009840.2010.00136 0.523 
24.9739940925.2-0.283  0.009680.1990.00988 -1.42  
24.9838022925.2-0.287  0.009580.1990.0101  -1.44  
24.8983439325.2-0.474  0.009440.1960.027   -2.38  
25  1012760225.2-0.215  0.009410.2  0.00556 -1.08  
25.21053022625.2-0.0873 0.009370.2010.000909-0.438 
25.3913775825.20.0191 0.0095 0.2014.41e-050.0958
25.21458588925.3-0.0557 0.009020.2010.000356-0.28  
25.41666259425.30.0749 0.008850.2010.0006320.376 
25.41790062525.30.0693 0.008760.2010.0005350.348 
25.61750512725.30.371  0.008790.1980.0154  1.86  
25  1904467925.3-0.326  0.008670.1990.0117  -1.64  
25.21944515825.3-0.118  0.008650.2  0.00153 -0.591 
25.22389764325.3-0.058  0.008360.2010.000357-0.291 
25.62262084025.30.318  0.008440.1990.0108  1.6   
25.42326044425.30.145  0.0084 0.2  0.00224 0.728 
25.22730593725.3-0.0835 0.008180.2010.000724-0.419 
25.62803259725.30.293  0.008140.1990.0089  1.47  
25.63087766125.30.26   0.008010.1990.0069  1.31  
25.53064194025.30.171  0.008020.2  0.00299 0.861 
25.43804830725.40.0477 0.007790.2010.0002250.24  
25.53479443325.40.153  0.007870.2  0.00233 0.766 
25.53632977125.40.176  0.007830.2  0.00307 0.881 
25.44010188525.40.0684 0.007750.2010.0004610.343 
25.84087768425.40.385  0.007740.1980.0146  1.93  
26.24286920425.40.766  0.007710.1890.0574  3.84  
25.55008387025.40.103  0.007710.2010.00104 0.517 
25.75426936425.40.214  0.007770.2  0.00452 1.07  
25.86578513425.50.262  0.008160.1990.0071  1.31  
25.97296101625.50.409  0.008580.1970.0183  2.05  
25.78671863725.60.147  0.009750.2  0.00267 0.737 
25.79063098825.60.0789 0.0102 0.2010.0008070.396 
25.711393366425.70.013  0.0135 0.2012.93e-050.0655
25.611344194225.7-0.0748 0.0134 0.2010.000963-0.376 
25.612208458625.8-0.184  0.015  0.2  0.00654 -0.927 
25.411423750025.7-0.338  0.0135 0.1980.0199  -1.7   
25.710829144025.7-0.0114 0.0126 0.2012.08e-05-0.0572
25.812620304925.80.0273 0.0158 0.2010.0001520.137 
25.613398816525.8-0.228  0.0175 0.2  0.0118  -1.15  
25.813511820825.8-0.0432 0.0177 0.2010.000429-0.218 
25.915652133225.9-0.0405 0.0232 0.2010.000497-0.205 
25.817110200426  -0.197  0.0276 0.2  0.0141  -0.997 
26  18133280026  -0.0132 0.031  0.2017.16e-05-0.0669
26.317878530826  0.288  0.0301 0.1990.0333  1.46  
25.918245536326  -0.178  0.0313 0.2  0.0133  -0.906 
25.920198518126.1-0.227  0.0386 0.2  0.0269  -1.16  
26.220071294426.10.0588 0.0381 0.2010.00178 0.3   
26.121593515426.2-0.13   0.0444 0.2  0.0103  -0.667 
26.221699352726.2-0.0282 0.0449 0.2010.00049 -0.144 
26.221521303426.20.0201 0.0441 0.2010.0002430.103 
26.224248743826.3-0.0637 0.0569 0.2010.00324 -0.328 
26.324483256726.3-0.0341 0.0581 0.2010.000954-0.176 
26.323184394026.30.065  0.0517 0.2010.00304 0.334 
26.425056064226.30.0202 0.061  0.2010.0003520.104 
26.425056132026.30.0604 0.061  0.2010.00315 0.311 
26.427222935326.40.002090.073  0.2014.64e-060.0109
26.527865925526.50.0131 0.0767 0.2010.0001920.0679
26.527260743426.40.0807 0.0732 0.2010.00694 0.419 

2.3.1.1 Visualization

ggplot(data1,aes(x=Annual_CO2_emissions,y=Annual_Mean_Tamperature))+
  geom_point()+
  stat_smooth(method="lm",col="red")+
  xlab("CO2 Emissions")+
  ylab("Temperature")+
  ggtitle("Line model visulizing model1")+
  theme(plot.title = element_text(hjust = 0.5))

2.3.2 Model 2

fit_B <- augment(model2)
fit_B
Annual_CO2_emissionsGDPPopulation.fitted.resid.hat.sigma.cooksd.std.resid
73288.44e+081533320-4.38e+064.38e+060.023  7.53e+060.00274 0.591 
183208.06e+081585304-4.33e+064.35e+060.0227 7.53e+060.00266 0.586 
476327.7e+08 1639062-4.29e+064.34e+060.0224 7.53e+060.0026  0.584 
109927.45e+081694657-4.24e+064.25e+060.022  7.53e+060.00246 0.573 
476327.29e+081752135-4.19e+064.23e+060.0217 7.53e+060.0024  0.57  
1099206.96e+081811562-4.14e+064.25e+060.0214 7.53e+060.00238 0.571 
1245766.73e+081873003-4.08e+064.2e+06 0.021  7.53e+060.00229 0.566 
916006.58e+081936527-4.02e+064.11e+060.0207 7.53e+060.00215 0.553 
1282406.29e+082002203-3.96e+064.09e+060.0203 7.53e+060.00209 0.55  
1209121.83e+092066448-3.48e+063.6e+06 0.0202 7.53e+060.00161 0.484 
586241.83e+092129124-3.41e+063.47e+060.0199 7.53e+060.00147 0.467 
476322.17e+092190086-3.23e+063.28e+060.0196 7.53e+060.0013  0.441 
512962.39e+092249181-3.09e+063.14e+060.0194 7.53e+060.00118 0.423 
512962.12e+092306249-3.13e+063.18e+060.019  7.53e+060.00118 0.427 
329762.03e+092364764-3.09e+063.13e+060.0187 7.53e+060.00113 0.42  
403042.19e+092424761-2.98e+063.02e+060.0185 7.53e+060.00103 0.406 
366402.07e+092486279-2.95e+062.99e+060.0182 7.53e+060.0009960.402 
879362.28e+092549355-2.82e+062.9e+06 0.0179 7.53e+060.0009250.39  
916002.7e+09 2614029-2.6e+06 2.69e+060.0177 7.54e+060.0007850.362 
1685442.59e+092677183-2.57e+062.74e+060.0174 7.53e+060.0007990.368 
2344963.11e+092738718-2.33e+062.57e+060.0172 7.54e+060.0006920.345 
183203.31e+092798533-2.2e+06 2.22e+060.0169 7.54e+060.0005090.298 
622883.97e+092856522-1.91e+061.97e+060.0168 7.54e+060.0003980.265 
1465604.54e+092912576-1.65e+061.8e+06 0.0166 7.54e+060.0003290.242 
2051844.94e+092969729-1.46e+061.66e+060.0164 7.54e+060.0002770.223 
5496004.54e+093028000-1.53e+062.08e+060.0161 7.54e+060.0004270.28  
6521925.73e+093087412-1.06e+061.71e+060.0161 7.54e+060.0002870.23  
6705126.11e+093147988-8.67e+051.54e+060.0159 7.54e+060.0002290.206 
7694406.52e+093212623-6.58e+051.43e+060.0157 7.54e+060.0001950.192 
11101926.33e+093279664-6.52e+051.76e+060.0154 7.54e+060.0002910.236 
14143047.86e+093349182-5.3e+04 1.47e+060.0153 7.54e+060.0002010.197 
20042088.89e+0934212533.78e+051.63e+060.0151 7.54e+060.0002440.218 
25831207.75e+0934959556.5e+04 2.52e+060.0147 7.54e+060.0005690.338 
25061767.9e+09 35733671.96e+052.31e+060.0145 7.54e+060.00047 0.31  
29641766.84e+093652493-8.42e+043.05e+060.0141 7.53e+060.0007970.409 
34001928.8e+09 37333716.74e+052.73e+060.014  7.54e+060.0006340.366 
33782088.46e+0938160406.44e+052.73e+060.0137 7.53e+060.0006240.367 
38068961.07e+1039005401.51e+062.3e+06 0.0137 7.54e+060.0004390.308 
41146721.39e+1039869102.7e+06 1.42e+060.0137 7.54e+060.0001670.19  
36713281.39e+1040742672.79e+068.86e+050.0134 7.54e+066.39e-050.119 
27076961.39e+1041626112.87e+06-1.62e+050.0131 7.54e+062.08e-06-0.0217
17953601.29e+1042519452.62e+06-8.2e+05 0.0128 7.54e+065.22e-05-0.11  
15901761.31e+1043422712.8e+06 -1.21e+060.0126 7.54e+060.000111-0.162 
17330721.34e+1044335903.01e+06-1.28e+060.0123 7.54e+060.000122-0.171 
18063521.31e+1045268293e+06       -1.2e+06 0.0121 7.54e+060.000105-0.161 
20481761.45e+1046220293.57e+06-1.53e+060.0119 7.54e+060.000168-0.204 
23449601.61e+1047192314.23e+06-1.89e+060.0118 7.54e+060.000254-0.253 
19126081.31e+1048184773.29e+06-1.38e+060.0114 7.54e+060.000132-0.185 
17037601.36e+1049198103.59e+06-1.88e+060.0113 7.54e+060.000241-0.252 
25684641.11e+1050191762.81e+06-2.38e+050.011  7.54e+063.78e-06-0.0319
21837441.1e+10 51164622.87e+06-6.87e+050.0109 7.54e+063.1e-05 -0.092 
6595201.49e+1052115544.33e+06-3.67e+060.0108 7.53e+060.00088 -0.492 
13190401.49e+1053043344.43e+06-3.11e+060.0107 7.53e+060.000623-0.417 
11028641.49e+1053946794.53e+06-3.42e+060.0105 7.53e+060.000745-0.458 
6118881.49e+1054865644.62e+06-4.01e+060.0104 7.53e+060.00101 -0.537 
6045601.49e+1055800134.72e+06-4.12e+060.0103 7.53e+060.00105 -0.551 
6851689.86e+0956750543.08e+06-2.39e+060.0103 7.54e+060.000356-0.321 
11902621.12e+1057717143.64e+06-2.45e+060.0102 7.54e+060.000369-0.328 
12201121.48e+1058700204.98e+06-3.76e+060.0101 7.53e+060.00086 -0.503 
36553481.6e+10 59825545.51e+06-1.86e+060.009997.54e+060.000208-0.249 
52203671.51e+1061099025.34e+06-1.22e+050.009977.54e+068.98e-07-0.0164
63708861.58e+1062712285.76e+066.11e+050.009917.54e+062.23e-050.0818
64988991.59e+1064496095.97e+065.25e+050.009917.54e+061.65e-050.0703
70041641.69e+1066394186.52e+064.86e+050.009917.54e+061.41e-050.065 
74297261.7e+10 68366396.76e+066.66e+050.01   7.54e+062.68e-050.0892
75625291.8e+10 70389107.33e+062.33e+050.0101 7.54e+063.31e-060.0313
35133311.79e+1072456807.51e+06-4e+06       0.0103 7.53e+060.000995-0.536 
38320991.79e+1074580157.73e+06-3.9e+06 0.0106 7.53e+060.000973-0.523 
31832441.92e+1076783818.39e+06-5.2e+06 0.0108 7.52e+060.00177 -0.697 
42014502.06e+1079101889.11e+06-4.91e+060.0111 7.53e+060.00161 -0.658 
46849192.2e+10 81563429.86e+06-5.18e+060.0114 7.52e+060.00185 -0.693 
47765492.32e+1084178211.06e+07-5.79e+060.0118 7.52e+060.0024  -0.776 
58279682.43e+1086923371.12e+07-5.41e+060.0124 7.52e+060.00219 -0.724 
73994092.59e+1089737911.21e+07-4.66e+060.013  7.53e+060.00171 -0.625 
83802292.77e+1092538271.3e+07 -4.62e+060.0136 7.53e+060.00176 -0.619 
98343932.91e+1095265581.38e+07-3.93e+060.0143 7.53e+060.00134 -0.527 
101276022.96e+1097900831.42e+07-4.08e+060.0152 7.53e+060.00154 -0.548 
105302263.22e+10100463211.54e+07-4.84e+060.0157 7.53e+060.00225 -0.65  
91377583.41e+10102979831.63e+07-7.14e+060.0164 7.51e+060.00511 -0.958 
145858893.62e+10105493951.73e+07-2.67e+060.0171 7.54e+060.000745-0.359 
166625943.88e+10108041311.84e+07-1.78e+060.0177 7.54e+060.000342-0.239 
179006254.17e+10110624341.97e+07-1.82e+060.0182 7.54e+060.00037 -0.245 
175051274.78e+10113242772.21e+07-4.57e+060.018  7.53e+060.00231 -0.614 
190446795.14e+10115926382.36e+07-4.54e+060.0185 7.53e+060.00235 -0.61  
194451585.18e+10118711022.4e+07 -4.58e+060.02   7.53e+060.00258 -0.616 
238976435.82e+10121621892.65e+07-2.65e+060.0199 7.54e+060.000858-0.356 
226208406.3e+10 124686882.85e+07-5.89e+060.0203 7.52e+060.00434 -0.792 
232604446.85e+10127903133.07e+07-7.49e+060.0207 7.51e+060.00715 -1.01  
273059377.41e+10131228333.3e+07 -5.71e+060.0211 7.52e+060.00425 -0.768 
280325978.02e+10134600353.55e+07-7.45e+060.0215 7.51e+060.00737 -1     
308776618.59e+10137980943.78e+07-6.92e+060.022  7.51e+060.0065  -0.932 
306419409.1e+10 141340603.99e+07-9.27e+060.0226 7.49e+060.012   -1.25  
380483079.66e+10144712154.22e+07-4.13e+060.0232 7.53e+060.00245 -0.557 
347944331.04e+11148194304.51e+07-1.03e+070.0233 7.48e+060.0154  -1.39  
363297711.03e+11151923004.51e+07-8.82e+060.0263 7.5e+06 0.0127  -1.19  
401018851.04e+11155989244.59e+07-5.82e+060.029  7.52e+060.00615 -0.786 
408776841.1e+11 160437364.85e+07-7.57e+060.0304 7.51e+060.011   -1.02  
428692041.2e+11 165220045.24e+07-9.53e+060.0308 7.49e+060.0176  -1.29  
500838701.31e+11170224705.67e+07-6.62e+060.0311 7.51e+060.00858 -0.896 
542693641.43e+11175289616.14e+07-7.1e+06 0.0311 7.51e+060.00987 -0.961 
657851341.58e+11180298246.71e+07-1.28e+060.0301 7.54e+060.000309-0.173 
729610161.73e+11185199417.27e+072.2e+05 0.0292 7.54e+068.84e-060.0297
867186371.92e+11190026607.98e+076.93e+060.0273 7.51e+060.00818 0.936 
906309882.11e+11194849018.68e+073.79e+060.0256 7.53e+060.0023  0.512 
1139336642.33e+11199775089.49e+071.9e+07 0.0238 7.34e+060.0533  2.56  
1134419422.59e+11204876041.04e+089.03e+060.0219 7.49e+060.011   1.22  
1220845862.8e+11 210176191.12e+089.88e+060.0216 7.49e+060.013   1.33  
1142375002.61e+11215627901.06e+088e+06       0.0278 7.5e+06 0.0112  1.08  
1082914402.79e+11221146471.13e+08-4.72e+060.0277 7.53e+060.00385 -0.637 
1262030493.05e+11226612931.23e+083.67e+060.0264 7.53e+060.00222 0.495 
1339881653.09e+11231942521.24e+089.51e+060.0288 7.49e+060.0163  1.29  
1351182083.28e+11237091151.32e+083.56e+060.0288 7.53e+060.00229 0.481 
1565213323.5e+11 242083911.4e+08 1.69e+070.0288 7.38e+060.0513  2.28  
1711020043.76e+11246988211.49e+082.2e+07 0.029  7.27e+060.0879  2.97  
1813328003.98e+11251906471.57e+082.41e+070.0301 7.21e+060.11    3.26  
1787853084.23e+11256906151.66e+081.24e+070.0319 7.45e+060.0312  1.68  
1824553634.53e+11262019541.77e+085.25e+060.035  7.52e+060.00612 0.711 
2019851814.79e+11267203671.87e+081.53e+070.0388 7.41e+060.058   2.08  
2007129444.75e+11272360031.86e+081.48e+070.0387 7.42e+060.0545  2.02  
2159351545.14e+11277350382e+08       1.61e+070.045  7.39e+060.0757  2.2   
2169935275.45e+11282080282.11e+086.01e+060.0518 7.52e+060.0123  0.821 
2152130345.75e+11286509622.22e+08-6.57e+060.06   7.51e+060.0173  -0.902 
2424874386.02e+11290681892.32e+081.1e+07 0.0685 7.47e+060.0562  1.51  
2448325676.38e+11294689232.44e+085.02e+050.0826 7.54e+060.0001460.0699
2318439406.7e+11 298666062.56e+08-2.39e+070.0969 7.19e+060.402   -3.35  
2505606426.99e+11302709652.66e+08-1.56e+070.111  7.39e+060.203   -2.21  
2505613207.39e+11306846522.8e+08 -2.98e+070.135  6.97e+060.946   -4.27  
2722293537.74e+11311046552.93e+08-2.06e+070.157  7.26e+060.558   -2.99  

2.3.2.1 Visualization

GDP<-data1$GDP
Population<-data1$Population
CO2_emissions<-data1$Annual_CO2_emissions
s3d<-scatterplot3d(GDP,Population,CO2_emissions,
                   pch=16,
                   main="3D Scatter Plot visualizing model2")
s3d$plane3d(model2,lty.box='solid',col='blue',draw_polygon = T)

2.3.3 Criteria and Outcomes

  1. The variables are statistically significant as the p-value of each variable is less than the significant level (<0.05).

  2. The high value of the F statistic (1: 402.1, 2: 5921) with a low p-value (<2.2e-16), indicates the potential relationship between the independent variables and the outcome.

  3. The R-squared value (1:0.7567, 2:0.9894 ), is a corrected goodness-of-fit (model accuracy) measure for linear models. It identifies the percentage of variance in the target field that is explained by the input or inputs, both showing an approximately high value(76% and 99%).

Through the modelling analysis, both models produce statistically acceptable results, indicating that increasing GDP and population would bring more CO2 emissions. The increase in CO2 emissions directly affects temperature increase and results in climate change while indirectly leading to sea-level rise.

3 Limitation

The high-positive correlation between the variables is known from the data modelling analysis, but the data set is so limited that the p-value and confidence interval (CO2 emissions and temperature, GDP and CO2 emissions) are very low. Although we cannot directly explain the relationship between CO2 and the four variables, the reality is that as the CO2 emissions in Malaysia increased, the temperature also increased. In addition, the increase in population and economic activity does generate more CO2 emissions in the absence of external governance of carbon dioxide emissions. Also, the excess carbon dioxide emissions could lead to temperature increase. Strictly speaking, this study can only confirm the normal association of the data and cannot be used as an official document for commercial, academic, and governmental purposes, since there are still a lot of variables would affects CO2 emissions and climate change. i.e. energy factors, greenhouse gas emissions. Facts under data analysis have proved that the decentralized collection, quality and quantity of data directly affects the Accuracy and Rigidity of the results.

4 Conclusion

Excess carbon dioxide emissions in the atmosphere have irreversible consequences for Malaysia’s environment, the most notable of which is temperature rise. GDP growth and CO2 emissions are strongly correlated in the last 3 decades, and the positive relationship indicates that CO2 emissions rise as a result of economic growth, highlighting the significance and urgency of environmental governance. To protect our natural environment and the earth, individuals and groups should toward a low-carbon economy and CO2 emissions reduction.

5 References

  • Abbasi, Muhammad Ali, Parveen, Shabana, Khan, Saleem, & Kamal, Muhammad Abdul. (2020). Urbanization and energy consumption effects on carbon dioxide emissions: Evidence from Asian-8 countries using panel data analysis. Environmental Science and Pollution Research International, 27(15), 18029-18043.

  • Deviren, Seyma Akkaya, & Deviren, Bayram. (2016). The relationship between carbon dioxide emission and economic growth: Hierarchical structure methods. Physica A, 451, 429-439.

  • Dechezlepretre, A. (2020). THE ECONOMIC COST OF AIR POLLUTION: EVIDENCE FROM EUROPE. From https://www.oecd.org/officialdocuments/publicdisplaydocumentpdf/?cote=ECO/WKP(2019)54&docLanguage=En

  • Lindsey, R. (2021). Retrieved 7 April 2022, from https://drupal-www.climate.woc.noaa.gov/news-features/understanding-climate/climate-change-atmospheric-carbon-dioxide

  • Zhang, Lingyun, Li, Zecheng, Kirikkaleli, Dervis, Adebayo, Tomiwa Sunday, Adeshola, Ibrahim, & Akinsola, Gbenga Daniel. (2021). Modeling CO2 emissions in Malaysia: An application of Maki cointegration and wavelet coherence tests. Environmental Science and Pollution Research International, 28(20), 26030-26044.