MulaiMulai sekarang secara gratis

Peramalan permintaan listrik

Anda juga dapat memodelkan permintaan listrik harian sebagai fungsi dari suhu. Seperti yang mungkin Anda lihat pada tagihan listrik, listrik lebih banyak digunakan pada hari-hari panas karena pendingin udara dan pada hari-hari dingin karena pemanas.

Dalam latihan ini, Anda akan menyesuaikan model regresi kuadrat dengan galat ARMA. Satu tahun data harian disimpan sebagai elecdaily yang mencakup total permintaan harian, variabel indikator untuk hari kerja (hari kerja direpresentasikan dengan 1, dan hari non-kerja direpresentasikan dengan 0), serta suhu maksimum harian. Karena terdapat musiman mingguan, frequency telah disetel ke 7.

Mari kita lihat tiga baris pertama:

> elecdaily[1:3, ]
       Demand Temperature Workday
[1,] 174.8963        26.0       0
[2,] 188.5909        23.0       1
[3,] 188.9169        22.2       1

elecdaily telah dimuat sebelumnya ke ruang kerja Anda.

Latihan ini adalah bagian dari kursus

Peramalan di R

Lihat Kursus

Petunjuk latihan

  • Hasilkan plot runtun waktu untuk hanya permintaan harian dan suhu maksimum dengan facetting.
  • Siapkan matriks peubah penjelas yang mencakup MaxTemp untuk suhu maksimum, MaxTempSq yang merepresentasikan nilai kuadrat dari suhu maksimum, dan Workday, dalam urutan tersebut. Jelas bahwa argumen kedua cbind() memerlukan operator matematika sederhana.
  • Sesuaikan model regresi dinamis pada kolom permintaan dengan galat ARIMA dan beri nama fit.
  • Jika hari berikutnya adalah hari kerja (indikator 1) dengan prakiraan suhu maksimum 20°C, berapa prakiraan permintaannya? Isi nilai yang sesuai dalam cbind() untuk argumen xreg pada forecast().

Latihan interaktif praktis

Cobalah latihan ini dengan menyelesaikan kode contoh berikut.

# Time plots of demand and temperatures
autoplot(elecdaily[, c(___, ___)], facets = ___)

# Matrix of regressors
xreg <- cbind(MaxTemp = elecdaily[, "Temperature"],
              MaxTempSq = ___,
              Workday = ___)

# Fit model
fit <- auto.arima(___, xreg = xreg)

# Forecast fit one day ahead
forecast(___, xreg = cbind(___, ___, ___))
Edit dan Jalankan Kode