Hi Brigitte, >Did somebody know why asreml does not provide the same REML loglikehood >as coxme, lme4 or lmne.
I don't know the answer to this, but I'd guess it is either to do with the use of the average information REML algorithm or asreml-r is for some reason ending up with a different subset of the data. >If it was just a constant value between the two models (with or without >the fixed effect) it would not be important. But it is not. >I checked that the variance component estimators were equal. I'm still not clear that it is important (if the data subset analysed is the same). You would only use the REML likelihoods to compare models with different random effects and the same fixed effect structure (is there another use for the REML likelihood other than that?), so then it is really a question of whether for a given pair of random effect models and the same data the likelihood ratio test statistic changes across analysis methods. Unless for some reason you are comparing two random effect models fitted with different routines (one of which is asreml-r). Ron. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
