ประมาณค่าโมเดล White Noise
สำหรับอนุกรมเวลา y ที่กำหนด เราสามารถ fit โมเดล white noise (WN) ได้โดยใช้ฟังก์ชัน arima(..., order = c(0, 0, 0)) โดยโมเดล WN คือโมเดล ARIMA(0,0,0) การเรียกใช้ฟังก์ชัน arima() จะคืนค่าข้อมูลเกี่ยวกับโมเดลที่ประมาณได้ ซึ่งสำหรับโมเดล WN จะประกอบด้วยค่าเฉลี่ยที่ประมาณได้ (แสดงเป็น intercept) และค่าความแปรปรวนที่ประมาณได้ (แสดงเป็น sigma^2)
ในแบบฝึกหัดนี้ จะได้สำรวจคุณสมบัติของโมเดล WN ค่าเฉลี่ยที่ประมาณได้คือเท่าไร ลองเปรียบเทียบกับค่าเฉลี่ยของตัวอย่างโดยใช้ฟังก์ชัน mean() และค่าความแปรปรวนที่ประมาณได้คือเท่าไร ลองเปรียบเทียบกับค่าความแปรปรวนของตัวอย่างโดยใช้ฟังก์ชัน var()
อนุกรมเวลา y ถูกโหลดไว้แล้ว และแสดงอยู่ในกราฟด้านข้าง
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การวิเคราะห์อนุกรมเวลาใน R
คำแนะนำการฝึกหัด
- ใช้
arima()เพื่อประมาณค่าโมเดล WN สำหรับyโดยอย่าลืมระบุอาร์กิวเมนต์order = c(0, 0, 0)หลังจากระบุข้อมูลแล้ว - คำนวณค่าเฉลี่ยและค่าความแปรปรวนของ
yโดยใช้mean()และvar()ตามลำดับ จากนั้นเปรียบเทียบผลลัพธ์กับผลลัพธ์ที่ได้จากคำสั่งarima()
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Fit the WN model to y using the arima command
# Calculate the sample mean and sample variance of y