各位大侠, 现有一组土壤水分数据,包括5个肥料水平(0-180kg ha-1),每个水平3次重复,测定5年(y1-y5)。土壤水分变异一方面来自于肥料处理，另一方面也可能来自年份（降水），那么如何比较年份间的差异是否显著呢？
我用SAS的GLM过程语句如下，结果可以显示出肥料影响、年份影响及其交互作用（Ｆ值）是否显著，同时可以比较各年处理间是否显著。问题在于年份间差异如何显示（就像处理间比较时的显著性标注）？用当年的５个肥料处理数据作为重复进行方差分析似乎有问题，请教各位怎么才能实现，多谢多谢！！！
data SW;
input Nrate y1 y2 y3 y4 y5;
cards;
0 728 671 736 722 707
0 729 661 645 659 669
0 728 682 732 710 684
45 660 606 665 659 651
45 679 605 608 602 595
45 673 625 680 668 654
90 607 583 614 595 566
90 595 558 603 573 545
90 607 573 563 546 530
135 610 567 580 564 545
135 595 554 604 563 521
135 567 532 585 543 499
180 611 565 556 550 541
180 603 562 603 566 527
180 581 611 596 545 491
;
proc glm data=SW;
class Nrate;
model y1-y5=Nrate /nouni;
repeated time 5 /printe;
means Nrate /lsd;
run;
