可复制的表如下~~

贴上部分原始数据，谢谢~~

ID        Station        Year        Month        COD        AN
1333        20        2005        11        0.4        0.05
487        25        2007        10        0.5        0.14
525        21        2005        5        0.5        0.05
980        9        2003        10        0.5        0.05
1015        30        2007        4        0.5        0.05
1019        30        2004        5        0.5        0.05
1027        30        2007        6        0.5        0.05
1057        30        2006        11        0.5        0.05
1360        2        2010        6        0.5        0.06
1685        24        2007        10        0.5        0.18
1211        17        2005        3        0.6        0.05
1218        17        2008        4        0.6        0.29
1207        17        2009        2        0.7        0.05
1215        17        2009        3        0.7        0.05
1219        17        2009        4        0.7        0.05
614        4        2008        8        0.74        0.09
612        4        2008        7        0.8        0.07
973        9        2005        9        0.8        0.05
1202        17        2009        1        0.8        0.05
1239        17        2009        7        0.8        0.05
1245        17        2004        9        0.8        0.05
120        19        2010        11        0.9        0.308
987        9        2006        11        0.9        0.05
1201        17        2008        1        0.9        0.05
1205        17        2007        2        0.9        0.05
1209        17        2003        3        0.9        0.05
1212        17        2006        3        0.9        0.05
1220        17        2010        4        0.9        0.05
942        9        2009        3        1        0.05
946        9        2008        4        1        0.05
991        9        2010        11        1        0.107
1198        17        2005        1        1        0.05
1214        17        2008        3        1        0.07
1227        17        2009        5        1        0.05
1230        17        2008        6        1        0.05
1260        17        2007        11        1        0.05
1265        17        2008        12        1        0.05
111        19        2010        9        1.1        0.05
972        9        2004        9        1.1        0.1
994        9        2009        12        1.1        0.05
1210        17        2004        3        1.1        0.05
1222        17        2004        5        1.1        0.11
1234        17        2004        7        1.1        0.17
2067        18        2003        3        1.1        0.05
2080        18        2005        11        1.1        0.22

假设数据集名字是dt,那么用如下命令吧
ldt<-split(dt,dt$Station)
拆分后的数据，可以如下这么般用
ldt[[1]]
ldt[[2]]
……

这其实是一个在同一数据集分组（站点）似合模型的工作？！
若是，用dplyr包中的by_group与do 二函数，即可，相当简单。。具体可参考dplyr包的do 的说明例子