Try using mcmcsamp() to sample from the posterior distribution of the parameter estimates. You can calculate a p-value from that, if that is your desire. Instructions are in the R wiki: http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests
HTH, Simon. Dan Bebber wrote: > Hello, > I've just located the illuminating explanation by Douglas Bates on degrees > of freedom in mixed models. > The take-home message appears to be: don't trust the p-values from lme. > Questions: > Should I give up hypothesis testing for fixed effects terms in mixed models? > Has my time spent reading Pinheiro & Bates been in vain? > Is there a publication on this issue? > > Thanks, > Dan Bebber > > Department of Plant Sciences > University of Oxford > > ______________________________________________ > R-help@stat.math.ethz.ch mailing list > 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. > > -- Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat. Centre for Resource and Environmental Studies The Australian National University Canberra ACT 0200 Australia T: +61 2 6125 7800 email: Simon.Blomberg_at_anu.edu.au F: +61 2 6125 0757 CRICOS Provider # 00120C The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. - John Tukey. ______________________________________________ R-help@stat.math.ethz.ch mailing list 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.