X 轴标签未在 R 上的基础绘图系统上打印

如何解决X 轴标签未在 R 上的基础绘图系统上打印

“R 版本 4.0.3 (2020-10-10)”

我有一些时间序列数据，我试图在其中查看消费随时间的变化。这是一个reprex。

set.seed(1)
date <-  seq(as.POSIXct('2010-01-01'),as.POSIXct('2010-04-10'),by=86400)
Consumption <- rnorm(99)


Data <- data.frame(date,Consumption)


plot(Data$Consumption~Data$date,type='l') # X-axis labels and ticks not printing
par(yaxt='n')
axis(side = 2,at=seq(-3,3,0.5),labels = seq(-3,0.5)) # This works on the first  plot on the y axis



plot(Data$date~Data$Consumption,type='l') # X-axis refusing to print despite assigning it.
par(xaxt='n')
axis(side = 1,0.5)) # This works on the first  plot

初始 plot() 中输出的图形正是我想要的，除了它没有任何 x 轴标签。

我将 Base Plotting 用于作业而不是日常使用，并且通常会使用 ggplot。 我一直在试图弄清楚为什么 x 轴没有绘制。最初我认为问题出在日期变量上，并尝试用 lubridate::ymd() 清理它。然而，当我为了这个问题而开始制作上述 reprex 时，很明显 X 轴标签和刻度线没有打印。在第二个图中，我将消费变量放在 x 轴上。我很惊讶地发现日期在 Y 轴上整齐地打印出来。

我做错了什么？

解决方法

我很容易看到两个问题：

1. 变化：Consumption

2. 问题在于'par'。当一个块中有多个绘图时，与 ggplot 不同，绘图无法正确处理。删除标准并运行下面它应该可以工作

    set.seed(1)
    date <-  seq(as.POSIXct('2010-01-01'),as.POSIXct('2010-04-10'),by=86400)
    Consumption <- rnorm(100)
    Data <- data.frame(date,Consumption)
    plot(Data$Consumption~Data$date,type='l') 
    plot(Data$date~Data$Consumption,type='l') 

请注意，无论何时定义 par 以及在两个不同的块中运行每个图时，标签都会正确显示。你不会有任何问题。但是当你在一个块中绘制两个图表时，如果你有标准，你总是会遇到问题。

,

如果您想更好地控制轴标签和标题发生的情况，您可以手动创建它们。因此，首先生成一个没有标题和标签的图。然后，使用 axis() 和 mtext() 手动创建它们。在此过程中，您可以使用 par(mar=...) 增加图底部的空间。微调是通过 las、cex.axis 和 line 等参数完成的。最后，您将 mar 重置为其旧值。

您可以使用下面的代码获得更详细的 X 轴标签

### format to day (probably not the best way to do this) 
Data$date2 <-format(Data$date,"%d%b")
Data$date2 <- as.Date(Data$date2,format = "%d%b")

### increase room at bottom of the plot for X-axis title
### the labels will eat up space,if we do nothing it will be a mess
### set title with mtext later
par(mar = c(7,4,2))

### plot without X-axis labels (xaxt) and title (xlab)
### work with "at" and "labels" in axis-function
### rotate labels 90° (las) an reduce font (cex.axis)
### put title 5 lines below graph (line)
###
### Remark: the graph window has to be big enough
plot(Data$Consumption ~ Data$date,type= "l",xaxt = "n",xlab = NA) 
axis(side = 1,at = Data$date,labels =  Data$date2,las = 2,cex.axis=.75)
mtext(side = 1,text = "Date",line = 5)

这会产生以下图表：

替代方案 每 7 个项目的刻度和标签

per7 <- seq(1,99,7)
plot(Data$Consumption ~ Data$date,at = Data$date[per7],labels =  Data$date2[per7],line = 5)

### reset mar
par(mar = c(5,2))

它给出了以下图片：

请告诉我这是否是您想要的。