[R-sig-ME] extracting p values for main effects of binomial glmm
Ben Bolker
bbolker at gmail.com
Wed Mar 4 23:46:26 CET 2015
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
On 15-03-04 05:09 PM, Doogan, Nathan wrote:
> At last check, the summary() of a glmer() object does report
> p-values. It's output from lmer() that does not. I suppose I could
> be working with an older version...
>
> -Nate
>
Yes (and has more or less for forever), although these are based on
Wald Z-scores. Thanks to the reminder from Jonathan Baron about
?pvalues ...
>
>
> -- Nathan J. Doogan, Ph.D. Post-Doctoral Researcher College of
> Public Health The Ohio State University
>
> ________________________________________ From: R-sig-mixed-models
> [r-sig-mixed-models-bounces at r-project.org] on behalf of Megan
> Kutzer [makutzer at gmail.com] Sent: Wednesday, March 04, 2015 3:11 PM
> To: r-sig-mixed-models at r-project.org Subject: [R-sig-ME] extracting
> p values for main effects of binomial glmm
>
> Hi,
>
> I'm fairly new to mixed models and have done a lot of reading
> without much success. Unfortunately there is no one at my
> institution who is really familiar with them so I thought I would
> try this list.
>
> I'm running a binomial generalized linear mixed effects model and
> I need p-values for the main effects. I know this isn't entirely
> correct with this type of model but my supervisor wants the
> p-values!
>
> The model is:
>
> glmer (Proportion hatched ~ Diet * Infection status * Day +
> (1|SubjectID) + (1|Day), family=binomial)
>
> where,
>
> Proportion hatched = cbind(Offspring, Eggs-Offspring) Diet is a
> factor with 2 levels Infection status is a factor with 4 levels
> Day is a factor with 3 levels
>
> Using Subject ID number and Day as random effects is supposed to
> control for pseudoreplication in the model, although I am not
> entirely sure that this is specified in the correct way. I wanted
> to include experimental replicate here too but the model failed to
> converge.
>
> My question is: is there a way to get p-values for the main fixed
> effects of Diet, Infection and Day?
>
> If you need more specific model information or the model output I
> would be happy to provide it.
>
> Thanks, Megan
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.11 (GNU/Linux)
iQEcBAEBAgAGBQJU94tCAAoJEOCV5YRblxUHnPcIAKgay6VpUxGahQfi7JMsj43Z
iDPdR2d97ANFLqOXaZJpGktquGKtFbYGmua9XSyz4DSdl5bpX9n2hb+NYT4Secdo
1cVc28noVbP83Od6wHdJgyRWmzJ85onZ+m0t8CK2BtnJO4C/1GMUZwX9H4+F8g6N
ljbIaAD07AOmx14VF8l4XpW9guw3NNEOgKcupZBJ5K7zsLkkdtnLoAJZFCcxpypg
WdGcrcft82nJorqwcEYrjpqF/X9CzFvRIPAhf+cbZbDe06W98coM4RVNyLRXrLNE
2NKOKK2Pg5tHVAW6PnlaOBBuSCFgpSeWTQINDLda/Q3j5OUqmoCeLjpjjTb3A30=
=iWlw
-----END PGP SIGNATURE-----
More information about the R-sig-mixed-models
mailing list