如何在面板数据中为不同年份创建单独的直方图？

我有一个数据框报纸_年，其中包含每年一组报纸的发行量（$CIRC）。我想看看这些数字的分布如何随时间变化。所以我想为这些不同的年份创建多个单独的直方图。

我已经尝试过以下方法：

ggplot(newspaper_yearly,aes(x=CIRC))+geom_histogram()+facet_grid(~YEAR==2004)+theme_bw()

但这显示了两个直方图，一个其中 YEAR==2004 为真，另一个为 YEAR=2004 不为真。我只想查看 YEAR=2004 为 true 时的直方图。

编辑： 这是基本数据结构的清理示例：

      YEAR   CIRC
45938 1972  16557
10396 1900   2320
56311 2000   1195
1002  1872   1200
53335 1992  17764
7376  1896   1760
30101 1940 100651
18633 1916  11956
3171  1884   1900
54022 1992   5530
38751 1956   8006
42125 1964  10208
636   1872   1500
48706 1980  18830
22497 1924     NA
28024 1936   7211
7684  1896  21752
56087 2000 107129
43935 1968   9288
34692 1948   5083

我知道我可以制作这样的子集（这实际上是我想要的结果），但我想避免每年制作一个子集。

datahist2000 <- newspaper_yearly[ which(newspaper_yearly$YEAR == "2000"), ]
hist(datahist2000$CIRC)

像这样的事情可能会有所帮助。

par(mfrow=c(3,3))
for(i in levels(d$YEAR)){
datahist <- d[which(d$YEAR == i), ]
hist(datahist$CIRC)}

我使用你的子集方法来解决 for 循环的问题。我不太知道这是否是您想要实现的目标。我想每年“保监会”的条目应该不少吧？否则单独的图没有多大意义，至少对于您提供的数据没有意义。