用 R 中的因子绘制折线图

Question

我正在用 R 处理这个图表：

但是，正如您所见，我得到的是这些实线破折号，而不是漂亮的线条。这是我用来制作此图表的代码：

par(mar = c(5,5,2,5))
    with(bedtimes, plot(trap, funestus, type="l", col="red3", 
                 ylab=expression(italic(p))),
                 ylim=c(0,3))  

    par(new = T)
    with(bedtimes, plot(trap, bed, pch=16, axes=F, xlab=NA, ylab=NA, cex=1.2))
    axis(side = 4)
    mtext(side = 4, line = 3, 'Proportion in bed')
    polygon(bedtimes$bed,col=transp("gray", 0.3), border = NA)

这是我正在使用的数据的输出：

(removed)

我意识到这是因为我的 x 轴是一个因素而不是数字。然而，试图改变这一点（例如，使用as.POSIXct(paste0("2016-07-12",bedtimes$trap)) 给我带来了各种各样的问题，因为我必须确保 R 最初使用 bedtimes$trap <- factor(bedtimes$trap, levels = bedtimes$trap) 以正确的顺序绘制了这些因素

如何生成同样的图表，但使用线条而不是这些破折号？ 我最终想要一个看起来与此相似的图表，给你一个想法（虽然不完全相同）：

谢谢！

Answer 1

正如我在 cmets 中所说，您可以尝试使用因子变量的级别：

par(mar = c(5,5,2,5))
with(bedtimes, plot(as.numeric(trap), funestus, type="l", col="red3", 
                    ylab=expression(italic(p))),
                    ylim=c(0,3))  
par(new = T)
with(bedtimes, plot(as.numeric(trap), bed, type="l", axes=F, xlab="", ylab=NA, cex=1.2))
axis(4)
axis(side = 1, at=1:length(levels(bedtimes$trap)), levels(bedtimes$trap))
mtext(side = 4, line = 3, 'Proportion in bed')
polygon(bedtimes$bed,col="gray", border = NA)