Re: [R] Why are coefficient estimates using ML and REML are different in lme?

2012-10-29 Thread S Ellison
Yi) Different criteria would be _exepcted_ to give different estimates. ii) Look at teh standard errors on the coefficients. Essentially all of them are larger than the estimates for both fitting criteria. Essentially, both models are telling you that your coefficients are not significantly

Re: [R] Why are coefficient estimates using ML and REML are different in lme?

2012-10-29 Thread Houhou Li
Thank you so much, S. Ellison.   Maybe I am wrong since I am still learning the mixed model, my impression is that variance and standard error estimates from REML and ML are different, coefficient estimates for fixed effects from both ML and REML are not necessarily identical, but should be

[R] Why are coefficient estimates using ML and REML are different in lme?

2012-10-28 Thread Houhou Li
Hi, All,   My data collection is from 4 regions (a, b, c, d). Within each region, it has 2 or 3 units. Within each unit, it has measurement from about 25 sample site. I was trying to use lme function to discribe relationship between y and a few covariates. Both y and covariates were measured at