2016-06-14 2 views
1

제 계절 시계열 데이터 주위에 TBATS 모델을 장착하고 forecast 패키지를 사용하여 예측을 얻었습니다. 내 R 코드는 다음과 같습니다R의 'forecast'라이브러리를 사용하여 예측 값을 CSV 파일에 저장하는 방법은 무엇입니까?

library("forecast") 
data = read.csv("data.csv") 
season_info <- msts(data,seasonal.periods=c(24,168)) 
model <- tbats(season_info) 
forecasted <- forecast.tbats(best_model,h=24,level=90) 
forecasted 

지금, 나는 같은 출력이라는 변수가 '예상'이 있습니다

> forecasted 
     Point Forecast Lo 90 Hi 90 
6.940476  5080.641 4734.760 5426.523 
6.946429  5024.803 4550.111 5499.496 
6.952381  4697.625 4156.516 5238.733 
6.958333  4419.105 3832.765 5005.446 
6.964286  4262.782 3643.528 4882.037 
6.970238  4187.629 3543.062 4832.196 
6.976190  4349.196 3684.444 5013.947 
6.982143  4484.108 3802.574 5165.642 
6.988095  4247.858 3551.955 4943.761 
6.994048  3851.379 3142.831 4559.927 
7.000000  3575.951 2855.962 4295.941 
7.005952  3494.943 2764.438 4225.449 
7.011905  3501.354 2760.968 4241.739 
7.017857  3445.563 2695.781 4195.345 

내가 열 '예측'에서 예측 된 값을 수집하고 저장해야 CSV 파일에 저장합니다. TBATS의 페이지와 '예측'방법을 온라인으로 읽으려고했지만 'Hi' 'Lo'와 'Point'와 같은 다른 열은 무시하고 예상 값의 특정 열을 추출 할 수있는 방법을 언급하지 않았습니다.

내 CSV에서이 출력을 찾고 있어요 :

hour,forecasted_value 
0,5080.641 
1,5024.803 
2,4697.625 
... 

답변

2

그들은 세 부분으로 저장됩니다. 당신은 str(ret)와 객체 구조를 볼 수 있습니다 :

library(forecast) 
fit <- tbats(USAccDeaths) 
ret <- forecast(fit) 
ret$upper # Upper interval 
ret$lower # Lower interval 
ret$mean # Point forecast 
+0

뜻을 의미한다. 0-24의 숫자를 생성하고 CSV에 저장하는 방법에 대한 아이디어가 있습니까? 예를 들어 위 질문에서 언급 한 출력 CSV 형식입니다. – learnerX

+0

정확히 24 개의 평균값이 있습니까? –

+0

예,'h = 24'이 줄에서 :'forecasted <- forecast.tbats (best_model, h = 24, level = 90)'은 정확히 24 가지 예측을 보장합니다. – learnerX

2
당신은 print()를 사용하여 표시되는 출력을 얻을 수 있습니다

:

library("forecast") 
data = read.csv("data.csv") 
season_info <- msts(data,seasonal.periods=c(24,168)) 
model <- tbats(season_info) 
forecasted <- forecast.tbats(best_model,h=24,level=90) 
dfForec <- print(forecasted) 

이 당신에게 data.frame을 줄 것이다, 지금 당신은 당신이 원하는 열을 선택할 수 있습니다 , 즉. 포인트 예측에 대해서만 dfForec[, 1]을 입력 한 다음 write.csv(dfForec[, 1, drop = FALSE], ...)을 사용하여 플랫 파일에 씁니다. MSTS (데이터, 계절 -

0

사용하여 포인트

라이브러리 ("예측")

데이터 = read.csv ("data.csv") < season_info

을 예측 얻기를위한 기능을 의미 .periods에서의 C = (24,168))

모델 < - tbats (season_info)

가 <을 예측 - (forecast.tbats은 (best_model, H = 24, 수준 = 90)) $

또는

은 $ 예상 당신은 작동이 맞아,

관련 문제