Salin nilai terakhir selama periode harian

Saya memiliki objek XTS multi-hari, dan saya mencoba membuat indikator yang jika benar, akan tetap benar sepanjang hari. Pendekatan yang saya coba (tetapi tidak berhasil) adalah menggabungkan fungsi na.locf dengan apply daily:

output <- apply.daily(x, na.locf)

Kode yang dapat direproduksi:

y <- as.xts(c(NA,NA,1,NA,NA,NA,NA,NA,NA),as.POSIXct(c(
                                   "2010-01-05 00:00:00", "2010-01-05 00:04:00", "2010-01-05 00:08:00", 
                                   "2010-01-05 00:12:00", "2010-01-05 00:16:00", "2010-01-05 00:20:00",
                                   "2010-01-06 00:00:00", "2010-01-06 00:04:00", "2010-01-06 00:08:00")))

Output yang diinginkan adalah menyalin '1' untuk sisa hari itu. Jadi:

y <- as.xts(c(NA,NA,1,1,1,1,NA,NA,NA),as.POSIXct(c(
                                   "2010-01-05 00:00:00", "2010-01-05 00:04:00", "2010-01-05 00:08:00", 
                                   "2010-01-05 00:12:00", "2010-01-05 00:16:00", "2010-01-05 00:20:00",
                                   "2010-01-06 00:00:00", "2010-01-06 00:04:00", "2010-01-06 00:08:00")))

r xts
person Ed Wilson    schedule 12.05.2016    source sumber

Jawaban (2)


arrow_upward
2
arrow_downward

Salah satu pilihannya adalah

y1 <- ave(y, as.Date(index(y)), FUN= function(x) na.locf(x, na.rm=FALSE))
y1
#                      [,1]
#2010-01-05 00:00:00   NA
#2010-01-05 00:04:00   NA
#2010-01-05 00:08:00    1
#2010-01-05 00:12:00    1
#2010-01-05 00:16:00    1
#2010-01-05 00:20:00    1
#2010-01-06 00:00:00   NA
#2010-01-06 00:04:00   NA
#2010-01-06 00:08:00   NA

str(y1)
# An ‘xts’ object on 2010-01-05/2010-01-06 00:08:00 containing:
#  Data: num [1:9, 1] NA NA 1 1 1 1 NA NA NA
#  Indexed by objects of class: [POSIXct,POSIXt] TZ: 
#  Original class: 'double'  
# xts Attributes:  
# NULL

str(y)
#An ‘xts’ object on 2010-01-05/2010-01-06 00:08:00 containing:
#  Data: num [1:9, 1] NA NA 1 NA NA NA NA NA NA
#  Indexed by objects of class: [POSIXct,POSIXt] TZ: 
#  Original class: 'double'  
#  xts Attributes:  
# NULL
person akrun    schedule 12.05.2016
comment
Ini berhasil dengan baik. Apakah Anda menggunakan fungsi ave hanya untuk membuat subset? Atau ada alasan lain? - person Ed Wilson; 12.05.2016
comment
@EdWilson Fungsi ave digunakan untuk menerapkan na.locf pada setiap 'Tanggal' - person akrun; 13.05.2016

arrow_upward
2
arrow_downward

Saya pikir period.apply dipanggil oleh apply.daily tidak menyukai nilai yang dikembalikan oleh na.locf. Belum menyelidiki secara menyeluruh alasannya. Bagaimanapun, saya mencoba melakukan apa yang seharusnya dilakukan secara berputar-putar. Saya dapat melihat bahwa jawaban Akrun jelas lebih unggul dari ini. Tinggalkan saja ini di sini.

R> y <- as.xts(c(NA,NA,1,NA,NA,NA,NA,NA,NA),
+             as.POSIXct(c("2010-01-05 00:00:00", "2010-01-05 00:04:00",
+                          "2010-01-05 00:08:00", "2010-01-05 00:12:00",
+                          "2010-01-05 00:16:00", "2010-01-05 00:20:00",
+                          "2010-01-06 00:00:00", "2010-01-06 00:04:00",
+                          "2010-01-06 00:08:00")))
R> endpoints(y, "days")
[1] 0 6 9
R> ep <- endpoints(y, "days")
R> diff(ep)
[1] 6 3
R> dep <- diff(ep)
R> rep.int(1:length(dep), times=dep)
[1] 1 1 1 1 1 1 2 2 2
R> runs <- rep.int(1:length(dep), times=dep)
R> lapply(split(y, runs), na.locf, na.rm=FALSE)
$`1`

2010-01-05 00:00:00 NA
2010-01-05 00:04:00 NA
2010-01-05 00:08:00  1
2010-01-05 00:12:00  1
2010-01-05 00:16:00  1
2010-01-05 00:20:00  1

$`2`

2010-01-06 00:00:00 NA
2010-01-06 00:04:00 NA
2010-01-06 00:08:00 NA

R> splits <- lapply(split(y, runs), na.locf, na.rm=FALSE)
R> do.call('rbind', splits)

2010-01-05 00:00:00 NA
2010-01-05 00:04:00 NA
2010-01-05 00:08:00  1
2010-01-05 00:12:00  1
2010-01-05 00:16:00  1
2010-01-05 00:20:00  1
2010-01-06 00:00:00 NA
2010-01-06 00:04:00 NA
2010-01-06 00:08:00 NA
R> ynew <- do.call('rbind', splits)
person asb    schedule 12.05.2016
comment
Pendekatan ini masuk akal dan merupakan pendekatan xts yang khas. Meskipun bisa mencapai hasil penuh hanya dengan menggunakan do.call(rbind, lapply(split(y, "days"), na.locf)) - person FXQuantTrader; 27.10.2016