Построение нескольких графиков временных рядов в R

Question

Может кто-нибудь помочь.У меня большой набор данных, и мне нужно построить несколько графиков временных рядов.

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

ts.plot(bc,gpars= list(col=rainbow(10)))

данные ниже, если кто-нибудь может помочь?

TIME      15         16           17           18          20          22          23
08:00:00 130.4905899    15.44164769 948.9185939 6211.837354 1071.730556 10.08920076 7.793301031
11:00:00 125.7301547    18.87143833 991.0009783 6304.471569 1082.126629 10.80475415 7.857773565
14:00:00 153.3779662    17.63335938 949.1741247 6209.524186 1079.102756 10.68438383 8.326058855
17:00:00 132.3256891    15.8961511  917.0452991 6123.402395 1081.166439 10.41007094 7.856372445
20:00:00 130.6405835    15.28122651 917.0229181 6135.679239 1084.589394 10.70688202 7.741277402
08:00:00 124.1484465    17.14357927 948.9060481 6126.791479 1085.907147 10.76713085 7.810187162
08:00:00 161.0455657    17.10409992 881.4517913 5839.73355  1073.585164 9.925269955 7.987206082
08:00:00 165.645823     16.45928764 860.4285647 5781.612679 1073.439013 10.01297791 7.983672272

Я искал и не могу найти то, что мне нужно сделать?

Заранее спасибо.

Код яиспользуется следующим образом:

library(tidyverse)
df <- read.csv("Strep_time_series_NBS.csv", header = TRUE, row.names = NULL)
timestamps <- lubridate::dmy_h(c(3008201808, 3008201811,     3008201814,    3008201817, 3008201820,3108201808, 0309201808, 1009201808))

df_fix_times <- df %>%
  mutate(TIME=timestamps)

df_tidy <- df_fix_times %>%
  gather(series, value, -TIME)

ggplot(df_tidy, aes(TIME, value, group =series))+
  geom_line() +
  facet_wrap(~series, scales = "free_y")

library(sqldf)

df2=sqldf("What needs to go in here?")   


library(officer);library(ggplot2);library(magrittr);library(rvg)


doc <- read_pptx()
for (i in 1:nrow(df))
{
  data_graph_ind=df_tidy[(df_tidy$series==df[i,1]),]
  doc <- add_slide(doc, layout = "Two Content", master = "Office Theme")
  doc <- doc %>% 
  ph_with_text(type = "title",df2[i,1]) 

  doc <- doc %>%
  ph_with_vg(doc, ggobj=ggplot(data = data_graph_ind, aes(x = TIME, y    =      value)) + 
             geom_line()+geom_point(),type="body")
  print(doc,target="C:/path/example.pptx")

  }

Спасибо

Answer 1

Другой подход, если вы хотите создавать графики на отдельных слайдах:

Шаг 1) Считать данные

    df1 <-
  read.table(header = TRUE, row.names = NULL, text = 
               "TIME    15  16  17  18  20  22  23
0 hrs   130.4905899 15.44164769 948.9185939 6211.837354 1071.730556 10.08920076 7.793301031
3 Hrs   125.7301547 18.87143833 991.0009783 6304.471569 1082.126629 10.80475415 7.857773565
6 Hrs   153.3779662 17.63335938 949.1741247 6209.524186 1079.102756 10.68438383 8.326058855
8 Hrs   132.3256891 15.8961511  917.0452991 6123.402395 1081.166439 10.41007094 7.856372445
12 Hrs  130.6405835 15.28122651 917.0229181 6135.679239 1084.589394 10.70688202 7.741277402
24 Hrs  124.1484465 17.14357927 948.9060481 6126.791479 1085.907147 10.76713085 7.810187162
96 Hrs  161.0455657 17.10409992 881.4517913 5839.73355  1073.585164 9.925269955 7.987206082
168 hrs 165.645823  16.45928764 860.4285647 5781.612679 1073.439013 10.01297791 7.983672272")

Шаг 2)

library(ggplot2);library(tidyr);library(dplyr)
df_tidy <- df1 %>%
  mutate(TIME = row.names %>% as.integer) %>%
  select(-row.names) %>%
  gather(series, value, -TIME)

Шаг 3) Создание графиков на слайдах один за другим в ppt-документе

library(sqldf)

df2=sqldf("select distinct(series) from df_tidy")   


library(officer);library(ggplot2);library(magrittr);library(rvg)


doc <- read_pptx()
for (i in 1:nrow(df2))
                  {
                     data_graph_ind=df_tidy[(df_tidy$series==df2[i,1]),]
                     doc <- add_slide(doc, layout = "Two Content", master = "Office Theme")
                     doc <- doc %>% 
                     ph_with_text(type = "title",df2[i,1]) 

                     doc <- doc %>%
                     ph_with_vg(doc, ggobj=ggplot(data = data_graph_ind, aes(x = TIME, y = value)) + 
                                    geom_line()+geom_point(),type="body")
                     print(doc,target="C:/path/example.pptx")

                  }

Answer 2

(РЕДАКТИРОВАТЬ: обновлено, чтобы использовать пересмотренные данные образца с часами.)

Вот один подход.Сначала загрузите прибор и прочитайте данные:

library(tidyverse)   

df <- read.table(
    header = TRUE, row.names = NULL, text =
      "TIME      15         16           17           18          20          22          23
08:00:00 130.4905899    15.44164769 948.9185939 6211.837354 1071.730556 10.08920076 7.793301031
             11:00:00 125.7301547    18.87143833 991.0009783 6304.471569 1082.126629 10.80475415 7.857773565
             14:00:00 153.3779662    17.63335938 949.1741247 6209.524186 1079.102756 10.68438383 8.326058855
             17:00:00 132.3256891    15.8961511  917.0452991 6123.402395 1081.166439 10.41007094 7.856372445
             20:00:00 130.6405835    15.28122651 917.0229181 6135.679239 1084.589394 10.70688202 7.741277402
             08:00:00 124.1484465    17.14357927 948.9060481 6126.791479 1085.907147 10.76713085 7.810187162
             08:00:00 161.0455657    17.10409992 881.4517913 5839.73355  1073.585164 9.925269955 7.987206082
             08:00:00 165.645823     16.45928764 860.4285647 5781.612679 1073.439013 10.01297791 7.983672272"
  )

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

timestamps <- lubridate::ymd_h(c(
  2018091308, 2018091311, 2018091314, 2018091317, 
  2018091320, 2018091408, 2018091708, 2018092008))

df_fix_times <- df %>%
  mutate(TIME = timestamps)

Теперь мы можем преобразовать в «приборку»формат, где каждое наблюдение находится в ряду.Это облегчит построение графика.

df_tidy <- df_fix_times %>%
  gather(series, value, -TIME)

И теперь мы можем использовать вызов facet_wrap здесь, чтобы поместить каждую серию в свой собственный график.

ggplot(df_tidy, aes(TIME, value, group = series)) +
  geom_line() +
  facet_wrap(~series, scales = "free_y")