Re: [R] glmmBUGS: logistic regression on proportional data

2009-02-08 Thread David Winsemius
Just make sure it IS help. Following Dieter's advice I checked to see what format V&R use to pass data as success/failure and now wonder if the the proper invocation following their examples on pg 49, 59 and 134 in would be: mod<-glmmBUGS(

Re: [R] glmmBUGS: logistic regression on proportional data

2009-02-08 Thread Dieter Menne
David Winsemius comcast.net> writes: > >> > >> Newdat<-data.frame(Newtree=rep(1:3, each=20), Newsect=rep(c("a","b"), > >> each=10), Newdist=rep(1:5, 2), > >> y=rpois(60,2), tot=rep(c(14,12,10,8,6), 12)) > >> > >> yseed<-cbind(Newdat$y, Newdat$tot) > >> > >> mod<-glmmBUGS(yseed~N

Re: [R] glmmBUGS: logistic regression on proportional data

2009-02-08 Thread David Winsemius
Installing the glmmBUGS package and a bit of experimentation produces this minor modification of your code that seems to run without error. It's back to you to see if the output is sensible. library(glmmBUGS) Newdat<-data.frame(Newtree=rep(1:3, each=20), Newsect=rep(c("a","b"), each=10), Newd

Re: [R] glmmBUGS: logistic regression on proportional data

2009-02-08 Thread David Winsemius
On Feb 8, 2009, at 10:46 AM, Dieter Menne wrote: John Poulsen zoo.ufl.edu> writes: I am trying to run a logistic regression with random effects on proportional data in glmmBUGS. I am a newcomer to this package, and wondered if anyone could help me specify the model correctly. I am trying t

Re: [R] glmmBUGS: logistic regression on proportional data

2009-02-08 Thread Dieter Menne
John Poulsen zoo.ufl.edu> writes: > I am trying to run a logistic regression with random effects on > proportional data in glmmBUGS. I am a newcomer to this package, and > wondered if anyone could help me specify the model correctly. > > I am trying to specify the response variable, /yseed/,

[R] glmmBUGS: logistic regression on proportional data

2009-02-08 Thread John Poulsen
Hello, I am trying to run a logistic regression with random effects on proportional data in glmmBUGS. I am a newcomer to this package, and wondered if anyone could help me specify the model correctly. I am trying to specify the response variable, /yseed/, as # of successes out of total obse