R: использование ddply в цикле над столбцами фрейма данных

Question

Мне нужно рассчитать и добавить к кадру данных несколько новых столбцов на основе значений в каждом столбце в поднаборе столбцов в кадре данных. В этих столбцах хранятся данные временных рядов (имеется общий столбец даты). Например, мне нужно рассчитать изменение за тот же месяц в предыдущем году для дюжины столбцов. Я мог бы указать их и рассчитать их отдельно, но это становится обременительным с большим количеством столбцов для преобразования, поэтому я пытаюсь автоматизировать процесс с циклом for.

Я делал все ОК, пока не попытался использовать ddply, чтобы создать столбец для текущей общей суммы за год. Что происходит, так это то, что ddply добавляет новые строки во время каждой итерации через цикл и включает в себя эти новые строки в вычислении cumsum. У меня есть два вопроса.

В. Как я могу получить ddply для вычисления правильной cumsum? В. Как я могу указать имя столбца во время вызова ddply, вместо использования фиктивного значения и переименования потом?

[Edit: я говорил слишком рано, обновленный код ниже не работает на данный момент, просто FYI]

require(lubridate) 
require(plyr) 
require(xts) 

set.seed(12345) 
# create dummy time series data 
monthsback <- 24 
startdate <- as.Date(paste(year(now()),month(now()),"1",sep = "-")) - months(monthsback) 
mydf <- data.frame(mydate = seq(as.Date(startdate), by = "month", length.out = monthsback), 
        myvalue1 = runif(monthsback, min = 600, max = 800), 
        myvalue2 = runif(monthsback, min = 200, max = 300)) 

mydf$year <- as.numeric(format(as.Date(mydf$mydate), format="%Y")) 
mydf$month <- as.numeric(format(as.Date(mydf$mydate), format="%m")) 
newcolnames <- c('myvalue1','myvalue2') 

for (i in seq_along(newcolnames)) { 
    print(newcolnames[i]) 
    mydf$myxts <- xts(mydf[, newcolnames[i]], order.by = mydf$mydate) 
    ## Calculate change over same month in previous year 
    mylag <- 12 
    mydf[, paste(newcolnames[i], "_yoy", sep = "", collapse = "")] <- as.numeric(diff(mydf$myxts, lag = mylag)/ lag(mydf$myxts, mylag)) 
    ## Calculate change over previous month 
    mylag <- 1 
    mydf[, paste(newcolnames[i], "_mom", sep = "", collapse = "")] <- as.numeric(diff(mydf$myxts, lag = mylag)/ lag(mydf$myxts, mylag)) 

    ## Calculate cumulative figure 
    #mydf$newcol <- as.numeric(mydf$myxts) 
    mydf$newcol <- 1 
    mydf <- ddply(mydf, .(year), transform, newcol = cumsum(as.numeric(mydf$myxts))) 
    colnames(mydf)[colnames(mydf)=="newcol"] <- paste(newcolnames[i], "_cuml", sep = "", collapse = "") 

} 

mydf 

Answer 1

В вашем цикле, так как myxts не является частью кадра данных, не разделились в заявлении ddply вместе со всем остальным. Измените его на:

mydf$myxts <- xts(mydf[, newcolnames[i]], order.by = mydf$mydate) 

Я не знаю ни одного способа, чтобы использовать динамически генерируемые имена с transform.