Lembaga : UIN Maulana Malik Ibrahim Malang
Jurusan : Teknik Informatika

1. Pengertian Regresi Linear

Regresi linier adalah model paling sederhana yang paling sering dijelaskan dalam statistik. Modelnya sangat sederhana dimana kita dapat mencoba membangun model dengan pendekatan linier menggunakan prinsip meminimalkan jumlah sisa kuadrat dalam data. Secara artian regresi merupakan suatu metode yang menentukan sebab-akibat antara variabel satu dengan lainnya.

2. Data Riwayat Covid-19 dan Google Mobility Index di Jakarta pada Bulan Juli 2020

library(readxl)
Data<- read_excel(path = "DataMobility1-7.xlsx")
Data
library(ggplot2)
library(reshape2)
x <- Data$Meninggal
retail <- Data$retail_and_recreation_percent_change_from_baseline
grocery <- Data$grocery_and_pharmacy_percent_change_from_baseline
park <- Data$parks_percent_change_from_baseline
station <- Data$transit_stations_percent_change_from_baseline
workplace <- Data$workplaces_percent_change_from_baseline
residental <- Data$residential_percent_change_from_baseline
df <- data.frame(x, retail, grocery, park, station, workplace,residental  )

# melt the data to a long format
df2 <- melt(data = df, id.vars = "x")

# plot, using the aesthetics argument 'colour'
ggplot(data = df2, aes(x = x, y = value, colour = variable))+
  geom_point() +
  geom_line() + 
  theme(legend.justification = "top") +
  labs(title = "GOOGLE MOBILITY INDEX", 
         subtitle = "Provinsi DKI Jakarta Indonesia Bulan April 2020", 
         y = "Mobility Index", x = "Data Meninggal") +
theme(axis.text.x = element_text(angle = -90))

3. Regresi Linear Berganda

model <- lm(Data$Meninggal~Data$retail_and_recreation_percent_change_from_baseline+Data$grocery_and_pharmacy_percent_change_from_baseline+Data$parks_percent_change_from_baseline+Data$transit_stations_percent_change_from_baseline+Data$workplaces_percent_change_from_baseline+Data$residential_percent_change_from_baseline)
model
## 
## Call:
## lm(formula = Data$Meninggal ~ Data$retail_and_recreation_percent_change_from_baseline + 
##     Data$grocery_and_pharmacy_percent_change_from_baseline + 
##     Data$parks_percent_change_from_baseline + Data$transit_stations_percent_change_from_baseline + 
##     Data$workplaces_percent_change_from_baseline + Data$residential_percent_change_from_baseline)
## 
## Coefficients:
##                                             (Intercept)  
##                                               1708.4186  
## Data$retail_and_recreation_percent_change_from_baseline  
##                                                 11.3836  
##  Data$grocery_and_pharmacy_percent_change_from_baseline  
##                                                -16.8611  
##                 Data$parks_percent_change_from_baseline  
##                                                  4.4519  
##      Data$transit_stations_percent_change_from_baseline  
##                                                 12.1658  
##            Data$workplaces_percent_change_from_baseline  
##                                                 -0.5654  
##           Data$residential_percent_change_from_baseline  
##                                                 -0.1557

4. Uji Asumsi Klasik

4.1 Uji Normalitas

Uji normalitas dalam contoh regresi dipakai untuk menguji apakah nilai residual yang didapatkan menurut regresi terdistribusi secara normal atau tidak. Model regresi yang baik mempunyai nilai residual yang terdistribusi secara normal. Pada contoh berikut akan dilakukakan uji normalitas menggunakan metode One-sample Kolmogorov-Smirnov

ks.test(model$residuals, ecdf(model$residuals))
## 
##  One-sample Kolmogorov-Smirnov test
## 
## data:  model$residuals
## D = 0.032258, p-value = 1
## alternative hypothesis: two-sided
  • Berdasarkan dasar teori atau syarat dari uji normalitas adalah jika nilai p-value > 0.05 maka data berdistribusi normal. Namun apabila sebaliknya maka data tidak berdstribusi normal. Sehingga dapat disimpulkan bahwa model regresi diatas berdistribusi normal dikarenakan nilai p-value sama dengan 1 dimana > 0.05.

4.2 Uji Multikolinearitas

library(olsrr)
ols_vif_tol(model)
  • Berdasarkan dasar teori uji multikolinearitas jika nilai tolerance > 0.1 dan nilai VIF < 10 maka tidak terjadi gejala multikolinearitas. Sehingga dapat disimpulkan pada regresi diatas terjadi gejala multikolinearitas dikarenakan dari semua variabel independent memiliki nilai tolerance < 0.1 dan nilai VIF > 10.

4.3 Uji Autokorelasi

library(lmtest)
library(car)
dwtest(model)
## 
##  Durbin-Watson test
## 
## data:  model
## DW = 1.1796, p-value = 0.002389
## alternative hypothesis: true autocorrelation is greater than 0
  • Berdasarkan dasar teori atau syarat dari uji autokorelasi adalah jika nilai p-value > 0.05 maka tidak terjadi autokorelasi. Namun apabila sebaliknya maka terjadi autokorelasi. Sehingga dapat disimpulkan pada model diatas terdapat gejala autokorelasi.

4.4 Uji Homogenitas

Uji homogenitas dimaksudkan untuk memperlihatkan bahwa dua atau lebih kelompok data sampel berasal dari populasi yang memiliki variansi yang sama. Pada contoh berikut akan dilakukakan uji homogenitas menggunakan metode studentized Breusch-Pagan.

bptest(model)
## 
##  studentized Breusch-Pagan test
## 
## data:  model
## BP = 5.8392, df = 6, p-value = 0.4414
  • Berdasarkan dasar teori atau syarat dari uji homogenitas jika nilai p-value > 0.05 maka variansi setiap sampel sama (homogen). Namun apabila sebaliknya nilai p-value< 0.05, maka variansi setiap sampel tidak sama (tidak homogen). Sehingga dapat disimpulkan bahwa data tersebut homogen.

5. Interpretasi Hasil Regresi Berganda

Dalam menampilkan hasil regresi kita dapat menggunakan fungsi summary.

summary(model)
## 
## Call:
## lm(formula = Data$Meninggal ~ Data$retail_and_recreation_percent_change_from_baseline + 
##     Data$grocery_and_pharmacy_percent_change_from_baseline + 
##     Data$parks_percent_change_from_baseline + Data$transit_stations_percent_change_from_baseline + 
##     Data$workplaces_percent_change_from_baseline + Data$residential_percent_change_from_baseline)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -40.248 -18.209   1.911  18.859  39.181 
## 
## Coefficients:
##                                                          Estimate Std. Error
## (Intercept)                                             1708.4186   177.0458
## Data$retail_and_recreation_percent_change_from_baseline   11.3836     3.1336
## Data$grocery_and_pharmacy_percent_change_from_baseline   -16.8611     2.7846
## Data$parks_percent_change_from_baseline                    4.4519     0.9526
## Data$transit_stations_percent_change_from_baseline        12.1658     2.2539
## Data$workplaces_percent_change_from_baseline              -0.5654     0.8870
## Data$residential_percent_change_from_baseline             -0.1557     0.5318
##                                                         t value       Pr(>|t|)
## (Intercept)                                               9.650 0.000000000982
## Data$retail_and_recreation_percent_change_from_baseline   3.633        0.00133
## Data$grocery_and_pharmacy_percent_change_from_baseline   -6.055 0.000002976404
## Data$parks_percent_change_from_baseline                   4.673 0.000095298066
## Data$transit_stations_percent_change_from_baseline        5.398 0.000015239078
## Data$workplaces_percent_change_from_baseline             -0.638        0.52983
## Data$residential_percent_change_from_baseline            -0.293        0.77229
##                                                            
## (Intercept)                                             ***
## Data$retail_and_recreation_percent_change_from_baseline ** 
## Data$grocery_and_pharmacy_percent_change_from_baseline  ***
## Data$parks_percent_change_from_baseline                 ***
## Data$transit_stations_percent_change_from_baseline      ***
## Data$workplaces_percent_change_from_baseline               
## Data$residential_percent_change_from_baseline              
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 24.59 on 24 degrees of freedom
##   (1048520 observations deleted due to missingness)
## Multiple R-squared:  0.8517, Adjusted R-squared:  0.8146 
## F-statistic: 22.97 on 6 and 24 DF,  p-value: 0.000000007693

  • Nilai Minimum adalah -40.248

  • Nilai Quartal ke-1 adalah -18.209

  • Nilai Tengah adalah 1.911

  • Nilai Quartal ke-3 adalah 263.8

  • Nilai Maksimum adalah 39.181

  • Dasar teori yang digunakan dalam interpretasi hasilregre adalah apabila nilai signifikansi (Pr(>|t|)) < 0.05 maka variabel independent (variabel x) secara parsial berpengaruh terhadap variable dependent (variable y). Sehingga dapat disimpulkan :

    • Variable retail_and_recreation_percent_change_from_baseline signifikan berpengaruh terhadap variable positif.
    • Variable grocery_and_pharmacy_percent_change_from_baseline tidak signifikan berpengaruh terhadap variable positif.
    • Variable parks_percent_change_from_baseline signifikan berpengaruh terhadap variable positif.
    • Variable transit_stations_percent_change_from_baselineidak signifikan berpengaruh terhadap variable positif.
    • Variable workplaces_percent_change_from_baseline signifikan berpengaruh terhadap variable positif.
    • Variable residential_percent_change_from_baseline signifikan berpengaruh terhadap variable positif.

*Selain itu kita dapat disimpulkan apakah seluruh variable independent berpengaruh secara simultan terhadap variable dependent dimana dasar dari penentuan tersebut yaitu apabila nilai p-value dari F-statistic < 0.05 maka keseluruahan variable independent signifikan berpengaruh secara simultan (bersama-sama) terhadap variable dependent (variable positif). Besar pengaruh tersebut dapat dilihat dari nilai R-squared dimana pada model diatas bernilai 0.7864. Sehingga dapat disimpulkan variable independent signifikan berpengaruh terhadap variable dependent sebesar 78.64 %.

6. Plot Hasil Regresi Berganda

6.1 Plot Pengaruh Variable Retail and Recreation Percent Change From Baseline terhadap Variable Positif

plot(Data$retail_and_recreation_percent_change_from_baseline, Data$Meninggal, col = "dodgerblue")

6.2 Plot Pengaruh Variable Grocery and Pharmacy Percent Change From Baseline terhadap Variable Positif

plot(Data$grocery_and_pharmacy_percent_change_from_baseline, Data$Meninggal, col = "red")

6.3 Plot Pengaruh Variable Parks Percent Change From Baseline terhadap Variable Positif

plot(Data$parks_percent_change_from_baseline, Data$Meninggal, col = "darkorange")

6.4 Plot Pengaruh Variable Stations Percent Change From Baseline terhadap Variable Positif

plot(Data$transit_stations_percent_change_from_baseline, Data$Meninggal, col = "darkgreen")

6.5 Plot Pengaruh Variable Workplaces Percent Change From Baseline terhadap Variable Positif

plot(Data$workplaces_percent_change_from_baseline, Data$Meninggal, col = "blueviolet")

6.6 Plot Pengaruh Variable Residential Percent Change From Baseline terhadap Variable Positif

plot(Data$residential_percent_change_from_baseline, Data$Meninggal, col = "darkcyan")

6.7 Plot Hasil Regresi Linear Berganda

plot(model)