[R-sig-ME] broom package and MCMCglmm estimates

dani orchidn at live.com
Tue Oct 24 23:56:49 CEST 2017


Hello Pierre,


Thank you so much for your prompt response, I will follow your advice and report back.


best,

daniNm

________________________________
From: Pierre de Villemereuil <pierre.de.villemereuil at mailoo.org>
Sent: Tuesday, October 24, 2017 1:47 PM
To: r-sig-mixed-models at r-project.org
Cc: dani
Subject: Re: [R-sig-ME] broom package and MCMCglmm estimates

Hi,

As far as I can see, there is no method to tidy up MCMCglmm output in broom. You can see that by looking at the "tidiers" functions in broom using:
library(help = broom) #No MCMCglmm

However, MCMCglmm output contains "mcmc" class objects that can be tidied up using tidyMCMC():
- fixed effects: tidyMCMC(model$Sol)
- random effects: tidyMCMC(model$VCV)

Look up ?tidyMCMC for more options. I'm afraid it's difficult to get the pMCMC back though. You can get 95% credible intervals are a replacement for those (set conf.int = TRUE).

I have also seen that Ben Bolker is working on a fork/split of broom for mixed models especially, maybe he plans on having a specific tidier for MCMCglmm?

Best,
Pierre.


On Tuesday, 24 October 2017 20:13:53 NZDT dani wrote:
> Hi everyone,
>
>
> I would like to visualize the estimates from a MCMC glmm model and I tried using broom but I received the following message:
>
> Error in as.data.frame.default(mj) :
>   cannot coerce class ""MCMCglmm"" to a data.frame
> I was wondering what should I do about this?
> Thank you in advance for all your help!
> Best,
> DaniNM
>
>        [[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 Info Page - stat.ethz.ch<https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models>
stat.ethz.ch
Your email address: Your name (optional): You may enter a privacy password below. This provides only mild security, but should prevent others from messing ...



>


	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list