Sai phân theo mùa để đạt tính dừng
Với dữ liệu có mùa vụ, người ta thường lấy sai phân giữa các quan sát ở cùng một mùa của các năm liên tiếp, thay vì giữa các kỳ liên tiếp. Ví dụ, với dữ liệu theo quý, bạn sẽ lấy chênh lệch giữa Q1 của một năm và Q1 của năm trước đó. Cách này gọi là sai phân theo mùa.
Đôi khi bạn cần áp dụng cả sai phân theo mùa và sai phân trễ 1 cho cùng một chuỗi, tức là tính sai phân của chính sai phân.
Trong bài tập này, bạn sẽ kết hợp sai phân và biến đổi để làm cho một chuỗi thời gian trông có tính dừng. Bộ dữ liệu ở đây là h02, chứa 17 năm doanh số thuốc corticosteroid theo tháng tại Australia. Dữ liệu đã được nạp vào không gian làm việc của bạn.
Bài tập này là một phần của khóa học
Dự báo bằng R
Hướng dẫn bài tập
- Vẽ biểu đồ dữ liệu để quan sát xu hướng và mùa vụ.
- Lấy
log()của dữ liệuh02rồi áp dụng sai phân theo mùa bằng cách dùng giá trịlagphù hợp trongdiff(). Gán kết quả vàodifflogh02. - Vẽ biểu đồ dữ liệu sau khi đã log và sai phân.
- Vì
difflogh02vẫn có vẻ chưa dừng, hãy lấy thêm một sai phân trễ 1 bằng cách áp dụngdiff()lên chính nó và lưu vàoddifflogh02. Vẽ chuỗi kết quả. - Vẽ ACF của chuỗi cuối cùng
ddifflogh02bằng hàm thích hợp.
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# Plot the data
___
# Take logs and seasonal differences of h02
difflogh02 <- diff(log(___), lag = ___)
# Plot difflogh02
___
# Take another difference and plot
ddifflogh02 <- ___
___
# Plot ACF of ddifflogh02
___