[R] Needs suggestions for choosing appropriate R packages

From: zhijie zhang <epistat_at_gmail.com>
Date: Wed, 26 Nov 2008 21:29:28 +0800


Dear all,
  I am thinking to fit a multilevel dataset with R. I have found several possible packages for my task, such as
glmmPQL(MASS),glmm(repeated),glmer(lme4), et al. I am a little confused by these functions.
  Could anybody tell me which function/package is the correct one to analyse my dataset?
My dataset is as follows:
the response variable P is binary variable (the subject is a patient or not);
two explanatory variables X1 (age) and X2 (sex); this dataset was sampled from three different level, district,school,individual, so this was regarded as a multilevel dataset. I hope to fit the 3-level model(Y is binary variable): Logit(Pijk)=(a0+b0k+c0jk)+b1*X1+b2*X2
  i-individual, first level; j-school, 2nd level; k-district,3rd level.  I know that the GLIMMIX procedure in the latest version SAS9.2 is a choice for that, but unfortunately we donot have the latest version.  R must have the similar functions to do that, can anybody give me some suggestions or help on analysing my dataset? Q1: Which package/functions is appropriate for my task? Could u show me some example codes if possible?
Q2: Logit(Pijk)=(a0+b0k+c0jk)+(b1+b1j)*X1+b2*X2   If the randome effect was also specified in X1 as above, Which package/functions is possible?
  Thanks a lot.

-- 
With Kind Regards,

oooO:::::::::
(..):::::::::

:\.(:::Oooo::
::\_)::(..):: :::::::)./::: ::::::(_/:::: ::::::::::::: [***********************************************************************] ZhiJie Zhang ,PhD Dept.of Epidemiology, School of Public Health,Fudan University Office:Room 443, Building 8 Office Tel./Fax.:+86-21-54237410 Address:No. 138 Yi Xue Yuan Road,Shanghai,China Postcode:200032 Email:epistat_at_gmail.com <Email%3Aepistat_at_gmail.com> Website: www.statABC.com [***********************************************************************] oooO::::::::: (..):::::::::
:\.(:::Oooo::
::\_)::(..):: :::::::)./::: ::::::(_/:::: ::::::::::::: [[alternative HTML version deleted]] ______________________________________________ R-help_at_r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Received on Wed 26 Nov 2008 - 15:14:18 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 26 Nov 2008 - 15:30:27 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive