在R中，如何在一个脚本/代码部分中为多个站点和参数创建多个时间序列图？

Question

我正在尝试为时间序列数据（多年来的年平均水质）创建数百个图（每个站点/程序一个图）。我需要绘制每个站点/程序与每个参数（ANC，DOC等）的组合。年限从1990年到2017年。如果未收集参数，则将其存储为0，但是如果需要，我可以将其更改为null。我想将其作为一种代码/脚本来节省时间。

我对ggplot2更为熟悉，但是我也对其他软件包开放。

这是我的数据示例：

PROGRAM <fctr> SITE <fctr> YEAR <int> ANC <dbl> DOC <dbl> NO3 <dbl> SBC_ALL <dbl> SBC.0 <dbl> SO4 <dbl>

1   VTSSS   VT64    1992    66.753  0   0.803   94.890  0   8.135
2   VTSSS   VT64    1998    64.778  0   0.803   92.783  0   8.230
3   VTSSS   VT43    1995    58.925  0   0.698   91.945  0   8.730
4   VTSSS   VT43    1990    63.240  0   0.885   91.618  0   8.733
5   VTSSS   VT43    1992    66.043  0   0.825   93.873  0   8.758
6   VTSSS   VT43    2004    66.298  0   0.688   95.358  0   8.813

我尝试过ggplot，但是它只带有一个空图。我在想我需要一个绘图软件包，可以在其中按PROGRAM和SITE进行子集。

我尝试过：

ggplot(data = mydata, mapping = aes(x = "year", y = "anc", color = "site")) +
 geom_line() +
 theme_bw()

预期结果：想要绘制一个图：VT43，参数：ANC随时间变化（1990年至2017年）。 VT43和DOC。 VT43和NO3。等等....

实际结果：

空图，x轴为年份，y轴为anc，网站上没有点或线。