[R] anova for the objects of list

[Hai Lin]

Dear R experts,

Sorry for the previous email. I just figured it out
myself. I would welcome other solutions.

Thanks.
Kevin

my.anova <- sapply(names(my.x),
                               function(x)
                               {
                               anova(my.x[[x]],
my.x1[[x]])
                               }
                                 )

--- Hai Lin <kevinvol2002 at yahoo.com> wrote:

> Dear R experts,
> 
> I have two list objects from runing lme model for
> each
> subjects. And my.x is the full model with all terms
> I
> am interested in and my.x1 is the null model with
> only
> intercept.
> 
> I am trying to compare two models at each subject
> level. (Doing anova(my.x[[1]], my.x1[[1]]) is to
> produce stats for the first subject). My goal is to
> obtain stats. for an entire set with hundreds of
> subjects. I did something like this, lapply(my.Tem,
> function(x)anova(x)), but my.Tem actually is not a
> list after combining. It did not work the way as I
> expected.
> 
> I am very thankful if anyone here could point me
> out.
> 
> Kevin
> 
> ####
> 5 subjects extracted as my pilot data.
> > class(my.x)
> [1] "list"
> > length(my.x)
> [1] 5
> > class(my.x1)
> [1] "list"
> > length(my.x1)
> [1] 5 
> 
> ####
> my.Tem <- cbind(my.x, my.x1)
> > class(my.Tem)
> [1] "matrix"
> > my.Tem
>    my.x    my.x1  
> s1 List,15 List,15
> s2 List,15 List,15
> s3 List,15 List,15
> s4 List,15 List,15
> s5 List,15 List,15
> 
> 
> __________________________________________________
> Do You Yahoo!?

> protection around 
> http://mail.yahoo.com 
>