library(qgam)
data(AUDem)
# Mean demand over the period
plot(AUDem$meanDem$dem, type = 'l')
# 20 quantiles of individual demand over 5 days
matplot(seq(0.01, 0.99, length.out = 20),
t(AUDem$qDem48[c(1, 50, 75, 100, 250), ]),
type = 'l',
ylab = "Electricity demand (KW)",
xlab = expression("Probability level " * "(p)"),
lty = 1)
Run the code above in your browser using DataLab