Re: [R] specification for glmmPQL

2005-09-14 Thread Andrew R. Criswell
Dear Prof. Bates and Group: I hope it is not to late to revisit this thread. My concern is with the difference in standard errors estimated from data that is arranged as grouped (data.1) and ungrouped (data.2). With the grouped data set, the effect of treatment is highly significant; with

[R] specification for glmmPQL

2005-09-04 Thread Andrew R. Criswell
Hello All, I have a question regarding how glmmPQL should be specified. Which of these two is correct? summary(fm.3 - glmmPQL(cbind(response, 100 - response) ~ expt, data = data.1, random = ~ 1 | subject, family = binomial)) summary(fm.4 -

Re: [R] specification for glmmPQL

2005-09-04 Thread Andrew R. Criswell
Hello Dr. Bates and group, I understand, the attached data file did not accompany my original message. I have listed below the code used to create that file. data.1 - data.frame(subject = factor(rep(c(one, two, three, four, five, six, seven,

Re: [R] specification for glmmPQL

2005-09-04 Thread Douglas Bates
On 9/4/05, Andrew R. Criswell [EMAIL PROTECTED] wrote: Hello Dr. Bates and group, I understand, the attached data file did not accompany my original message. I have listed below the code used to create that file. data.1 - data.frame(subject = factor(rep(c(one, two, three, four,