On Thu, 26 Aug 2004, Deepayan Sarkar wrote:

> Quoting David Scott <[EMAIL PROTECTED]>:
> 
> > 
> > I am trying to create a plot similar to Figure 3.2 in Bates and Pinheiro.
> > 
> > I have repeated measurements on about 80 subjects from 2 treatment groups.  
> > I would like to have the panels for the two treatment groups in separate
> > groups and within those groups have the panels ordered on maximum value
> > (as is the default).
> > 
> > I am ok with getting plots similar to Figs 3.1 and 3.2, but can't see how 
> > to change the ordering of the panels to what I want.
> > 
> > Here is the definition of my groupedData object
> > 
> > RAWlmeData <- groupedData(RAW~Elapsed|ID,
> >               data=RAWData,
> >               labels=list(x="Elapsed time",y="Airways resistance"),
> >               units=list(x="(hours)",y="cm H20/L/sec"))
> > 
> > I guess I could just plot the two treatment groups separately in turn but 
> > I feel there is something I am missing.
> 
> My reading of page 105 suggests that you need to specify 
> 'outer=<whatever your grouping factor is>'.
> Have you tried that?
> 
I did try that. It produces Fig 3.3 with multiple lines giving the results 
for each subject by treatment. With around 40 subjects per treatment 
group, that isn't very attractive in my case.

I should also mention that I looked for on-line complements to Pinheiro 
and Bates, and found MEMSS.tar.gz on 
http://cm.bell-labs.com/cm/ms/departments/sia/NLME/MEMSS/index.html
which promised scripts for the examples in the book but turned out to have 
only chapters 1 and 2, and for S-PLUS rather than R (going on the .q 
extension).

David Scott

_________________________________________________________________
David Scott     Department of Statistics, Tamaki Campus
                The University of Auckland, PB 92019
                Auckland        NEW ZEALAND
Phone: +64 9 373 7599 ext 86830         Fax: +64 9 373 7000
Email:  [EMAIL PROTECTED] 


Graduate Officer, Department of Statistics

______________________________________________
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Reply via email to