我有几个关于绘制数据的一般性问题。首先，我使用rbind整理我的所有数据 - 包括时间，动物的长度，地点，年份和湖泊。
time(days) L Site Year Loch
1 2.3 LM 2017 Leven
2 2.34 LM 2017 Leven
...
729 5.09 LM 2017 Leven
730 5.1 LM 2017 Leven
1 2.33 LM 2020 Leven
2 2.343 LM 2020 Leven
...
729 5.228 LM 2017 Leven
730 5.229 LM 2020 Leven
1 2.33 LM 2030 Leven
我使用模拟的气候变化温度来强制我的模型每十年一直到2060年。正如您所看到的，每个站点都有模拟数据，每个十年有730天。因此，我为每个站点提供了5套730天的数据集（2017年，2020年，2030年，2040年，2050年和2060年）。同样地，我有来自2个湖泊（leven和etive）和6个站点（每个湖泊3个）的数据，总共5840个观测值。
如何绘制图表以便按每个站点绘制模型及其对应的年份标签或图例？
现在我有一些看起来像这样的东西：
qplot(Time, Length, data=Future_Model_Data, colour=Year)
您建议使用哪种测试来显示时间序列数据之间的变化或差异？我正在研究格兰杰测试。