[ADMB Users] glmmadmb in R
Mollie Brooks
mbrooks at ufl.edu
Thu Nov 10 13:49:20 PST 2011
This is an example of the trouble that can happen when the same information is in multiple locations and people (understandably) only have the energy to update the information in a subset of the locations.
I guess the information at http://admb-project.org/examples/r-stuff/glmmadmb is not current. I'll try to fix this.
Mollie
On 10 Nov 2011, at 4:42 PM, Mollie Brooks wrote:
> Hi Chris,
> I get the same error as you when I run
> fit = glmmadmb(y~Base*trt+Age+Visit,random=~Visit,group="subject",data=epil2,family="nbinom")
> but this works
> fit = glmmadmb(y~Base*trt+Age+Visit+(Visit|subject),data=epil2,family="nbinom")
>
> So I guess you have the newer version of glmmADMB (like me) and I guess it requires this syntax.
> cheers,
> Mollie
>
> Mollie Brooks
> Ph.D. Candidate
> NSF IGERT Fellow
> Biology Department
> University of Florida
> mbrooks at ufl.edu
> http://people.biology.ufl.edu/mbrooks
>
>
>
>
> On 10 Nov 2011, at 4:01 PM, Christopher Smith wrote:
>
>> Hello,
>>
>> I am fairly new to R but am attempting to use GLMMs with neg. bin. distribution to analyze count data. (I have successful done this using glmer() and glm.nb() ). I am having trouble getting the function glmmabmb() to work, even for the example provided at: http://admb-project.org/examples/r-stuff/glmmadmb
>>
>> I am getting the error,
>>
>> > fit = glmmadmb(y~Base*trt+Age+Visit,random=~Visit,group="subject",data=epil2,family="nbinom")
>> Error in glmmadmb(y ~ Base * trt + Age + Visit, random = ~Visit, group = "subject", :
>> unused argument(s) (group = "subject")
>>
>>
>> And/or
>>
>> > fit = glmmadmb(y~Base*trt+Age+Visit,random=~Visit,data=epil2,family="nbinom")
>> Error in `colnames<-`(`*tmp*`, value = character(0)) :
>> attempt to set colnames on object with less than two dimensions
>>
>> The second error is the same error I get when trying to use the function with my data.
>>
>> E.g.,
>>
>> > glmmadmb(ANAM~tmin+CWD,family="nbinom",random=~Day)
>> Error in `colnames<-`(`*tmp*`, value = character(0)) :
>> attempt to set colnames on object with less than two dimensions
>>
>>
>> I am on a Mac (OS X 10.5.8) running R (2.12- 32bit). Any advice?
>>
>> Thanks,
>> Chris
>> _______________________________________________
>> Users mailing list
>> Users at admb-project.org
>> http://lists.admb-project.org/mailman/listinfo/users
>
> _______________________________________________
> Users mailing list
> Users at admb-project.org
> http://lists.admb-project.org/mailman/listinfo/users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20111110/4476739e/attachment.html>
More information about the Users
mailing list