Re: need help with setting logistic mixed model

["Stella Karuri" <swkaruri@insightful.com>]

The argument ‘random’ can accept a list (See Pinheiro, J. C. and Bates., D. M. (2000). Mixed-effects Models in S and S-PLUS)

 

CSO.glme.out2<-glme(use~ethn+sex+exposure+jobduration+age+time+time^2+dept+jobgroup,

            Random = list ( id =  ~ dept + job – 1,  # exclude intercept to avoid over fitting

                                   plant= ~ 1),  

                                   data = "">

 

Stella

---

From: s-news-owner@lists.biostat.wustl.edu [mailto:s-news-owner@lists.biostat.wustl.edu] On Behalf Of Ali & Hind Lazrak
Sent: Tuesday, February 12, 2008 9:05 PM
To: s-news@lists.biostat.wustl.edu
Subject: [S] need help with setting logistic mixed model

 

Dear listserve users,

 

I need help with the following command line where I need to introduce a random effect for department (variable name: dept) and for job (variable name: job)at the subject (variable name: studyno1) level

I could not find any references on how to do that, so your help would be greatly appreciated

CSO.glme.out2<-glme(use~ethn+sex+exposure+jobduration+age+time+time^2+dept+jobgroup,random=~1| plant/id,data="" CSOgee$coefficients, method="REPQL")

 

Thank you so much in advance.