[R-sig-ME] mcmcsamp in lme4

[Robert A LaBudde]

At 10:55 AM 6/9/2009, De Smedt Sebastiaan wrote:
>I fitted a mixed-effects model with a Variance Structure in nlme
>(because Variance Structures are not yet implemented in lme4). After
>several model fittings (with method="ML"), I found the optimal fixed
>structure of my explanatory variables.
>
>I understood from this mailing list that the p-values of the fixed
>effects, given by summary(namemodel) are not trustable and thus I want
>to use a Monte Carlo simulation in order to estimate the significance of
>the fixed terms. This command doesn't exist in nlme, and I can't fit my
>model in lme4. Is there a way to overcome this problem? Or another
>function in nlme that gives comparable results?

Try intervals(fit) in 'nlme'.

================================================================
Robert A. LaBudde, PhD, PAS, Dpl. ACAFS  e-mail: ral at lcfltd.com
Least Cost Formulations, Ltd.            URL: http://lcfltd.com/
824 Timberlake Drive                     Tel: 757-467-0954
Virginia Beach, VA 23464-3239            Fax: 757-467-2947

"Vere scire est per causas scire"