Ось X для раз

Question

Я пытаюсь создать серию графиков, на которых один и тот же пациент принимает напитки и мочится в разное время.Каждый сюжет представляет один день.Я хочу сравнить дни и, следовательно, мне нужно убедиться, что все построенные графики имеют одинаковую ось X.Мой код ниже которого я набрал из Как указать фактические значения оси X для построения в виде отметок оси X в R

### Data Input

time_Thurs <- c("01:10", "05:50", "06:00","06:15", "06:25", "09:35", "10:00", "12:40",
              "14:00", "17:20", "18:50", "19:10", "20:10", "21:00", "22:05", "22:35")
event_Thurs <- c("u", "u", "T", "T", "u", "u", "T","T","u", "u", "T", "T", "T", "T", "u", "W")
volume_Thurs <- c(NA, NA, 0.25, 0.25, NA, NA, 0.125, 0.625, NA, NA, 0.25, 0.25, 0.25, 0.25,
            NA, 0.25)
total_liquids_Thurs <- sum(volume_Thurs, na.rm=TRUE)
time_Thurs <- paste("04/04/2019", time_Thurs, sep=" ")

time_Fri <- c("01:15", "06:00", "06:10", "06:25", "06:30", "07:10", "08:40", "09:20",
              "12:45", "13:45")
event_Fri <- c("u","u", "T","T","u","uu","T", "u", "T", "u")
volume_Fri <- c(NA, NA, 0.25, 0.25, NA, NA, 0.125, NA, 0.625, NA)
total_liquids_Fri <- sum(volume_Fri, na.rm=TRUE)
time_Fri <- paste("05/04/2019", time_Fri, sep=" ")

### Collect all data together

event <- c(event_Thurs, event_Fri)
Volume <- c(volume_Thurs, volume_Fri)
time_log <- c(time_Thurs, time_Fri)
time_log <- strptime(time_log, format = "%d/%m/%Y %H:%M")
time_view <- format(time_log, "%H:%M")

### Put into Dataframe

patient_data <- data.frame(time_log, time_view, event, Volume)

# write.csv(patient_data, file="patient_data.csv", row.names = FALSE)

daily_plot <- function(x, day) {
  #    x patient data - a data.frame with four columns:
  #    POSIXct time, time, event and Volume
  #    date number of day of month
  #    y volume of liquid
  #    TotVol total volume of intake over week
  #    Event - drink or otherwise

  x <- x[as.numeric(format(x[,1], "%d")) == day, ]
  TotVol <- sum(x[,4], na.rm = TRUE)
  DayOfWeek <- weekdays(x[1,1], abbreviate = FALSE)

  plot(x[,1],x[,4], 
       xlim = c(x[1,1],x[length(x[,1]),1]),
       xlab="Hours of Study", ylab = "Volume of Liquid Drank /L",
       main = paste("Total Liquids Drank = ", TotVol, " L on ", DayOfWeek, "Week 1, Apr 2019"),
       sub = "dashed red line = urination", pch=16,
       col = c("black", "yellow", "green", "blue")[as.numeric(x[,3])],
       xaxt = 'n'
  )
  xAxis_hrs <- seq(as.POSIXct(x[1,1]), as.POSIXct(x[length(x[,1]),1]), by="hour")
  axis(1, at = xAxis_hrs, las = 2)
  abline( v = c(x[x[,3] == "u",1]), lty=3, col="red")
}

Когда я запускаю функцию,

daily_plot(patient_data, 4)

Я хочу распечатать свою ось X с поправками в виде часов, представляющих события за 24-часовой период.

Когда я обертываю свой вектор xAxis_hrs в strptime (xAxis_hrs, format = "% H ") код вылетает - то есть ось x не распечатывается, и я вижу, Ошибка в оси (1, at = xAxis_hrs, las = 2): (список) объект не может быть приведен к типу 'double',Любая помощь?

Answer 1

Проблема заключается в том, что вы передаете метки неверно названному аргументу, а именно at (это должны быть числовые позиции меток). Вместо этого используйте следующее:

axis(1, at = xAxis_hrs, labels = strptime(xAxis_hrs, format = "%H"), las = 2)

К сожалению, это не меняет того факта, что метки осей не вписываются в график и конфликтуют с заголовком оси. Первый может быть исправлен путем настройки полей графика. Я не знаю хорошего решения для последнего, хотя изменение формата времени могло бы помочь: вероятно, нет необходимости / полезно печатать полные минуты и секунды (которые всегда равны 0). На самом деле, вы хотели использовать format вместо strptime?

Кроме того, я принципиально согласен с другим ответом, рекомендующим ggplot2 в долгосрочной перспективе. Это делает такие вещи менее болезненными.

Answer 2

Если вы открыты для решения ggplot:

library(tidyverse)
library(lubridate)

daily_ggplot <- function(df, selected_day) {
  df_day <- filter(df, day(time_log) == selected_day)
  df_urine <- filter(df_day, event == "u")
  df_drink <- filter(df_day, event != "u")
  TotVol <- sum(df_day$Volume, na.rm = TRUE)
  Date <- floor_date(df_day$time_log[1], 'days')
  DayOfWeek <- weekdays(Date, abbreviate = F)

  plot_title <- paste0("Total drank = ", TotVol, "L on ", DayOfWeek, " Week 1, Apr 2018")

  ggplot(df_drink) +
    aes(time_log, Volume, color = event) +
    geom_point(size = 2) +
    geom_vline(data = df_urine, aes(xintercept = time_log), color = "red", linetype = 3) +
    labs(x = "Hours of Study", ylab = "Volume of Liquid Drank (L)",
         title = plot_title, subtitle = "lines = urination") +
    theme_bw() +
    scale_x_datetime(date_labels = "%H:%M", limits = c(Date, Date + days(1)))
}

daily_ggplot(patient_data, 4)