Я работаю с несколькими временными рядами и их характеристиками в таблице data.table в длинном формате, и я хотел бы построить несколько различных типов объекта xts
, используя синтаксис data.table
. Вот что я имею в виду:
Некоторые предварительные данные
library(data.table)
set.seed(1)
DT <- data.table(
dat = as.Date("2013-01-01") + rep(1:5, 2),
a = c(1, -1),
x = i <- rnorm(10),
y = 2 * i + rnorm(10),
z = 3 * i + rnorm(10))
Давайте (попробуем) построить 3 разных временных ряда
DT[a == 1,{
x.ts <- xts(x, order.by = dat)
y.ts <- xts(y, order.by = dat)
}]
DT[
, list(sz = mead(z/x)), by = dat][ # operate on columns of the datatable
, z.ts <- xts(sz, order.by = dat)] # construct xts
У этого кода две проблемы:
- он выводит на экран
y.ts
иz.ts
; - объекты создаются только в среде
j
, а не в глобальной среде.
Чтобы решить (2), можно назвать объекты в глобальной среде и использовать <<-
:
x.ts <- y.ts <- z.ts <- NA
DT[a == 1,{
x.ts <<- xts(x, order.by = dat)
y.ts <<- xts(y, order.by = dat)
}
]
DT[
, list(sz = mean(z/x)), by = dat][
, z.ts <<- xts(sz, order.by = dat)]
И действительно
> z.ts
[,1]
2013-01-02 2.300730
2013-01-03 4.969685
2013-01-04 1.959377
2013-01-05 1.961270
2013-01-06 3.256254
Как я могу сделать этот тип присваивания в таблице данных бесшумным? Оправдано ли такое использование <<-
или есть лучший способ сделать это?