使用 ggplot 在 R 中绘制多个大型时间序列时出现问题

Question

我有一个包含 16 个独立时间序列的大型数据集。我想将它们绘制在 3x7 网格中，顶行是以 IN 结尾的每个时间序列，底行是以 OUT 结尾的每个时间序列。在中间行，我将重复两个以 RN 结尾的时间序列中的每一个，它们对应于每个 IN/OUT 对。我已经设置了我认为完成的必要逻辑和语法，但我继续与数据类作斗争，并且无法让脚本运行足够长的时间来调试其余部分。我是 R 的新手，所以我不确定如何继续。

library(ggplot2)

numbers <- read.csv("./AllData.csv", header=TRUE) 
num <- (ncol(numbers) - 4)/4*3 #converts 36 columns to 24 plots involving 8 timeseries

par(mfrow=c(3,7))
for(i in 1:num){
  if (i < 8) {
    qplot(as.POSIXct(as.Date(numbers[4*i+1])), numbers[4*i+2], data="numbers", xlab="Date", ylab="Feet", main=numbers[4*i+2,1])
  } else if (i < 15) {
    qplot(as.POSIXct(as.character(numbers[4*i-6]), format="%m/%d/%Y %H:%M"), numbers[4*i-5], data="numbers", xlab="Date", ylab="Feet", main=numbers[4*i-5,1])
  } else {
    qplot(as.POSIXct(as.character(numbers[4*i-13]), format="%m/%d/%Y %H:%M"), numbers[4*i-12], data="numbers", xlab="Date", ylab="Feet", main=numbers[4*i-12,1])
  }

代替“as.POSIXct”，我还尝试了“as.Date”，“as.POSIXlt”，并简单地按原样绘制数据集，但都会导致类似的错误，通知我 ggplot 无法使用无论我使用什么课程。我尝试使用下面的基本绘图功能，但随后收到错误，我的 x 和 y 长度不同，我已经检查过并且长度相同。

par(mfrow=c(3,7))
for(i in 1:num){
  if (i < 8) {
    x <- as.POSIXct(as.character(numbers[4*i+1]), format="%m/%d/%Y %H:%M")
    x <- numbers[4*i+1]
    y <- numbers[4*i+2]
    main <- numbers[4*i+2,1]
  } else if (i < 15) {
    x <- as.POSIXct(as.character(numbers[4*i-6]), format="%m/%d/%Y %H:%M")
    y <- numbers[4*i-5]
    main <- numbers[4*i-5,1]
  } else {
    x <- as.POSIXct(as.character(numbers[4*i-13]), format="%m/%d/%Y %H:%M")
    y <- numbers[4*i-12]
    main <- numbers[4*i-12,1]
  }
  plot(x, y, type="l", xlab="Date", ylab="Feet", heading=main)
}

任何智慧或指导将不胜感激。这是我的数据集：https://drive.google.com/file/d/0BzEFYP8cpjHpUlluSlNZaTN1RDA/view?usp=sharing

Answer 1

这里有几处我已更改，[（子集）与[[（索引）混淆，ggplot 数据框被指定为字符变量（导致错误消息ggplot 无法处理字符变量），列名（读入names 属性）与第一行等混淆。

说句公道话，这首先发生在每个人身上。我强烈建议您先阅读 Hadley Wickham 的“Advanced-R”编程书籍（或类似的）的前两章，然后再尝试其他如此雄心勃勃的东西。你会明白我的意思，如果一个人不能很好地理解名称（和其他属性）、子集、数据框和列表（数据框是 R 中的专用列表），那么使用 R 是非常困难的。

在线版链接在这里：http://adv-r.had.co.nz/

还可以先尝试为一个小示例工作。这太大了，无法开始。

无论如何这里是更正的代码：

library(ggplot2)

numbers <- read.csv("./AllData.csv", header=TRUE) 
num <- (ncol(numbers) - 4)/4*3 #converts 36 columns to 24 plots involving 8 timeseries

par(mfrow=c(3,7))
for(i in 1:1){
  if (i < 8) {
    qplot(as.Date(numbers[[4*i+1]],"%m/%d/%Y %H:%M"), numbers[[4*i+2]], xlab="Date", ylab="Feet", main=names(numbers)[4*i+2])
#    qplot(as.POSIXct(as.character(numbers[4*i+1]), format="%m/%d/%Y %H:%M"), numbers[[4*i+2]],  xlab="Date", ylab="Feet", main=numbers[4*i+2,1])
  } else if (i < 15) {
    qplot(as.Date(numbers[[4*i-6]],"%m/%d/%Y %H:%M"), numbers[[4*i-5]],  xlab="Date", ylab="Feet", main=names(numbers)[4*i-5])
#    qplot(as.POSIXct(as.character(numbers[4*i-6]), format="%m/%d/%Y %H:%M"), numbers[[4*i-5]],  xlab="Date", ylab="Feet", main=numbers[4*i-5,1])
  } else {
    qplot(as.Date(numbers[[4*i-13]],"%m/%d/%Y %H:%M"), numbers[[4*i-12]], xlab="Date", ylab="Feet", main=names(numbers)[4*i-12])
#    qplot(as.POSIXct(as.character(numbers[4*i-13]), format="%m/%d/%Y %H:%M"), numbers[[4*i-12]], xlab="Date", ylab="Feet", main=numbers[4*i-12,1])
  }
}

这是情节（我将其缩减为一个情节，因为它们需要很长时间）：