3

我有 30 次运行数据，每次都存储在一个单独的 CSV 文件中，runi.csv，i = 0:29。

假设我想将它们全部收集到一个列表中。我知道如何做到这一点的最好方法是

runs = list()
for (i in 1:30) { runs[[i]] = read.csv(paste("run", i-1, ".csv")); }

现在让我们进一步说，存储在列表中的每个数据帧都具有相同的列布局，并且我对由“x”标识的列和由“y”标识的列感兴趣。

绘制所有 30 次运行的 (x, y) 对的最简单方法是什么？这是我目前的做法（我觉得必须有更好的方法）：

xList = list()
yList = list()
for (i in 1:30) { xList[[i]] = runs[[i]]$x; yList[[i]] = runs[[i]]$y; }
matplot(x=as.data.frame(xList), y=as.data.frame(yList))

当我尝试对数据进行转换时，这会变得更加痛苦；我不知道如何将函数应用于存储在列表中的每个数据框的特定列。

这里的任何帮助都会非常有帮助。

4

2 回答 2

5

使用所有数据创建一个数据框可能会好得多。例如导入时添加运行号（runs[[i]] = data.frame(read.csv(paste("run", i-1, ".csv")), Run=i)），然后执行alldata <- do.call(rbind, runs)。

现在您可以使用lattice或ggplot2来制作绘图。例如，要通过运行获取使用不同颜色的所有运行的散点图，请执行以下操作：

library(ggplot2)
qplot(x, y, colour=Run, data=alldata, geom="point")
于 2010-02-12T19:53:31.350 回答
3

在处理这样的列表时，最好使用 l*ply 函数（来自 plyr）或 lapply。

最简单的导入方法可能是这样的：

library(plyr)
runs <- llply(paste("run",1:30,".csv",sep=""), read.csv)

这是绘制它们的一种方法：

# some dummy data
runs <- list(a=data.frame(x=1:5, y=rnorm(5)), b=data.frame(x=1:5, y=rnorm(5)))
par(mfrow=c((length(runs)/2),2));
l_ply(1:length(runs), function(i) { plot(runs[[i]]$x, runs[[i]]$y) })

当然，您也可以将其输出到另一台设备（例如 pdf）而不使用par().

于 2010-02-12T19:49:58.233 回答