我正在尝试确定各个季节之间站点之间的累积生长天数是否存在差异。相信我需要使用固定效果模型吗？我有5个站点和三个不同季节的数据，如下所示：
2 0017-07-01 10.724167 10.72417 spring
5 0017-07-02 9.899083 20.62325 fall
1 0017-07-03 10.782583 31.40583 spring
2 0017-07-04 11.719833 43.12567 spring
2 0017-07-05 11.534500 54.66017 spring
3 0017-07-06 11.845917 66.50608 fall
1 0017-07-07 11.653250 78.15933 spring
4 0017-07-08 10.768000 88.92733 winter
1 0017-07-09 10.584167 99.51150 winter
5 0017-07-10 12.258833 111.77033 spring
我正在使用看起来像这样的代码，但是我一直使用不同的迭代并获得不同的结果，却不知道它们的含义。即GDD〜网站+季节+网站*季节，或GDD〜网站+季节等。
results.lme1 <- lme(GDD~Site,random=~1|Season, data=data)
anova(results.lme1)
summary(results.lme1)
lmemodel2 <- confint(model1 <- glht(results.lme1, linfct = mcp(Site =
"Tukey")))
summary(lmemodel2)
任何帮助，指导和解释将不胜感激。谢谢！