[prev in list] [next in list] [prev in thread] [next in thread] 

List:       r-help
Subject:    [R] help with lme()
From:       "Bill Shipley" <bill.shipley () usherbrooke ! ca>
Date:       2003-11-04 15:42:45
[Download RAW message or body]

Hello. I am trying to determine whether I should be using ML or REML
methods to estimate a linear mixed model.   In the book by Pinheiro &
Bates (Mixed-effects models in S and S-PLUS, page 76) they state that
one difference between REML and ML is that « LME models with different
fixed-effects structures fit using REML cannot be compared on the basis
of their restricted likelihoods.  In particular, likelihood ratio tests
are not valid under these circumstances.”

I am not sure what “fixed-effects structures” means.  Does it mean that,
as long as the types of contrasts are the same between two models, they
ARE comparable, but are NOT comparable if the types of contrasts are
changes?  Or rather, does it simply mean that one should use t or F
tests for the fixed effects, and restrict the likelihood ratio tests to
the random effects only if using REML?

 

Bill Shipley

Associate Editor, Ecology

North American Editor, Annals of Botany

Département de biologie, Université de Sherbrooke,

Sherbrooke (Québec) J1K 2R1 CANADA

Bill.Shipley@USherbrooke.ca

 <http://callisto.si.usherb.ca:8080/bshipley/>
http://callisto.si.usherb.ca:8080/bshipley/

 


	[[alternative HTML version deleted]]


[prev in list] [next in list] [prev in thread] [next in thread]