Good point about statistical software. This can lead to a lot of confusion in comparing output from different programs with the same data and same model. In MIXED and GLIMMIX, however, one does get the complete likelihood (or restricted likelihood) in the information-criterion statistics, including the term(s) that does not depend on the the parameters. For REML estimation, only variance-covariance parameters are used in the calculation because the fixed effect are eliminated at each ML step. I hope to see a compete output from the OP to help figure out this issue.
... View more
