[R] lmer and estimation of p-values: error with mcmcpvalue()

[Gregor Gorjanc]

Christoph Scherber <Christoph.Scherber <at> agr.uni-goettingen.de> writes:
> 
> Dear Henric,
> 
> Thanks, now it works; but how reliable are these estimates? Especially 
> with p-values close to 0.05 it is of course important that the range of 
> the estimates is not too large. I´ve just run several simulations, each 
> of which yielding sometimes quite different p-values.

I have not been using mcmcsamp(), but try increasing number of simulations.
Additionally, since you have only "single" parameters, you can directly look
at their posterior distributions and compute the areas under the curve.

Gregor