불규칙한 시계열 M-D-Y hh : mm : ss를 일반 TS로 채우기 NA

Question

이전 게시물을 읽었지만 원하는 것을 얻을 수 없습니다. 나는 16 일 간격으로 세리를 얻을 필요가있다. (최소한 첫 번째와 마지막 날,이 경우 첫 번째/마지막 관측과 함께 간격을 시작/끝내야한다.) 나는 관측 된 변수가 상응하는 inteval과 NA에 위치하도록하고 싶다.

mdyhms     Ya Yb 
Mar-27-2009 19:56:47 25 58.25 
Mar-27-2009 20:38:59 9 81.25 
Mar-28-2009 08:00:30 9 88.75 
Mar-28-2009 09:26:29 0 89.25 
Mar-28-2009 11:57:01 8.5 74.25 
Mar-28-2009 12:19:10 7.5 71.00 
Mar-28-2009 14:17:05 1.5 70.00 
Mar-28-2009 15:13:14 NA NA 
Mar-28-2009 17:09:53 4 85.50 
Mar-28-2009 18:37:24 0 86.00 
Mar-28-2009 19:19:23 0 50.50 
Mar-28-2009 20:45:50 0 36.25 
Mar-29-2009 08:44:16 4.5 34.50 
Mar-29-2009 10:35:12 8.5 39.50 
Mar-29-2009 11:09:13 3.67 69.00 
Mar-29-2009 12:40:07 0 54.25 
Mar-29-2009 14:31:48 5.33 35.75 
Mar-29-2009 16:19:27 6.33 71.75 
Mar-29-2009 16:43:20 7.5 64.75 
Mar-29-2009 18:37:42 8 83.75 
Mar-29-2009 20:01:26 6.17 93.75 
Mar-29-2009 20:43:53 NA NA 
Mar-30-2009 08:42:05 12.67 88.50 
Mar-30-2009 09:52:57 4.33 75.50 
Mar-30-2009 12:01:32 1.83 70.75 
Mar-30-2009 12:19:40 NA NA 
Mar-30-2009 14:23:37 3.83 86.75 
Mar-30-2009 16:00:59 37.33 80.25 
Mar-30-2009 17:19:28 10.17 77.75 
Mar-30-2009 17:49:12 9.83 73.00 
Mar-30-2009 20:06:00 11.17 76.75 
Mar-30-2009 21:40:35 20.33 68.25 
Mar-31-2009 08:11:12 18.33 69.75 
Mar-31-2009 09:51:29 14.5 65.50 
Mar-31-2009 11:10:41 NA NA 
Mar-31-2009 13:27:09 NA NA 
Mar-31-2009 13:44:35 NA NA 
Mar-31-2009 16:01:23 NA NA 
Mar-31-2009 16:56:14 NA NA 
Mar-31-2009 18:27:28 NA NA 
Mar-31-2009 19:17:46 NA NA 
Mar-31-2009 21:12:22 NA NA 
Apr-01-2009 08:35:24 2.33 60.25 
Apr-01-2009 09:24:49 1.33 71.50 
Apr-01-2009 11:28:34 5.67 62.00 
Apr-01-2009 13:31:48 NA NA 
Apr-01-2009 14:52:18 NA NA 
Apr-01-2009 15:11:44 1.5 71.50 
Apr-01-2009 17:00:53 3.17 84.00 

감사합니다 [아와 YB는 관찰 된 변수] : 다음과 같이

내 데이터는 봐!

Answer 1

데이터 프레임을 "데이터"라고 가정하면 xts package을 사용합니다. 작업하기가 훨씬 쉽습니다.

#Conversion of dates 
Data$time <- as.POSIXct(Data$mdyhms,format="%b-%d-%Y %H:%M:%S") 

#conversion to time series 
library(xts) 
TimeSeries <- xts(Data[,c("Ya","Yb")],Data[,"time"]) 

그런 다음 TimeSeries를 연속적으로 사용할 수 있습니다. 정기적 인 시계열이 없으므로 표준 TS를 사용할 수 없습니다. 지구상에서 어떤 방식으로도 당신의 관찰 사이의 시간 간격은 같다고 당신이 변호 할 수는 없습니다.

편집 :

코멘트에서 발언 관련하여

, 다음 시도 할 수 있습니다 :

#Calculate the period they're into 
#This is based on GMT and the fact that POSIXct gives the number of seconds 
#passed since the origin. 5400 is 1/16 of 86400 seconds in a day 

Data$mdyhms <- as.POSIXct(Data$mdyhms,format="%b-%d-%Y %H:%M:%S",tz="GMT") 
Data$Period <- as.numeric(Data$mdyhms) %/% 5400 * 5400 

#Make a new data frame with all periods in the range of the dataframe 

Date <- as.numeric(trunc(Data$mdyhms,"day")) 
nData <- data.frame(
    Period = seq(min(Date),max(Date)+86399,by=5400) 
) 
# Merge both dataframes and take the mean of values within a dataframe 

nData <- merge(Data[c('Ya','Yb','Period')],nData,by="Period",all=T) 
nData <- ddply(nData,"Period",mean,na.rm=T) 

#Make the time series and get rid of the NaN values 
#These come from averaging vectors with only NA 
TS <- ts(nData[c('Ya','Yb')],frequency=16) 
TS[is.nan(TS)] <- NA