Re: [R] help with winbugs glm

2012-02-24 Thread ilai
My apologies to R-list members, this discussion diverged and is no longer R related. Probably more fitting in one of the BUGS forums, but I keep cc ing for any future interested reader that stumbled upon this post in r-help. Now to the point. Thank you so much for the reference, I was not aware of

Re: [R] help with winbugs glm

2012-02-24 Thread Adan Jordan-Garza
Ohh! about the reference code for the categorical predictor, from the McCarthy book, Bayesian Methods for Ecology (ISBN978-0-521-85057-5) the following ANCOVA model has a Refuge as a categorical predictor with 3 levels coded as "1", "2" and "3". This is the cited model: model { a ~ dnorm(1.01, 24

Re: [R] help with winbugs glm

2012-02-24 Thread Adan Jordan-Garza
Ok this makes a lot of sense, thank you very much Ilai! Cheers Guillermo On Fri, Feb 24, 2012 at 12:16 AM, ilai wrote: > On Thu, Feb 23, 2012 at 8:32 PM, Adan Jordan-Garza > wrote: > > Hello Ilai, > > thank you very much for your response, > > can I bother you a little further? > > What do you

Re: [R] help with winbugs glm

2012-02-23 Thread ilai
Adan, How many levels does Depth have? my wild guess: 3 and your bugs model is not identifiable. Second, I think you may have a critical error in the way you formatted the data for the bugs model. From your code it looks like you are just using the factor Depth and not a design matrix of dummy var

[R] help with winbugs glm

2012-02-23 Thread Adan Jordan-Garza
Hi, I am running a model with count data and one categorical predictor (simple model for me to understand it fully), I did in R a glm like this: glm(Recruitment~Depth, family=poisson). I get the coefficientes and confidence intervals and all is ok. But then I want to do the same model with Bayesian