Объединение векторов неравной длины в кадр данных

Question

У меня есть список векторов, которые являются временными рядами неравной длины. Моя конечная цель - построить временной ряд на графике ggplot2. Полагаю, мне лучше сначала объединить векторы в кадре данных (где более короткие векторы будут расширены с помощью NA), также потому, что я хочу экспортировать данные в табличный формат, такой как .csv, для просмотра другими людьми.

У меня есть список, содержащий имена всех векторов. Хорошо, что заголовки столбцов задаются первым вектором, который самый длинный. E.g.:

> mylist
[[1]]
[1] "vector1"

[[2]]
[1] "vector2"

[[3]]
[1] "vector3"

и т.д.

Я знаю, что нужно использовать пакет plyr Хэдли, но я думаю, проблема в том, что мой список содержит имена векторов, а не сами векторы, поэтому, если я наберу:

do.call(rbind, mylist)

Я получаю df с одним столбцом, содержащим имена dfs, которые я хотел объединить.

> do.call(rbind, actives)
      [,1]           
 [1,] "vector1" 
 [2,] "vector2" 
 [3,] "vector3" 
 [4,] "vector4" 
 [5,] "vector5" 
 [6,] "vector6" 
 [7,] "vector7" 
 [8,] "vector8" 
 [9,] "vector9" 
[10,] "vector10"

и т.д.

Даже если я создаю список с самим объектом, я получаю пустой фрейм данных:

mylist <- list(vector1, vector2)
mylist
[[1]]
        1         2         3         4         5         6         7         8         9        10        11        12 
0.1875000 0.2954545 0.3295455 0.2840909 0.3011364 0.3863636 0.3863636 0.3295455 0.2954545 0.3295455 0.3238636 0.2443182 
       13        14        15        16        17        18        19        20        21        22        23        24 
0.2386364 0.2386364 0.3238636 0.2784091 0.3181818 0.3238636 0.3693182 0.3579545 0.2954545 0.3125000 0.3068182 0.3125000 
       25        26        27        28        29        30        31        32        33        34        35        36 
0.2727273 0.2897727 0.2897727 0.2727273 0.2840909 0.3352273 0.3181818 0.3181818 0.3409091 0.3465909 0.3238636 0.3125000 
       37        38        39        40        41        42        43        44        45        46        47        48 
0.3125000 0.3068182 0.2897727 0.2727273 0.2840909 0.3011364 0.3181818 0.2329545 0.3068182 0.2386364 0.2556818 0.2215909 
       49        50        51        52        53        54        55        56        57        58        59        60 
0.2784091 0.2784091 0.2613636 0.2329545 0.2443182 0.2727273 0.2784091 0.2727273 0.2556818 0.2500000 0.2159091 0.2329545 
       61 
0.2556818 

[[2]]
        1         2         3         4         5         6         7         8         9        10        11        12 
0.2824427 0.3664122 0.3053435 0.3091603 0.3435115 0.3244275 0.3320611 0.3129771 0.3091603 0.3129771 0.2519084 0.2557252 
       13        14        15        16        17        18        19        20        21        22        23        24 
0.2595420 0.2671756 0.2748092 0.2633588 0.2862595 0.3549618 0.2786260 0.2633588 0.2938931 0.2900763 0.2480916 0.2748092 
       25        26        27        28        29        30        31        32        33        34        35        36 
0.2786260 0.2862595 0.2862595 0.2709924 0.2748092 0.3396947 0.2977099 0.2977099 0.2824427 0.3053435 0.3129771 0.2977099 
       37        38        39        40        41        42        43        44        45        46        47        48 
0.3320611 0.3053435 0.2709924 0.2671756 0.2786260 0.3015267 0.2824427 0.2786260 0.2595420 0.2595420 0.2442748 0.2099237 
       49        50        51        52        53        54        55        56        57        58        59        60 
0.2022901 0.2251908 0.2099237 0.2213740 0.2213740 0.2480916 0.2366412 0.2251908 0.2442748 0.2022901 0.1793893 0.2022901 

но

do.call(rbind.fill, mylist)
data frame with 0 columns and 0 rows

Я пытался преобразовать векторы в кадры данных, но функции cbind.fill нет, поэтому plyr жалуется, что кадры данных имеют разную длину.

Итак, мои вопросы:

• Это лучший подход? Имейте в виду, что целями являются: а) график ggplot2 и б) таблица с временным рядом, которую нужно просматривать за пределами R

• Как лучше всего получить список объектов, начинающийся со списка имен этих объектов?

• Какой график лучше всего подходит для выделения шаблонов 60 временных рядов? Масштаб тот же, но я предсказываю, что будет много перепланировок. Так как это когортный анализ, может быть полезно использовать цвет, чтобы выделить разные когорты с точки зрения свежести (как непрерывную переменную). Но как избежать перелетов? Различия будут минимальными, поэтому огранка может сделать зрителя неспособным понять разницу.

Answer 1

Я думаю, что вы, возможно, подходите к этому неверно:

Если у вас есть временные ряды неравной длины, то самое лучшее, что нужно сделать, это сохранить их как временные ряды и merge их.Большинство пакетов временных рядов позволяют это.Таким образом, вы получите многовариантный временной ряд, и каждое значение будет правильно связано с одной и той же датой.

Итак, поместите ваш временной ряд в zoo объекты, merge их, затем используйте моя qplot.zoo функция для их построения.Это будет иметь дело с переключением с zoo на длинный фрейм данных.

Вот пример:

> z1 <- zoo(1:8, 1:8)
> z2 <- zoo(2:8, 2:8)
> z3 <- zoo(4:8, 4:8)
> nm <- list("z1", "z2", "z3")
> z <- zoo()
> for(i in 1:length(nm)) z <- merge(z, get(nm[[i]]))
> names(z) <- unlist(nm)
> z
  z1 z2 z3
1  1 NA NA
2  2  2 NA
3  3  3 NA
4  4  4  4
5  5  5  5
6  6  6  6
7  7  7  7
8  8  8  8
> 
> x.df <- data.frame(dates=index(x), coredata(x))
> x.df <- melt(x.df, id="dates", variable="val")
> ggplot(na.omit(x.df), aes(x=dates, y=value, group=val, colour=val)) + geom_line() + opts(legend.position = "none")

Answer 2

Если вы делаете это только потому, что ggplot2 (как и многие другие), например, фреймы данных, то вам не хватает данных в длинных форматах.Да, вы просто складываете все свои переменные ответа в один столбец, соединенные вместе.Тогда у вас будет 1 или более других столбцов, которые определяют, что отличает эти ответы.Это лучший способ настроить его на такие вещи, как ggplot.

Answer 3

Вы не можете.A data.frame() должен быть прямоугольным;но правила переработки гарантируют, что более короткие векторы будут расширены.

Таким образом, у вас может быть другая ошибка - данные, которые вы хотите rbind, не подходят, может быть?- но трудно сказать, поскольку вы не предоставили воспроизводимый пример.

Редактировать Учитывая ваше обновление, вы получаете точно то, что вы просили: списокимена объединяются на rbind.Если вы хотите, чтобы отображались базовые данные, вам необходимо задействовать get() или другой метод доступа к данным.