绘制具有不同y轴的两个时间序列：一个作为点图（或条形图），另一个作为线

Question

我有两个时间序列的数据，每个数据都有不同的值范围。我想将一个绘制为一个dotplot，另一个绘制为dotplot上的一条线。 （我会选择一个看起来很体面的条形图和条形图上的一条线，但我的偏好是一个点图。）

#make some data
require(lubridate)
require(ggplot)

x1 <- sample(1990:2010, 10, replace=F)
x1 <- paste(x1, "-01-01", sep="")
x1 <- as.Date(x1)
y1 <- sample(1:10, 10, replace=T)
data1 <- cbind.data.frame(x1, y1)

year <- sample(1990:2010, 10, replace=F)
month <- sample(1:9, 10, replace=T)
day <- sample(1:28, 10, replace=T)

x2 <- paste(year, month, day, sep="-")
x2 <- as.Date(x2)
y2 <- sample(100:200, 10, replace=T)
data2 <- cbind.data.frame(x2, y2)
data2 <- data2[with(data2, order(x2)), ]

# frequency data for dot plot
x3 <- sample(1990:2010, 25, replace=T)
data2 <- as.data.frame(x3)

我可以用ggplot中的一个数据集制作一个点图或条形图：

ggplot() + geom_dotplot(data=data2, aes(x=x3))

ggplot() + geom_bar(data=data, aes(x=x1, y=y1), stat="identity")

但我无法覆盖第二个数据集，因为ggplot不允许第二个y轴。

我无法弄清楚如何使用barplot（）绘制时间序列。

我可以使用plot（）将第一组数据绘制为“h”类型的绘图，并将第二组数据添加为一条线，但我不能使条形任何更粗，因为每一条对应于一年多的一天，我认为这很难看。

plot(data$x1, data$y1, type="h")
par(new = T)
plot(data2$x2, data2$y2, type="l", axes=F, xlab=NA, ylab=NA)
axis(side=4)

有什么想法吗？我唯一的想法是制作两个单独的图并将它们叠加在图形程序中。 ：/

Answer 1

一个简单的解决方法是遵循base绘图本能并为lwd加强type='h'。请务必设置lend=1以防止圆线：

par(mar=c(5, 4, 2, 5) + 0.1)
plot(data1, type='h', lwd=20, lend=1, las=1, xlab='Date', col='gray',
     xlim=range(data1$x1, data2$x2))
par(new=TRUE)
plot(data2, axes=FALSE, type='o', pch=20, xlab='', ylab='', lwd=2,
     xlim=range(data1$x1, data2$x2))
axis(4, las=1)
mtext('y2', 4, 3.5)

Answer 2

我删除了原来的答案。

---

要回答有关制作点图的问题，您可以重新排列数据，以便使用基本绘图功能。一个例子：

使用chron包进行绘图：

library(chron)

虚拟数据：

count.data <- data.frame("dates" = c("1/27/2000", "3/27/2000", "6/27/2000", "10/27/2000"), "counts" = c(3, 10, 5, 1), stringsAsFactors = F) 

复制列表中的日期：

rep.dates <- sapply(1:nrow(count.data), function(x) rep(count.data$dates[x], count.data$counts[x]))

将计数转换为序列：

seq.counts <- sapply(1:nrow(count.data), function(x) seq(1, count.data$counts[x], 1))

策划：

plot(as.chron(rep.dates[[1]]), seq.counts[[1]], xlim = c(as.chron("1/1/2000"), as.chron("12/31/2000")),
 ylim = c(0, 20), pch = 20, cex = 2)
for(i in 2:length(rep.dates)){
   points(as.chron(rep.dates[[i]]), seq.counts[[i]], pch = 20, cex = 2)
}