Как преобразовать вложенный список

У меня есть вложенный список и мне нужно вычислить корреляцию frt и srt

   $`bs. bs`
    fapp  frt sapp  srt
1   bs 2280   bs 0.25
2   bs 2287   bs 0.25
3   bs 2288   bs 0.25
4   bs 2289   bs 0.25

$`bs. lhc`
    fapp  frt sapp      srt
5   bs 2320  lhc 0.250000
6   bs 2333  lhc 0.250214
7   bs 2524  lhc 0.316449

источник:

    structure(list(`bs. bs` = structure(list(fapp = structure(c(1L,
    1L, 1L, 1L), .Label = "bs", class = "factor"), frt = c(2280L,
    2287L, 2288L, 2289L), sapp = structure(c(1L, 1L, 1L, 1L), .Label = c(" bs",
    " lhc"), class = "factor"), srt = c(0.25, 0.25, 0.25, 0.25)), .Names = c("fapp",
    "frt", "sapp", "srt"), row.names = c(NA, 4L), class = "data.frame"),
   `bs. lhc` = structure(list(fapp = structure(c(1L, 1L, 1L), .Label = "bs", class=         "factor"),
    frt = c(2320L, 2333L, 2524L), sapp = structure(c(2L,
    2L, 2L), .Label = c(" bs", " lhc"), class = "factor"),
    srt = c(0.25, 0.250214, 0.316449)), .Names = c("fapp",
"frt", "sapp", "srt"), row.names = 5:7, class = "data.frame")), .Names = c("bs. bs",
   "bs. lhc"))

что-то типа

    ddply(y,.(fapp + sapp),cor)

or

    ddply(y,.(fapp,sapp),cor)

не работает


r plyr
person nairboon    schedule 15.10.2012    source источник
comment
Чтобы получить наилучший ответ, создайте воспроизводимый пример...   -  person Paul Hiemstra    schedule 15.10.2012

Ответы (3)


arrow_upward
3
arrow_downward 

> ldply(y, function(x) { x$corr <- cor(x$frt, x$srt); x })
      .id fapp  frt sapp      srt      corr
1  bs. bs   bs 2280   bs 0.250000        NA
2  bs. bs   bs 2287   bs 0.250000        NA
3  bs. bs   bs 2288   bs 0.250000        NA
4  bs. bs   bs 2289   bs 0.250000        NA
5 bs. lhc   bs 2320  lhc 0.250000 0.9985343
6 bs. lhc   bs 2333  lhc 0.250214 0.9985343
7 bs. lhc   bs 2524  lhc 0.316449 0.9985343
Warning message:
In cor(x$frt, x$srt) : the standard deviation is zero
Calls: ldply -> llply -> structure -> lapply -> FUN -> cor

или сохранить результат в виде списка

> llply(y, function(x) { x$corr <- cor(x$frt, x$srt); x })
$`bs. bs`
  fapp  frt sapp  srt corr
1   bs 2280   bs 0.25   NA
2   bs 2287   bs 0.25   NA
3   bs 2288   bs 0.25   NA
4   bs 2289   bs 0.25   NA

$`bs. lhc`
  fapp  frt sapp      srt      corr
5   bs 2320  lhc 0.250000 0.9985343
6   bs 2333  lhc 0.250214 0.9985343
7   bs 2524  lhc 0.316449 0.9985343

Warning message:
In cor(x$frt, x$srt) : the standard deviation is zero
Calls: llply -> structure -> lapply -> FUN -> cor
person Maiasaura    schedule 15.10.2012

arrow_upward
1
arrow_downward

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

person Charlie    schedule 15.10.2012
comment
Вы имели в виду ldply наверное :). Хотя получить такой же компактный синтаксис с вложенным списком несколько сложнее. - person Paul Hiemstra; 15.10.2012

arrow_upward
1
arrow_downward

ddply явно предназначен для работы с data.frame. Итак, я бы сначала поместил все элементы списка в data.frame:

dat = do.call("rbind", nested_list)

а затем используйте ddply:

ddply(dat, .(fapp, sapp), corr)
person Paul Hiemstra    schedule 15.10.2012