R语言 如何在ggplot2中绘制自定义的每小时气象数据

svmlkihl  于 2023-07-31  发布在  其他
关注(0)|答案(1)|浏览(73)

我做了一个星期的实验。为了确定每个星期，我使用了一周的第一个日期。例如，下图显示20周实验数据。

我的天气数据是这样的

我希望每小时的天气数据有完全相同的图表（而不是每周一列，我有更多的列，每列代表一小时的天气数据）。为了使x轴类似于上面的实验图，我使用mutate(date2 = first(date))在date列的基础上制作了一个date2列，以确保在x轴上使用一周的第一个日期/天。问题是现在我不能绘制每小时的天气数据，因为scale_x_datetime()只适用于日期。我的date2对象是一个因子。
有没有办法解决这个问题？下面是可重复的示例

df <-
  structure(
    list(
      date = c(
        "05/26/2014 11:00 AM",
        "05/26/2014 12:00 PM",
        "05/26/2014 01:00 PM",
        "06/09/2014 12:00 PM",
        "06/09/2014 01:00 PM",
        "06/09/2014 02:00 PM",
        "06/16/2014 12:00 PM",
        "06/16/2014 01:00 PM",
        "06/16/2014 02:00 PM",
        "06/23/2014 12:00 PM",
        "06/23/2014 01:00 PM",
        "06/23/2014 02:00 PM",
        "06/30/2014 12:00 PM",
        "06/30/2014 01:00 PM",
        "06/30/2014 02:00 PM",
        "07/07/2014 12:00 PM",
        "07/07/2014 01:00 PM",
        "07/07/2014 02:00 PM",
        "07/14/2014 12:00 PM",
        "07/14/2014 01:00 PM",
        "07/14/2014 02:00 PM",
        "07/21/2014 12:00 PM",
        "07/21/2014 01:00 PM",
        "07/21/2014 02:00 PM",
        "07/28/2014 12:00 PM",
        "07/28/2014 01:00 PM",
        "07/28/2014 02:00 PM",
        "08/04/2014 12:00 PM",
        "08/04/2014 01:00 PM",
        "08/04/2014 02:00 PM",
        "08/11/2014 12:00 PM",
        "08/11/2014 01:00 PM",
        "08/11/2014 02:00 PM",
        "09/08/2014 12:00 PM",
        "09/08/2014 01:00 PM",
        "09/08/2014 02:00 PM",
        "09/15/2014 12:00 PM",
        "09/15/2014 01:00 PM",
        "09/15/2014 02:00 PM",
        "09/22/2014 12:00 PM",
        "09/22/2014 01:00 PM",
        "09/22/2014 02:00 PM",
        "09/29/2014 12:00 PM",
        "09/29/2014 01:00 PM",
        "09/29/2014 02:00 PM",
        "10/06/2014 12:00 PM",
        "10/06/2014 01:00 PM",
        "10/06/2014 02:00 PM",
        "10/13/2014 12:00 PM",
        "10/13/2014 01:00 PM",
        "10/13/2014 02:00 PM",
        "10/20/2014 12:00 PM",
        "10/20/2014 01:00 PM",
        "10/20/2014 02:00 PM",
        "10/27/2014 12:00 PM",
        "10/27/2014 01:00 PM",
        "10/27/2014 02:00 PM",
        "11/03/2014 12:00 PM",
        "11/03/2014 01:00 PM",
        "11/03/2014 02:00 PM"
      ),
      date2 = structure(
        c(
          1401104700,
          1401104700,
          1401104700,
          1402317900,
          1402317900,
          1402317900,
          1402920000,
          1402920000,
          1402920000,
          1403524800,
          1403524800,
          1403524800,
          1404129600,
          1404129600,
          1404129600,
          1404734400,
          1404734400,
          1404734400,
          1405339200,
          1405339200,
          1405339200,
          1405944000,
          1405944000,
          1405944000,
          1406548800,
          1406548800,
          1406548800,
          1407153600,
          1407153600,
          1407153600,
          1407758400,
          1407758400,
          1407758400,
          1410177600,
          1410177600,
          1410177600,
          1410782400,
          1410782400,
          1410782400,
          1411387200,
          1411387200,
          1411387200,
          1411992000,
          1411992000,
          1411992000,
          1412596800,
          1412596800,
          1412596800,
          1413201600,
          1413201600,
          1413201600,
          1413806400,
          1413806400,
          1413806400,
          1414411200,
          1414411200,
          1414411200,
          1415016000,
          1415016000,
          1415016000
        ),
        tzone = "UTC",
        class = c("POSIXct",
                  "POSIXt")
      ),
      week = structure(
        c(
          1L,
          1L,
          1L,
          2L,
          2L,
          2L,
          3L,
          3L,
          3L,
          4L,
          4L,
          4L,
          5L,
          5L,
          5L,
          6L,
          6L,
          6L,
          7L,
          7L,
          7L,
          8L,
          8L,
          8L,
          9L,
          9L,
          9L,
          10L,
          10L,
          10L,
          11L,
          11L,
          11L,
          12L,
          12L,
          12L,
          13L,
          13L,
          13L,
          14L,
          14L,
          14L,
          15L,
          15L,
          15L,
          16L,
          16L,
          16L,
          17L,
          17L,
          17L,
          18L,
          18L,
          18L,
          19L,
          19L,
          19L,
          20L,
          20L,
          20L
        ),
        levels = c(
          "1",
          "2",
          "3",
          "4",
          "5",
          "6",
          "7",
          "8",
          "9",
          "10",
          "11",
          "12",
          "13",
          "14",
          "15",
          "16",
          "17",
          "18",
          "19",
          "20",
          "21",
          "22",
          "23",
          "24",
          "25",
          "26"
        ),
        class = "factor"
      ),
      mean_temperature = c(
        25,
        25.25,
        26.075,
        26.2,
        26.35,
        27.65,
        29.95,
        29.4,
        30.5,
        22.75,
        24,
        25.275,
        25.6,
        26.8,
        27.125,
        27.5,
        27.525,
        26.1,
        30.65,
        30.65,
        30.425,
        26.15,
        24.5,
        23.7,
        25.65,
        25.85,
        25.2,
        26.975,
        27.125,
        28.275,
        18.55,
        18.775,
        19.125,
        18.8,
        20,
        19.375,
        17.825,
        17.925,
        18.25,
        18.975,
        20.825,
        20.1,
        18.1,
        18.4,
        19.2,
        16.325,
        18.4,
        18.75,
        14.675,
        14.925,
        15.3,
        15.9,
        17.175,
        17.525,
        22.45,
        23,
        23.475,
        17.8,
        18.05,
        17.725
      )
    ),
    class = c("grouped_df", "tbl_df", "tbl",
              "data.frame"),
    row.names = c(NA,-60L),
    groups = structure(
      list(
        week = structure(
          1:20,
          levels = c(
            "1",
            "2",
            "3",
            "4",
            "5",
            "6",
            "7",
            "8",
            "9",
            "10",
            "11",
            "12",
            "13",
            "14",
            "15",
            "16",
            "17",
            "18",
            "19",
            "20",
            "21",
            "22",
            "23",
            "24",
            "25",
            "26"
          ),
          class = "factor"
        ),
        .rows = structure(
          list(
            1:3,
            4:6,
            7:9,
            10:12,
            13:15,
            16:18,
            19:21,
            22:24,
            25:27,
            28:30,
            31:33,
            34:36,
            37:39,
            40:42,
            43:45,
            46:48,
            49:51,
            52:54,
            55:57,
            58:60
          ),
          ptype = integer(0),
          class = c("vctrs_list_of",
                    "vctrs_vctr", "list")
        )
      ),
      row.names = c(NA,-20L),
      .drop = TRUE,
      class = c("tbl_df",
                "tbl", "data.frame")
    )
  )

基本上，我试图修改下图以生成顶部给出的图（我的实验数据图），因此我将使x-axis与我的实验数据图中所示的相同。我会有每周每小时的天气数据列，而不是每天。

r

1条答案

按热度按时间
x0fgdtte

x0fgdtte1#

像这样吗？
从date中提取hour，将前者强制为因子，将后者强制为日期对象。然后按hour进行地块分组。y轴按周中断。
请注意，如果需要，我一次处理date一个步骤，以使其更可读，更易于调试。

suppressPackageStartupMessages({
  library(dplyr)
  library(lubridate)
  library(ggplot2)
})

df %>% 
  ungroup() %>%
  mutate(
    date = mdy_hm(date),
    hour = hour(date),
    date = ymd(format(date, "%Y-%m-%d")),
    hour = factor(hour)
  ) %>% 
  ggplot(aes(date, mean_temperature, group = hour)) +
  geom_col(position = position_dodge(), fill = "steelblue", color = "white") +
  scale_x_date(date_breaks = "1 week") +
  coord_flip() +
  theme_bw()

字符串


的数据
创建日期：2023年7月30日，使用reprex v2.0.2

赞(0）分享  回复(0）举报  2023-07-31
我来回答

相关问题

    查看更多

    热门标签

    更多
    JavaquerypythonNode开发语言requestUtil数据库Table后端算法LoggerMessageElementParser

    最新问答

    更多
    • 蜀ICP备13028337号-1 大数据知识库 https://www.saoniuhuo.com © All rights reserved
    • 本站内容来源互联网,如果侵犯您的权益请联系我们删除, 联系方式：448109455@qq.com