I am a first year graduate student in biostatistics, and I have somewhat of an idea of the difference between REML and ML. However, I want a more in-depth understanding of each estimation method, especially when it comes to Likelihood ratio tests.
I know that you cannot use REML for likelihood ratio tests, and it has something to do with using the residuals of the means. If you take the residuals of the means, the models are no longer nested (nested models are needed for LRT). However, I really want a more thorough explanation as to why we cannot use REML for LRT's.
I will appreciate it if someone can explain this to me.