MCMCglmm prior specification
?? Ok, thank you very much. I am actually asking all this because I am trying to write a very quick guide to MCMCglmm so that people can achieve in a pair of hours what cost me more than two weeks, this is, being able to specify an extremely simple model that at least runs and gives more or less reasonable results, including things like categorical data and bivariate responses. Obviously I don't want to do it too complicated (I wouldn't be able anyway, I have so much to learn myself yet!), so I suppose I will just give a pair of examples of uninformative priors for univariate and multivariate responses.?? I will post the guide here when it's finished so that you guys can point at the multiple errors it will surely have. ?? Thank you very much,?? Iker ?__________________________________________________________________ ?? Iker Vaquero-Alba ?? Visiting Postdoctoral Research Associate ?? Laboratory of Evolutionary Ecology of Adaptations ?? Joseph Banks Laboratories ?? School of Life Sciences ?? University of Lincoln?? Brayford Campus, Lincoln ?? LN6 7DL ?? United Kingdom ?? https://eric.exeter.ac.uk/repository/handle/10036/3381 De: Jarrod Hadfield <j.hadfield at ed.ac.uk> Para: Iker Vaquero Alba <karraspito at yahoo.es> CC: "r-sig-mixed-models at r-project.org" <r-sig-mixed-models at r-project.org> Enviado: Lunes 12 de octubre de 2015 12:40 Asunto: Re: [R-sig-ME] MCMCglmm prior specification Hi, Usually people think the default priors for the fixed effects (zero? mean, high variance) are reasonable. However, there are cases where? stronger priors are useful. For example, a) you might actually have? some prior information or b) you might have near or complete? separation in a GLMM (usually with categorical data) and you might? want to constrain the fixed effects so they don't result in extreme? predictions. Cheers, Jarrod ? Quoting Iker Vaquero Alba <karraspito at yahoo.es> on Mon, 12 Oct 2015? 11:05:58 +0000 (UTC):
?? Hi Jarrod,?? Thanks a lot for your reply. I had actually found? information about it in another post and wanted to "fix" my own? question, but for some reason individual posts don't arrive to my? email even if I have that option activated, so I couldn't reply to? myself.?? My question, then, would be: most of the times I see? priors with just R and G elements, but not B. Is it not necessary to? specify B? Also, if my model does not have random effects, then? would it be enough with a prior for R element? ?? Thank you very much.?? Iker.
__________________________________________________________________ ?? Iker Vaquero-Alba ?? Visiting Postdoctoral Research Associate ?? Laboratory of Evolutionary Ecology of Adaptations ?? Joseph Banks Laboratories ?? School of Life Sciences ?? University of Lincoln?? Brayford Campus, Lincoln ?? LN6 7DL ?? United Kingdom ?? https://eric.exeter.ac.uk/repository/handle/10036/3381 ? ? ? De: Jarrod Hadfield <j.hadfield at ed.ac.uk> ? Para: Iker Vaquero Alba <karraspito at yahoo.es> CC: "r-sig-mixed-models at r-project.org" <r-sig-mixed-models at r-project.org> ? Enviado: Lunes 12 de octubre de 2015 11:56 ? Asunto: Re: [R-sig-ME] MCMCglmm prior specification Hi, R defines the prior for the residual covariance matrices. B specifies? the prior for the `fixed' effects. Cheers, Jarrod Quoting Iker Vaquero Alba <karraspito at yahoo.es> on Mon, 12 Oct 2015? 10:34:01 +0000 (UTC): ?? Hello everyone, ?? I am reading a lot of documentation at the moment about prior? specification in MCMCglmm, and there is still something I have not? very clear. According to some sources, there are mainly two elements? to take into account when defining a prior:?? - An R structure that? needs to be specified for each fixed effect. And ?? - A G structure for each random effect. ?? However, in other sources another element, B, is introduced as? well, which refers to fixed effects too. In Jarrod's Course Notes, I? see that both R and B have to do with fixed effects: R specifies V? and nu arguments for the variance, and B specifies V and mu elements? for the mean. ?? ?? My confusion comes from the fact that almost every time I see an? example of a prior, it just has two elements, R (for fixed effects)? and G (for random effects). Why is this? Is it not so important to? define a prior for the mean? Is it enough with a prior specification? for the variance of fixed effects and another one for all the random? effects? ?? Thank you very much in advance.?? Iker ? __________________________________________________________________ ?? Iker Vaquero-Alba ?? Visiting Postdoctoral Research Associate ?? Laboratory of Evolutionary Ecology of Adaptations ?? Joseph Banks Laboratories ?? School of Life Sciences ?? University of Lincoln?? Brayford Campus, Lincoln ?? LN6 7DL ?? United Kingdom ?? https://eric.exeter.ac.uk/repository/handle/10036/3381 ??? [[alternative HTML version deleted]] _______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.
The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. [[alternative HTML version deleted]]