Блог

У меня есть фрейм данных с несколькими столбцами:

• состояние
• округ
• год

Тогда x, y и z, где x, y и z — наблюдения, уникальные для триплета, перечисленного выше. Я ищу разумный способ сохранить это во временных рядах, и xts мне этого не позволит, поскольку для каждого временного индекса существует несколько наблюдений. Я просмотрел пакет hts, но у меня возникли проблемы с пониманием того, как перенести в него мои данные из dataframe.

(Да, я опубликовал тот же вопрос на Quora, и мне посоветовали перенести его сюда!)

Одним из вариантов является изменение формы ваших данных, чтобы у вас был столбец для каждой комбинации штата и округа. Это позволяет вам построить матрицу xts :

 require(reshape)
Opt1 <- as.data.frame(cast(Data, Date ~ county   State, value="Val"))
rownames(Opt1) <- Opt1$Date
Opt1$Date <- NULL
as.xts(Opt1)
  

В качестве альтернативы вы могли бы работать со списком объектов xts, каждый раз проверяя, что у вас правильный формат, запрошенный xts. То же самое относится к любому из других пакетов timeseries. Возможным решением было бы :

 Opt2 <- 
  with(Data,
    by(Data,list(county,State,year),
      function(x){
        rownames(x) <- x$Date
        x <- x["Val"]
        as.xts(x)
      }
    )
  )
  

Что позволило бы сделать что-то вроде :

 Opt2[["d","b","2012"]]
  

чтобы выбрать конкретный временной ряд. Для этого можно использовать все опции xts. Вы можете перебирать округа, штаты и годы для построения графиков, подобных этому :

Код для построения графика :

 counties <- dimnames(Opt2)[[1]]
states <- dimnames(Opt2)[[2]]
years <- dimnames(Opt2)[[3]]

op <- par(mfrow=c(3,6))
apply(
  expand.grid(counties,states,years),1,
  function(i){
    plot(Opt2[[i[1],i[2],i[3]]],main=paste(i,collapse="-"))
    invisible()
  }
)
par(op)
  

Тестовые данные :

 Data <- data.frame( State = rep(letters[1:3],each=90),
            county = rep(letters[4:6],90),
            Date = rep(seq(as.Date("2011-01-01"),by="month",length.out=30),each=3),
            Val = runif(270)
)
Data$year <- as.POSIXlt(Data$Date)$year   1900