Berikut data volume produksi air minum PDAM Kabupaten Serang, Januari 2022–Desember 2024 (ribu m³/bulan):
| Tahun | Jan | Feb | Mar | Apr | Mei | Jun | Jul | Agu | Sep | Okt | Nov | Des |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 2022 | 410 | 415 | 405 | 420 | 430 | 445 | 455 | 440 | 425 | 420 | 418 | 435 |
| 2023 | 440 | 445 | 438 | 452 | 462 | 478 | 488 | 472 | 458 | 452 | 450 | 468 |
| 2024 | 470 | 476 | 468 | 483 | 494 | 510 | 521 | 504 | 489 | 483 | 480 | 499 |
# 1. Buat objek ts dari data di atas, lalu lakukan dekomposisi additive menggunakan decompose().
baja <- c(410, 415, 405, 420, 430, 445, 455, 440, 425, 420, 418, 435,
440, 445, 438, 452, 462, 478, 488, 472, 458, 452, 450, 468, 470, 476, 468, 483, 494, 510, 521, 504, 489, 483, 480, 499)
# Membuat objek time series
baja_ts <- ts(baja,start = c(2022, 1),frequency = 12)
print(baja_ts)
## Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
## 2022 410 415 405 420 430 445 455 440 425 420 418 435
## 2023 440 445 438 452 462 478 488 472 458 452 450 468
## 2024 470 476 468 483 494 510 521 504 489 483 480 499
# Dekomposisi Additive
dekom_add <- decompose(baja_ts, type = "additive")
# Menampilkan hasil dekomposisi
dekom_add
## $x
## Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
## 2022 410 415 405 420 430 445 455 440 425 420 418 435
## 2023 440 445 438 452 462 478 488 472 458 452 450 468
## 2024 470 476 468 483 494 510 521 504 489 483 480 499
##
## $seasonal
## Jan Feb Mar Apr May
## 2022 -4.53993056 -1.74826389 -11.91493056 -0.06076389 7.83506944
## 2023 -4.53993056 -1.74826389 -11.91493056 -0.06076389 7.83506944
## 2024 -4.53993056 -1.74826389 -11.91493056 -0.06076389 7.83506944
## Jun Jul Aug Sep Oct
## 2022 21.21006944 27.75173611 9.73090278 -7.35243056 -15.47743056
## 2023 21.21006944 27.75173611 9.73090278 -7.35243056 -15.47743056
## 2024 21.21006944 27.75173611 9.73090278 -7.35243056 -15.47743056
## Nov Dec
## 2022 -20.12326389 -5.31076389
## 2023 -20.12326389 -5.31076389
## 2024 -20.12326389 -5.31076389
##
## $trend
## Jan Feb Mar Apr May Jun Jul Aug
## 2022 NA NA NA NA NA NA 427.7500 430.2500
## 2023 443.7083 446.4167 449.1250 451.8333 454.5000 457.2083 459.8333 462.3750
## 2024 475.4583 478.1667 480.7917 483.3750 485.9167 488.4583 NA NA
## Sep Oct Nov Dec
## 2022 432.8750 435.5833 438.2500 440.9583
## 2023 464.9167 467.4583 470.0833 472.7500
## 2024 NA NA NA NA
##
## $random
## Jan Feb Mar Apr May Jun
## 2022 NA NA NA NA NA NA
## 2023 0.83159722 0.33159722 0.78993056 0.22743056 -0.33506944 -0.41840278
## 2024 -0.91840278 -0.41840278 -0.87673611 -0.31423611 0.24826389 0.33159722
## Jul Aug Sep Oct Nov Dec
## 2022 -0.50173611 0.01909722 -0.52256944 -0.10590278 -0.12673611 -0.64756944
## 2023 0.41493056 -0.10590278 0.43576389 0.01909722 0.03993056 0.56076389
## 2024 NA NA NA NA NA NA
##
## $figure
## [1] -4.53993056 -1.74826389 -11.91493056 -0.06076389 7.83506944
## [6] 21.21006944 27.75173611 9.73090278 -7.35243056 -15.47743056
## [11] -20.12326389 -5.31076389
##
## $type
## [1] "additive"
##
## attr(,"class")
## [1] "decomposed.ts"
# Plot dekomposisi
plot(dekom_add)
Hasil dekomposisi additive menunjukkan bahwa volume produksi air minum PDAM Kabupaten Serang memiliki tren yang meningkat selama periode 2022–2024. Selain itu, terdapat pola musiman yang berulang setiap 12 bulan, dengan produksi cenderung meningkat pada pertengahan tahun, terutama sekitar Juni–Juli, kemudian menurun pada beberapa bulan berikutnya. Komponen residual menunjukkan variasi yang tidak dapat dijelaskan oleh tren dan pola musiman.
# 2. Tampilkan dan gambarkan komponen tren. Apakah produksi air minum menunjukkan tren naik?
# Menampilkan komponen trend
dekom_add$trend
## Jan Feb Mar Apr May Jun Jul Aug
## 2022 NA NA NA NA NA NA 427.7500 430.2500
## 2023 443.7083 446.4167 449.1250 451.8333 454.5000 457.2083 459.8333 462.3750
## 2024 475.4583 478.1667 480.7917 483.3750 485.9167 488.4583 NA NA
## Sep Oct Nov Dec
## 2022 432.8750 435.5833 438.2500 440.9583
## 2023 464.9167 467.4583 470.0833 472.7500
## 2024 NA NA NA NA
# Grafik trend
plot(dekom_add$trend,
main = "Trend Volume Produksi PDAM",
ylab = "Trend",
xlab = "Tahun"
)
Komponen trend menunjukkan kecenderungan peningkatan volume produksi dari tahun 2022 hingga 2024. Artinya, secara jangka panjang volume produksi PDAM Kabupaten Serang mengalami pertumbuhan.
# 3. Tampilkan 12 indeks musiman. Bulan apa yang cenderung tertinggi dan terendah? Berikan dugaan penjelasannya.
dekom_add$seasonal[1:12]
## [1] -4.53993056 -1.74826389 -11.91493056 -0.06076389 7.83506944
## [6] 21.21006944 27.75173611 9.73090278 -7.35243056 -15.47743056
## [11] -20.12326389 -5.31076389
#Hasil lebih rapih
bulan <- c(
"Jan", "Feb", "Mar", "Apr", "Mei", "Jun",
"Jul", "Agu", "Sep", "Okt", "Nov", "Des"
)
indeks_add <- data.frame(
Bulan = bulan,
Indeks_Musiman = round(dekom_add$seasonal[1:12], 2)
)
indeks_add
## Bulan Indeks_Musiman
## 1 Jan -4.54
## 2 Feb -1.75
## 3 Mar -11.91
## 4 Apr -0.06
## 5 Mei 7.84
## 6 Jun 21.21
## 7 Jul 27.75
## 8 Agu 9.73
## 9 Sep -7.35
## 10 Okt -15.48
## 11 Nov -20.12
## 12 Des -5.31
Indeks musiman menunjukkan bahwa Juli memiliki efek musiman tertinggi, yaitu sekitar +27,75 ribu m³, sehingga produksi pada bulan Juli cenderung lebih tinggi sekitar 27,75 ribu m³ dibandingkan nilai trendnya.
Sebaliknya, November memiliki efek musiman terendah, yaitu sekitar −20,12 ribu m³, sehingga produksi pada bulan November cenderung lebih rendah sekitar 20,12 ribu m³ dibandingkan nilai trendnya.
Dengan demikian, Juni–Juli merupakan periode dengan kecenderungan produksi tinggi, sedangkan Oktober–November cenderung rendah.
# 4.Periksa komponen residu (ringkasan statistik dan simpangan bakunya). Apakah residu tampak berfluktuasi acak di sekitar nol?
# Menampilkan residual
dekom_add$random
## Jan Feb Mar Apr May Jun
## 2022 NA NA NA NA NA NA
## 2023 0.83159722 0.33159722 0.78993056 0.22743056 -0.33506944 -0.41840278
## 2024 -0.91840278 -0.41840278 -0.87673611 -0.31423611 0.24826389 0.33159722
## Jul Aug Sep Oct Nov Dec
## 2022 -0.50173611 0.01909722 -0.52256944 -0.10590278 -0.12673611 -0.64756944
## 2023 0.41493056 -0.10590278 0.43576389 0.01909722 0.03993056 0.56076389
## 2024 NA NA NA NA NA NA
# Ringkasan residual
summary(dekom_add$random, na.rm = TRUE)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## -0.9184 -0.4184 -0.0434 -0.0434 0.3316 0.8316 12
# Simpangan baku residual
sd(dekom_add$random, na.rm = TRUE)
## [1] 0.4862621
# Grafik residual
plot(
dekom_add$random,
main = "Residual Dekomposisi Additive",
ylab = "Residual",
xlab = "Tahun"
)
# Garis nol
abline(h = 0)
Residual merupakan bagian variasi data yang tidak dapat dijelaskan oleh trend dan pola musiman. Dari grafik, residual berfluktuasi di sekitar nilai nol dan relatif kecil. Hal ini menunjukkan bahwa sebagian besar pola sistematis pada data telah dapat dijelaskan oleh komponen trend dan seasonal.
Kesimpulan: Berdasarkan dekomposisi additive, volume produksi air minum PDAM Kabupaten Serang periode Januari 2022–Desember 2024 menunjukkan tren meningkat dan memiliki pola musiman tahunan. Produksi cenderung tinggi pada Juni–Juli dan rendah pada Oktober–November. Residual relatif kecil dan berfluktuasi di sekitar nol. Dengan demikian, pola utama volume produksi dapat dijelaskan oleh komponen trend dan musiman.
# Dekomposisi Multipicative
dekom_mult <- decompose(baja_ts, type = "multiplicative")
print(dekom_mult)
## $x
## Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
## 2022 410 415 405 420 430 445 455 440 425 420 418 435
## 2023 440 445 438 452 462 478 488 472 458 452 450 468
## 2024 470 476 468 483 494 510 521 504 489 483 480 499
##
## $seasonal
## Jan Feb Mar Apr May Jun Jul
## 2022 0.9901875 0.9962547 0.9744168 0.9999039 1.0166776 1.0449005 1.0625939
## 2023 0.9901875 0.9962547 0.9744168 0.9999039 1.0166776 1.0449005 1.0625939
## 2024 0.9901875 0.9962547 0.9744168 0.9999039 1.0166776 1.0449005 1.0625939
## Aug Sep Oct Nov Dec
## 2022 1.0218486 0.9835709 0.9656814 0.9556379 0.9883262
## 2023 1.0218486 0.9835709 0.9656814 0.9556379 0.9883262
## 2024 1.0218486 0.9835709 0.9656814 0.9556379 0.9883262
##
## $trend
## Jan Feb Mar Apr May Jun Jul Aug
## 2022 NA NA NA NA NA NA 427.7500 430.2500
## 2023 443.7083 446.4167 449.1250 451.8333 454.5000 457.2083 459.8333 462.3750
## 2024 475.4583 478.1667 480.7917 483.3750 485.9167 488.4583 NA NA
## Sep Oct Nov Dec
## 2022 432.8750 435.5833 438.2500 440.9583
## 2023 464.9167 467.4583 470.0833 472.7500
## 2024 NA NA NA NA
##
## $random
## Jan Feb Mar Apr May Jun Jul
## 2022 NA NA NA NA NA NA 1.0010461
## 2023 1.0014693 1.0005740 1.0008342 1.0004650 0.9998269 1.0005500 0.9987391
## 2024 0.9983158 0.9992111 0.9989510 0.9993202 0.9999583 0.9992351 NA
## Aug Sep Oct Nov Dec
## 2022 1.0007953 0.9982074 0.9984910 0.9980699 0.9981398
## 2023 0.9989899 1.0015778 1.0012941 1.0017152 1.0016454
## 2024 NA NA NA NA NA
##
## $figure
## [1] 0.9901875 0.9962547 0.9744168 0.9999039 1.0166776 1.0449005 1.0625939
## [8] 1.0218486 0.9835709 0.9656814 0.9556379 0.9883262
##
## $type
## [1] "multiplicative"
##
## attr(,"class")
## [1] "decomposed.ts"
plot(dekom_mult)
Hasil dekomposisi multiplicative juga menunjukkan adanya tren meningkat dan pola musiman tahunan. Namun, komponen seasonal dinyatakan dalam bentuk indeks atau proporsi, bukan dalam satuan ribu m³.
# Trend Multipicative
dekom_mult$trend
## Jan Feb Mar Apr May Jun Jul Aug
## 2022 NA NA NA NA NA NA 427.7500 430.2500
## 2023 443.7083 446.4167 449.1250 451.8333 454.5000 457.2083 459.8333 462.3750
## 2024 475.4583 478.1667 480.7917 483.3750 485.9167 488.4583 NA NA
## Sep Oct Nov Dec
## 2022 432.8750 435.5833 438.2500 440.9583
## 2023 464.9167 467.4583 470.0833 472.7500
## 2024 NA NA NA NA
plot(
dekom_mult$trend,
main = "Trend Volume Produksi PDAM - Multiplicative",
ylab = "Trend",
xlab = "Tahun"
)
Trend menunjukkan adanya peningkatan volume prediksi PDAM dari tahun 2022 hingga 2024.
# 12 indeks musiman multipicative
indeks_mult <- data.frame(
Bulan = bulan,
Indeks_Musiman = round(dekom_mult$seasonal[1:12], 4)
)
indeks_mult
## Bulan Indeks_Musiman
## 1 Jan 0.9902
## 2 Feb 0.9963
## 3 Mar 0.9744
## 4 Apr 0.9999
## 5 Mei 1.0167
## 6 Jun 1.0449
## 7 Jul 1.0626
## 8 Agu 1.0218
## 9 Sep 0.9836
## 10 Okt 0.9657
## 11 Nov 0.9556
## 12 Des 0.9883
Pada Juli indeks musiman sebesar 1,0626 menunjukkan produksi sekitar 6,26% lebih tinggi dari trend. Pada November menjadi 0,9556 yang menunjukkan produksi sekitar 4,44% lebih rendah dari trend. Selanjutnya, Pada Mei–Agustus cenderung berada di atas trend dan September–Maret cenderung berada di bawah trend.
# Residual Multipicative
dekom_mult$random
## Jan Feb Mar Apr May Jun Jul
## 2022 NA NA NA NA NA NA 1.0010461
## 2023 1.0014693 1.0005740 1.0008342 1.0004650 0.9998269 1.0005500 0.9987391
## 2024 0.9983158 0.9992111 0.9989510 0.9993202 0.9999583 0.9992351 NA
## Aug Sep Oct Nov Dec
## 2022 1.0007953 0.9982074 0.9984910 0.9980699 0.9981398
## 2023 0.9989899 1.0015778 1.0012941 1.0017152 1.0016454
## 2024 NA NA NA NA NA
summary(dekom_mult$random, na.rm = TRUE)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 0.9981 0.9989 0.9999 0.9999 1.0009 1.0017 12
sd(dekom_mult$random, na.rm = TRUE)
## [1] 0.001249183
plot(
dekom_mult$random,
main = "Residual Dekomposisi Multiplicative",
ylab = "Residual",
xlab = "Tahun"
)
abline(h = 1)
Multipicative, residual idealnya berfluktuasi disekitar 1, bukan 0.
Kesimpulan: Berdasarkan dekomposisi multiplicative, volume produksi PDAM Kabupaten Serang menunjukkan tren peningkatan dari 2022 hingga 2024 dan memiliki pola musiman tahunan. Indeks musiman tertinggi terjadi pada Juli sebesar 1,0626, yang berarti produksi sekitar 6,26% lebih tinggi dari trend. Indeks terendah terjadi pada November sebesar 0,9556, yang berarti produksi sekitar 4,44% lebih rendah dari trend.
perbandingan <- data.frame(
Bulan = bulan,
Additive = round(dekom_add$seasonal[1:12], 2),
Multiplicative = round(dekom_mult$seasonal[1:12], 4)
)
perbandingan
## Bulan Additive Multiplicative
## 1 Jan -4.54 0.9902
## 2 Feb -1.75 0.9963
## 3 Mar -11.91 0.9744
## 4 Apr -0.06 0.9999
## 5 Mei 7.84 1.0167
## 6 Jun 21.21 1.0449
## 7 Jul 27.75 1.0626
## 8 Agu 9.73 1.0218
## 9 Sep -7.35 0.9836
## 10 Okt -15.48 0.9657
## 11 Nov -20.12 0.9556
## 12 Des -5.31 0.9883
Berdasarkan hasil perbandingan, metode additive menunjukkan indeks musiman dalam satuan ribu m³, sedangkan metode multiplicative menunjukkan indeks dalam bentuk rasio terhadap trend.
Pada metode additive, bulan Juli memiliki indeks musiman tertinggi sebesar 27,75, yang berarti volume produksi pada bulan Juli cenderung 27,75 ribu m³ lebih tinggi dari nilai trendnya. Sebaliknya, bulan November memiliki indeks terendah sebesar −20,12, yang berarti produksi pada November cenderung 20,12 ribu m³ lebih rendah dari trendnya. Bulan Juni juga menunjukkan kenaikan yang cukup besar sebesar 21,21 ribu m³.
Pada metode multiplicative, bulan Juli juga memiliki indeks tertinggi sebesar 1,0626. Nilai tersebut menunjukkan bahwa produksi pada bulan Juli cenderung 6,26% lebih tinggi dari trendnya. Sementara itu, bulan November memiliki indeks terendah sebesar 0,9556, yang berarti produksi pada bulan November cenderung 4,44% lebih rendah dari trendnya.
Menariknya, kedua metode memberikan pola musiman yang sama. Bulan-bulan sekitar Mei–Agustus cenderung memiliki produksi di atas trend, sedangkan September–Maret cenderung berada di bawah trend. Puncak musiman terjadi pada Juli, sedangkan titik terendah terjadi pada November.