2. 지수평활화

1. 단순지수 평활화

 

시계열용.xlsx

0.01MB

R코드

 

library(readxl)
library(forecast)

# 1. 데이터 불러오기
data <- read_excel("d:/mydata/시계열용.xlsx")

# 2. 2003년부터 2009년까지의 wRC+ 추출 (7기)
wrc_data <- data[data$Time >= 2003 & data$Time <= 2009, "wRC+"]
wrc_ts <- ts(wrc_data$`wRC+`)

# 3. 단순지수평활 모형 적합
model <- ses(wrc_ts, h = 1)  # h = 1: 1기 예측

# 4. 예측 결과 출력
forecast_value <- model$mean
print(paste("2010년(10기) 예측 wRC+:", round(forecast_value, 2)))

# 시각화 (옵션)
plot(model, main = "단순지수평활에 의한 이치로 wRC+ 예측")

 

2. 이중지수 평활화

<여기서 예측식은 중요하지는 않다, 왜 저리 도출되는 지 알려하지 말자>

 

R코드

# 1. 데이터 불러오기
data <- read_excel("d:/mydata/시계열용.xlsx")

# 2. 2003~2009년까지 wRC+만 추출 (7기)
wrc <- data[data$Time >= 2003 & data$Time <= 2009, "wRC+"]
wrc_ts <- ts(wrc$`wRC+`, start = 2003)

# 3. 이중지수평활 (Holt's Linear Method) 적용
holt_model <- holt(wrc_ts, h = 1)

# 4. 예측 결과 출력
print(paste("예측된 2010년 wRC+:", round(holt_model$mean, 2)))

# 5. 시각화
plot(holt_model, main = "Holt 이중지수평활화: 2010년 wRC+ 예측", ylab = "wRC+", xlab = "Year")
points(2010, holt_model$mean, col = "red", pch = 19)
text(2010, holt_model$mean + 2, labels = round(holt_model$mean, 1), col = "red")

 

 

 

Cf) 지수평활화에서 보통 S0 = y1(실제값)으로 설정한다