[R-sig-ME] p-value for glmer with quasipoisson family
CL Pressland
Kate.Pressland at bristol.ac.uk
Tue Sep 29 15:52:38 CEST 2009
Please check the archive - there is a lot written on here explaining why
lme4 does not give p values, unless you're using a MCMC approach which is
better able to estimate the parameters.
In particular see:
https://stat.ethz.ch/pipermail/r-sig-mixed-models/2009q3/002833.html
--On 29 September 2009 15:09 +0200 Raldo Kruger <raldo.kruger at gmail.com>
wrote:
> Hi all,
>
> How come there is no p-value displayed in the summary of the results
> of glmer when using the quasipoisson family? Is there a way of
> extracting the p-values from the results?
>
> Thanks,
> Raldo
>
>> summary(ex5m_r)
> Generalized linear mixed model fit by the Laplace approximation
> Formula: Plants ~ Treat + Year + (1 | Site)
> Data: ex5m
> AIC BIC logLik deviance
> 4685 4719 -2336 4671
> Random effects:
> Groups Name Variance Std.Dev.
> Site (Intercept) 3.550 1.8841
> Residual 12.957 3.5996
> Number of obs: 928, groups: Site, 8
>
> Fixed effects:
> Estimate Std. Error t value
> (Intercept) 1.02627 0.68130 1.506
> TreatGr -0.38019 0.17078 -2.226
> TreatK -0.84274 0.19838 -4.248
> TreatN 0.04644 0.15215 0.305
> Yeartwo 0.66848 0.12931 5.170
>
> Correlation of Fixed Effects:
> (Intr) TretGr TreatK TreatN
> TreatGr -0.102
> TreatK -0.088 0.350
> TreatN -0.114 0.456 0.392
> Yeartwo -0.125 0.000 0.000 0.000
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
More information about the R-sig-mixed-models
mailing list