ggplot绘图时x轴用日期字符串出现UseMethod("depth")错误求助

我之前也碰到过一模一样的错误，咱们先理清楚问题出在哪，再一步步解决：

你的问题场景 #

你要绘制以localminute为X轴、meter_value为Y轴的散点图，数据集有11.9K行，数据片段如下：

localminute,dataid,meter_value
2016-01-01 05:18:49.497023-06,35,98178
2016-01-01 05:45:43.51577-06,35,98182
2016-01-01 05:57:45.024395-06,35,98184
2016-01-01 06:08:41.532267-06,35,98186
2016-01-01 06:36:50.053034-06,35,98188
2016-01-01 06:56:53.066606-06,35,98192
2016-01-01 08:57:53.65758-06,35,98212

使用这段ggplot代码时触发了错误：

test <- read.csv("test_35.csv") #read data in csv file, to plot data
ggplot(test, aes(x=localminute, y=meter_value)) + 
  geom_point() + 
  ggtitle("meter value for dataID=35") + 
  theme(axis.text.x = element_text(angle=90, hjust=1))

错误信息：

Error in UseMethod("depth") : no applicable method for 'depth' applied to an object of class "NULL"

问题根源&解决方案 #

这个错误最常见的两个原因：要么数据根本没读取成功（test变成了NULL），要么时间列的格式不对导致ggplot无法处理。咱们逐一排查：

1. 先确认数据是否正确加载 #

首先运行这两行代码，看看test里有没有数据：

# 打印前6行数据，确认读取成功
print(head(test))
# 查看各列的数据类型
str(test)

• 如果输出显示test是NULL，那就是文件路径错误！你得确保test_35.csv在当前R的工作目录里，或者直接写完整路径，比如read.csv("C:/Users/xxx/Documents/test_35.csv")（Windows）或者read.csv("~/Documents/test_35.csv")（Mac/Linux）。
• 如果能看到数据，但localminute的类型是chr（字符型），那继续下一步。

2. 把字符型时间转换成时间格式 #

read.csv会自动把时间字符串当成普通字符，ggplot处理字符型X轴时很容易出问题，而且时间也不会按顺序排列。咱们用lubridate包来解析带时区的时间：

# 先安装lubridate（如果还没装过）
install.packages("lubridate")
library(lubridate)

# 转换localminute列为POSIXct时间类型
test$localminute <- ymd_hms(test$localminute)

# 再检查类型，应该显示POSIXct/POSIXt
str(test$localminute)

3. 重新绘制图表 #

现在时间列是正确的时间类型了，再运行你的ggplot代码就没问题了。另外给你加个小优化：数据量很大，把点调小一点避免重叠：

library(ggplot2)
ggplot(test, aes(x=localminute, y=meter_value)) + 
  geom_point(size=1) + # 缩小点的尺寸
  ggtitle("meter value for dataID=35") + 
  theme(axis.text.x = element_text(angle=90, hjust=1)) +
  xlab("Local Time") + ylab("Meter Value") # 给坐标轴加清晰的标签

额外优化：处理大数据量的视觉问题 #

11.9K个点堆在一起可能会看不清趋势，你可以试试这两个方案：

• 抖动散点图：给点加一点随机偏移，减少重叠

ggplot(test, aes(x=localminute, y=meter_value)) + 
  geom_jitter(size=0.5, alpha=0.5) + # alpha设置透明度
  ggtitle("meter value for dataID=35 (Jittered)") + 
  theme(axis.text.x = element_text(angle=90, hjust=1))

• 时间聚合：按小时/分钟聚合计算均值，看整体趋势

# 按小时聚合，计算每个小时的meter_value均值
test$hour <- floor_date(test$localminute, "hour")
ggplot(test, aes(x=hour, y=meter_value)) + 
  stat_summary(fun="mean", geom="line") +
  ggtitle("Average Meter Value per Hour (dataID=35)") + 
  theme(axis.text.x = element_text(angle=90, hjust=1))

内容的提问来源于stack exchange，提问作者SKmay