У меня есть следующий фрейм данных с курсами времени Y1 для 10 предметов. Каждый вложенный элемент имеет одну числовую, одну нечисловую метку
Sub Time Y1 lab1 lab2
1 1 0.0 0.0 40 weak
2 1 0.5 1.1 40 weak
...
10 1 63 2.0 40 weak
...
91 10 0.0 0.0 23 strong
92 10 0.5 1.2 23 strong
...
100 10 63 1.6 23 strong
Я хотел бы рассчитать медиану Y1 для всех предметов, по группам, то есть для слабых / сильных подгрупп, и построить их вместе с временными рядами Y1 на одной фигуре.
Однако основной проблемой является медиана. Если я распаковываю / выкладываю фрейм данных, используя
DF2 <- unstack(DF1, Y1~Sub)
Затем я могу легко вычислить медиану с помощью rowMedians
, но не могу объединить ее с временными рядами ggplot для отдельных объектов из-за различий в размерах. Чтобы быть более точным, если я сделаю это
p1 <- ggplot(data = DF1,
mapping=aes(x = DF1$Time, y = DF1$Y1, group=DF1$Index)) +
geom_line(aes(color=DF1$lab2)) +
geom_line(data = medianDF, aes(x=medianDF$time,y=medianDF$median))
Я получаю сообщение об ошибке: Error: Aesthetics must be either length 1 or the same as the data (10): group
.
Любые идеи будут оценены.