Penelitian ini bertujuan untuk menganalisis pengaruh harga emas dan kurs jual terhadap inflasi bulanan Indonesia tahun 2025 menggunakan regresi linear berganda.
Variabel yang digunakan: * Y : Inflasi (%) * X1 : Harga Emas * X2 : Kurs Jual
library(readxl)
data_regresi <- read_excel("D:\\SEMESTER 6\\KOMLAN\\data regresi.xlsx")
data_regresi
## # A tibble: 12 × 4
## Bulan Inflasi `Harga Emas` `Kurs Jual`
## <chr> <dbl> <dbl> <dbl>
## 1 Januari 0.76 2845. 16343.
## 2 Februari -0.09 2896. 16420.
## 3 Maret 1.03 3150. 16538.
## 4 April 1.95 3319. 16904.
## 5 Mei 1.6 3344. 16523.
## 6 Juni 1.87 3308. 16392.
## 7 Juli 2.37 3349. 16358.
## 8 Agustus 2.31 3512 16390.
## 9 September 2.65 3873. 16595.
## 10 Oktober 2.86 4013. 16690.
## 11 November 2.72 4238. 16787.
## 12 Desember 2.92 4341. 16783.
summary(data_regresi)
## Bulan Inflasi Harga Emas Kurs Jual
## Length:12 Min. :-0.090 Min. :2845 Min. :16343
## Class :character 1st Qu.: 1.458 1st Qu.:3268 1st Qu.:16392
## Mode :character Median : 2.130 Median :3346 Median :16531
## Mean : 1.913 Mean :3516 Mean :16560
## 3rd Qu.: 2.667 3rd Qu.:3908 3rd Qu.:16713
## Max. : 2.920 Max. :4341 Max. :16904
Berdasarkan hasil statistik deskriptif, data yang digunakan terdiri dari 12 observasi bulanan selama tahun 2025. Variabel inflasi memiliki nilai minimum sebesar -0,09 dan maksimum sebesar 2,92. Nilai rata-rata inflasi selama tahun 2025 adalah 1,913 dengan median sebesar 2,130. Hal ini menunjukkan bahwa secara umum inflasi berada pada kisaran positif dan relatif stabil, meskipun sempat terjadi deflasi pada salah satu bulan.
Variabel harga emas memiliki nilai minimum sebesar 2.845 dan maksimum sebesar 4.341 dengan rata-rata sebesar 3.516. Median harga emas sebesar 3.346 menunjukkan bahwa sebagian besar harga emas berada di sekitar kisaran tersebut, dengan kecenderungan meningkat di akhir tahun.
Variabel kurs jual memiliki nilai minimum sebesar 16.343 dan maksimum sebesar 16.904 dengan rata-rata sebesar 16.560. Median sebesar 16.531 menunjukkan bahwa nilai tukar relatif stabil meskipun terdapat fluktuasi sepanjang tahun.
Bentuk umum persamaan regresi linier: \[ y = \beta_0 + \beta_1 X_1 + \beta_2 X_2 + \epsilon \]
model <- lm(Inflasi ~ `Harga Emas` + `Kurs Jual`, data = data_regresi)
summary(model)
##
## Call:
## lm(formula = Inflasi ~ `Harga Emas` + `Kurs Jual`, data = data_regresi)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.96325 -0.26388 0.04325 0.25100 0.63919
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 5.4480046 15.9789345 0.341 0.7410
## `Harga Emas` 0.0018127 0.0003944 4.596 0.0013 **
## `Kurs Jual` -0.0005983 0.0010169 -0.588 0.5707
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4932 on 9 degrees of freedom
## Multiple R-squared: 0.7742, Adjusted R-squared: 0.724
## F-statistic: 15.43 on 2 and 9 DF, p-value: 0.001235
Berdasarkan Tabel Coefficients didapatkan model regresi :
\[ Inflasi = 5.448 + 0.0018 Harga\ Emas -6\times 10^{-4} Kurs\ Jual \]
shapiro.test(residuals(model))
##
## Shapiro-Wilk normality test
##
## data: residuals(model)
## W = 0.95841, p-value = 0.7609
Jika P-Value > 0.05, maka residual berdistribusi normal. Berdasarkan test normalitas menggunakan Shapiro-Wilk didapatkan nilai P-Value = 0.7609 > 0.05, sehingga residual berdistribusi normal.
qqnorm(residuals(model))
qqline(residuals(model), col="darkblue")
Berdasarkan Q - Q Plot, plot – plot secara umum mengikuti garis lurus, sehingga dapat disimpulkan bahwa residual data berdistribusi normal.
plot(fitted(model), residuals(model))
abline(h = 0, col = "darkblue")
Berdasarkan gafik Scatterplot, dapat dilihat bahwa plot – plot menyebar secara acak disekitar nilai nol dan tidak membentuk pola tertentu, sehingga asumsi linieritas terpenuhi.
res <- residuals(model)
abs_res <- abs(res)
glejser <- lm(abs_res ~ `Harga Emas` + `Kurs Jual`, data = data_regresi)
summary(glejser)
##
## Call:
## lm(formula = abs_res ~ `Harga Emas` + `Kurs Jual`, data = data_regresi)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.33582 -0.16096 -0.02039 0.10922 0.54217
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -5.7972266 9.4033221 -0.617 0.553
## `Harga Emas` -0.0002347 0.0002321 -1.011 0.338
## `Kurs Jual` 0.0004201 0.0005984 0.702 0.500
##
## Residual standard error: 0.2902 on 9 degrees of freedom
## Multiple R-squared: 0.1024, Adjusted R-squared: -0.0971
## F-statistic: 0.5132 on 2 and 9 DF, p-value: 0.6151
Jika p-value > 0.05, tidak ada hubungan signifikan antara nilai prediktor dan absolut residual atau homoskedastisitas terpenuhi. Pada Tabel Coefficients, didapatkan nilai P-Value (0.6151) > 0.05, sehingga homoskedastisitas terpenuhi.
library(lmtest)
## Warning: package 'lmtest' was built under R version 4.3.3
## Loading required package: zoo
## Warning: package 'zoo' was built under R version 4.3.3
##
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
##
## as.Date, as.Date.numeric
dwtest(model)
##
## Durbin-Watson test
##
## data: model
## DW = 1.3879, p-value = 0.04274
## alternative hypothesis: true autocorrelation is greater than 0
Hasil uji Durbin-Watson menunjukkan nilai DW sebesar 1.3879 dengan p-value 0.04274, yang artinya terdapat autokorelasi positif pada residual.
library(car)
## Warning: package 'car' was built under R version 4.3.3
## Loading required package: carData
## Warning: package 'carData' was built under R version 4.3.3
vif(model)
## `Harga Emas` `Kurs Jual`
## 1.717215 1.717215
Nilai Variance Inflation Factor diatas pada Harga Emas sebesar 1.717215, Variance Inflation Factor dari Kurs Jual sebesar 1.717215. Berdasarkan nilai tersebut, semua Variance Inflation Factor kurang dari 10 maka dapat disimpulkan bahwa tidak terjadi multikolinearitas.
f_stat <- summary(model)$fstatistic[1]
df1 <- summary(model)$fstatistic[2]
df2 <- summary(model)$fstatistic[3]
p_value_f <- pf(f_stat, df1, df2, lower.tail = FALSE)
f_stat
## value
## 15.42921
p_value_f
## value
## 0.001235239
\[H_1 : 𝛽 ≠ 0 (model\ signifikan)\] * Taraf signifikansi : \[α = 5\% \] * Statistik Uji Berdasarkan hasil estimasi regresi diperoleh: \[Fhitung=15.42921\] dengan \[p-value= 0.001235239\] * Daerah Kritis \[Tolak\ H_0\ jika\ p-value < 0.05\] * Keputusan Karena p-value (0.001235) < 0.05, maka \[H_0\ ditolak\] * Kesimpulan Pada taraf signifikansi 5%, model regresi signifikan secara simultan. Artinya, variabel Harga Emas dan Kurs Jual secara bersama-sama berpengaruh terhadap Inflasi, sehingga model layak digunakan dalam analisis.
coef_summary <- summary(model)$coefficients
coef_summary
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 5.4480046031 1.597893e+01 0.3409492 0.740966076
## `Harga Emas` 0.0018126918 3.944473e-04 4.5955237 0.001298772
## `Kurs Jual` -0.0005983178 1.016858e-03 -0.5883987 0.570733637
Berdasarkan hasil analisis awal, variabel Kurs Jual tidak signifikan terhadap Inflasi, sehingga model belum bisa digunakan. Selain itu, uji Durbin-Watson menunjukkan adanya autokorelasi positif pada residual, yang berarti asumsi non-autokorelasi tidak terpenuhi. Hal ini menunjukkan bahwa model regresi awal belum memenuhi asumsi klasik dan hasil uji t. Oleh karena itu, perlu dilakukan uji ulang model dengan menghapus variabel Kurs Jual untuk mendapatkan model yang lebih sederhana, signifikan, dan memenuhi asumsi non-autokorelasi sehingga analisis regresi menjadi lebih valid.
model2 <- lm(Inflasi ~ `Harga Emas`, data = data_regresi)
summary(model2)
##
## Call:
## lm(formula = Inflasi ~ `Harga Emas`, data = data_regresi)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.97204 -0.29752 0.04656 0.31852 0.73525
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -3.932965 1.032274 -3.810 0.003429 **
## `Harga Emas` 0.001663 0.000291 5.714 0.000195 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4768 on 10 degrees of freedom
## Multiple R-squared: 0.7655, Adjusted R-squared: 0.7421
## F-statistic: 32.65 on 1 and 10 DF, p-value: 0.0001946
residuals_model2 <- residuals(model2)
shapiro.test(residuals_model2)
##
## Shapiro-Wilk normality test
##
## data: residuals_model2
## W = 0.96823, p-value = 0.8914
Jika P-Value > 0.05, maka residual berdistribusi normal. Berdasarkan test normalitas menggunakan Shapiro-Wilk didapatkan nilai P-Value = 0.7609 > 0.05, sehingga residual berdistribusi normal.
plot(fitted(model2), residuals(model2))
abline(h = 0, col = "darkblue")
Berdasarkan gafik Scatterplot, dapat dilihat bahwa plot – plot menyebar secara acak disekitar nilai nol dan tidak membentuk pola tertentu, sehingga asumsi linieritas terpenuhi.
abs_res <- abs(residuals_model2)
glejser <- lm(abs_res ~ `Harga Emas`, data = data_regresi)
summary(glejser)
##
## Call:
## lm(formula = abs_res ~ `Harga Emas`, data = data_regresi)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.38114 -0.16433 -0.03344 0.11522 0.55876
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.7328066 0.6173085 1.187 0.263
## `Harga Emas` -0.0001103 0.0001740 -0.634 0.540
##
## Residual standard error: 0.2851 on 10 degrees of freedom
## Multiple R-squared: 0.03865, Adjusted R-squared: -0.05749
## F-statistic: 0.402 on 1 and 10 DF, p-value: 0.5403
Jika p-value > 0.05, tidak ada hubungan signifikan antara nilai prediktor dan absolut residual atau homoskedastisitas terpenuhi. Pada Tabel Coefficients, didapatkan nilai P-Value (0.5403) > 0.05, sehingga homoskedastisitas terpenuhi.
library(lmtest)
dwtest(model2)
##
## Durbin-Watson test
##
## data: model2
## DW = 1.1695, p-value = 0.02456
## alternative hypothesis: true autocorrelation is greater than 0
Hasil uji Durbin-Watson menunjukkan nilai DW sebesar 1.3879 dengan p-value 0.04274, yang artinya terdapat autokorelasi positif pada residual.
f_stat <- summary(model2)$fstatistic[1]
df1 <- summary(model2)$fstatistic[2]
df2 <- summary(model2)$fstatistic[3]
p_value_f <- pf(f_stat, df1, df2, lower.tail = FALSE)
f_stat
## value
## 32.6466
p_value_f
## value
## 0.0001946052
\[H_1 : 𝛽 ≠ 0 (model\ signifikan)\] * Taraf Signifikansi \[α = 5\% \] * Statistik Uji Berdasarkan hasil estimasi regresi diperoleh: \[Fhitung=32.6466 \] dengan \[p-value= 0.0001946052\] * Daerah Kritis \[Tolak\ H_0\ jika\ p-value < 0.05\] * Keputusan Karena P-Value 0.0001946052 < 0.05, maka \[H_0\ ditolak\]. * Kesimpulan Pada taraf signifikansi 5%, model regresi signifikan secara simultan.Artinya, variabel Harga Emas berpengaruh terhadap Inflasi, sehingga model layak digunakan dalam analisis.
\[ Inflasi = -3.933 + 0.0017 \ Harga\ Emas \]