2014-12-30 2 views
1

내 측정 값의 원시 값을 포함하는 데이터 파일이 있습니다. 상단 입력 내용은 다음과 같습니다.15 분마다 데이터 포인트가있는 시계열

495.08 
1117.728 
872.712 
665.632 
713.296 
1172.44 
1302.544 
1428.832 
1413.536 
1361.896 
1126.656 
644.776 
1251.616 
1252.824 
[...] 

측정은 15 분 간격으로 수행됩니다. 이것을 R 시계열 객체에 넣기 위해 다음 코드를 사용합니다. 데이터가 일년에서만이지만

R output

:

data <- scan("./data__2.dat",skip=1) 
datats <- ts(data, frequency=24*60/15, start=c(2014,1)) 

는하지만이 날 떠난다. 따라서 빈도가 잘못된 것처럼 보입니다. 이 문제를 해결하는 방법에 대한 의견이 있으십니까? 이렇게함으로써

+0

내 생각 엔 당신이 주파수 내가'을 의심하지만 60 * 24 * 15분의 365 = 35040. ts'가 대처할 수 있어야한다이다 365보다 높은 빈도 ... –

+0

@ 데이비드 : 이것은 내 첫 추측이었습니다 ... 어떤 제안이 문제를 해결하는 방법? – Coxer

+0

는'xts'로 변환합니까? – Khashaa

답변

3

:

library(xts) 
library(lubridate) 

df <- data.frame(interval = seq(ymd_hms('2014-01-21 00:00:00'), 
           by = '15 min',length.out=(60*24*365/15)), 
       data = rnorm(60*24*365/15)) 
ts <- xts(df, order.by=df$interval) 

를 당신이 얻을 :

   interval  data 
1 2014-01-21 00:00:00 -1.3975823 
2 2014-01-21 00:15:00 -0.4710713 
3 2014-01-21 00:30:00 0.9149273 
4 2014-01-21 00:45:00 -0.3053136 
5 2014-01-21 01:00:00 -1.2459707 
6 2014-01-21 01:15:00 0.4749215