[英]Getting data from time series as per start and end time
我有100個文本文件,其中包含在不同時間點開始和結束的時間序列。 我想提取系列中常見時間段的值。 使用以下代碼生成示例數據:
set.seed(1)
D1 = data.frame(time = seq(ISOdatetime(2012, 6, 26, 3, 15, 00),
length = 500, by = 900),
value = rnorm(500))
D2 = data.frame(time = seq(ISOdatetime(2012, 6, 24, 5, 30, 00),
length = 541, by = 900),
value = rnorm(541))
D3 = data.frame(time = seq(ISOdatetime(2012, 6, 23, 5, 45, 00),
length = 700, by = 900),
value = rnorm(700))
這些數據將為您提供3個時間序列的開始和結束以及不同的時間。 我希望只保留常用時間段的值並刪除其余時間。 即如果,
然后我希望保留三個時間序列的交集數據即數據對應: -
我搜索了SO和其他網站,但沒有找到任何解決方案。 需要幫助。 我如何實現這一目標?
您似乎需要熟悉xts
包。 將數據框轉換為xts
時間序列對象並使用merge
。 merge
將合並所有值,因此如果您希望所有值都出現,您也可以使用na.omit
。
require(xts)
D1 = xts(d1$Value, d1$Time)
D2 = xts(d2$Value, d2$Time)
D3 = xts(d3$Value, d3$Time)
temp = merge(D1, D2, D3)
這是一些示例輸出。 對於head
和tail
,請注意NA
值的存在。
head(temp)
# D1 D2 D3
# 2012-06-26 13:15:19 -0.50219235 NA NA
# 2012-06-26 13:30:19 0.13153117 NA NA
# 2012-06-26 13:45:19 -0.07891709 NA NA
# 2012-06-26 14:00:19 0.88678481 NA NA
# 2012-06-26 14:15:19 0.11697127 NA NA
# 2012-06-26 14:30:19 0.31863009 NA NA
tail(temp)
# D1 D2 D3
# 2012-07-04 05:45:19 NA NA 1.4799645
# 2012-07-04 06:00:19 NA NA -0.3942801
# 2012-07-04 06:15:19 NA NA -0.6767234
# 2012-07-04 06:30:19 NA NA -0.2425192
# 2012-07-04 06:45:19 NA NA 0.4547177
# 2012-07-04 07:00:19 NA NA 1.1712661
head(na.omit(temp))
# D1 D2 D3
# 2012-06-27 14:15:19 -0.3329234 -1.63230970 0.75619287
# 2012-06-27 14:30:19 1.3631137 -0.06299626 -1.36131851
# 2012-06-27 14:45:19 -0.4691473 -0.70544686 -0.60876462
# 2012-06-27 15:00:19 0.8428756 -0.31417818 -0.21174696
# 2012-06-27 15:15:19 -1.4579937 -0.26694627 -0.67847242
# 2012-06-27 15:30:19 -0.4003059 0.15315947 0.06665787
tail(na.omit(temp))
# D1 D2 D3
# 2012-07-01 16:45:19 -0.49419020 1.1911322 2.73143169
# 2012-07-01 17:00:19 -1.71111303 0.7613245 0.57057667
# 2012-07-01 17:15:19 0.04005805 -0.1210687 1.32083870
# 2012-07-01 17:30:19 -0.56114348 -1.2250590 0.09951626
# 2012-07-01 17:45:19 -2.55736206 -0.1637461 -0.39435301
# 2012-07-01 18:00:19 -0.69677881 -1.3138963 0.63649492
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.