Question

我有一个大型数据文件（58列，27,000行），我尝试按一个变量（年份--30个级别）进行子集，然后为30个类别中的每个类别直方图另一个变量（日历）。我尝试了几种选择。

首先，我尝试了“应用”＃39;函数，定义＆＃34; hist＆＃34;作为我的乐趣;

    >tapply(calenday,Year,hist(calenday,main=unique(Year)))

但这只留下了一张图。

然后我试着循环它。我尝试了几种方法对我的数据帧进行子集化，然后对直方图函数应用for循环，但它没有工作。我最好只得到一张图。

然后我尝试了格子condiotonal＆＃34;直方图＆＃34;命令：

     >histogram(~calenday|factor(Year))

但堆叠在一起的30个图表有点太多了。

所以我尝试循环它：

     for (i in 1984:2013) {   
         histogram(~calenday | Year, subset=(Year==i), br=28,    
         main=paste("plot of", i), type="count")}

但是没有任何结果，尽管没有给出错误信息。一旦你把它作为一个脚本函数编写，我就读到某个地方存在问题，但我无法弄清楚如何解决它（没有解决方案）。

最后，我只是尝试了：

    tapply(calenday,Year,hist)

这确实有效，但制作的图表非常基本，标题为数字，间隔太大等等。

有什么建议吗？谢谢！

Answer 1

在不知道你的数据和情节的情况下，从我这边做了一些猜测。

但是，我会使用ggplot及其出色的facet_grid()和facet_wrap()功能。

请查看此文档，了解图表是否与您的需求相似。

祝你好运