3

我想自动生成几百个散点图，但我在处理子集的紧凑方法方面没有取得多大成功。数据框具有三个感兴趣的字段：

site           code            value

1              x                30
1              y                14
1              z                12
2              x                35
2              y                19

ETC，

其中站点 = 数字，值 = 数字，代码 = 分组变量。

大约有 400 个站点和 300 个代码。

我想要的是：

  1. 根据字段 =code 对 df 进行子集化
  2. 绘制 x=site vs y=value，作为散点图
  3. 将每个散点图写入 .pdf 中的单独页面
  4. 自动化流程

任何意见表示赞赏。例如，这种方法会生成错误“子集”必须评估为逻辑。lapply 调用也可能不起作用。

dfarray<-list(subset(df,df[,2]))
pdf(file="test.pdf")
figures<-lapply(dfarray, function(i) {(plot(i$[,1],i$[,3]))
})
dev.off()

非常感谢

合并 d_ply 效果很好。我对这个结果很满意，但是谁能指出 X、Y 轴的标签声明有什么问题？（注意：此代码块导致每个子集 1 个 pdf）

names(df)<-tolower(names(df))
names(df)
d_ply(df,"code",function(x){
pdf(file=paste(unique(x$code),".pdf"))
plot(data=x,value~site, xlab=names(x[1]),ylab=as.name(x$code))
dev.off()
})

更新：假设数据框 df 存在，以下是此代码的两个功能版本。这个将每个组的散点图写入单独的.pdf：

names(df)<-tolower(names(df)
names(df)
d_ply(df,"code",function(x){
pdf(file=paste(unique(x$code),".pdf"))
plot(data=x, value~site, xlab=names(df[1]),ylab=unique(x$code))
dev.off()
})

这个将所有散点图（每组一个）写入单个 .pdf：

names(df)<-tolower(names(df)
names(df)
pdf(file="df_all.pdf")
d_ply( df, "code", function(x) {
plot(data=x, value~site, xlab=names(df[1]), ylab=unique(x$code))
})
dev.off()
4

1 回答 1

3

我会去plyr::d_ply：

library(plyr)
pdf(file="test.pdf")
d_ply( df.1, "code", function(x) plot( data = x, value~site, main = unique(x$code) ) )
dev.off()

您可能想大致了解split-apply-(combine) strategy。

于 2013-05-29T23:46:30.410 回答