Skip to content

broken example: lme() + multcomp() Tukey on repeated measures design

4 messages · Dieter Menne, David Winsemius, William Simpson

William Simpson
# 
I am trying to do Tukey HSD comparisons on a repeated measures expt.
I found the following example on r-help and quoted approvingly elsewhere.
It is broken. Can anyone please tell me how to get it to work?

I am using R 2.4.1.
Error in eval(expr, envir, enclos) : object "Y" not found
Error in factor_contrasts(model) : no 'model.matrix' method for 'model' found!

It all went ok until the last line, where the actual Tukey multiple
comparisons was requested.

Thanks very much for any help!

Bill
Dieter Menne
# 
William Simpson <william.a.simpson <at> gmail.com> writes:
That's the problem. Not the R-Version, but the likelihood that the
multcomp package was not updated also, and lme definitively not
in early version.s



other attached packages:
[1] multcomp_1.0-7 mvtnorm_0.9-5  nlme_3.1-90    MASS_7.2-46
General Linear Hypotheses

Multiple Comparisons of Means: Tukey Contrasts


Linear Hypotheses:
                              Estimate
Marvellous - Golden.rain == 0    5.292
Victory - Golden.rain == 0      -6.875
Victory - Marvellous == 0      -12.167



Dieter
David Winsemius
#
On Apr 21, 2009, at 10:12 AM, William Simpson wrote:

            

      
      
      
    

        
How many years have you had that version without updating?

            

      
      
      
    

        
It is not broken when run in version 2.8.1.


David Winsemius, MD
Heritage Laboratories
West Hartford, CT

  
  


  


      

    

    

      
      

        


  

        

      William Simpson
    

    

      
      #
    

  


  

      
Thanks for the helpful replies.
I will try updating R and see how it goes.

MASS, nlme, multcomp are all the most recent.

Bill

On Tue, Apr 21, 2009 at 3:12 PM, William Simpson

        
<william.a.simpson at gmail.com> wrote: