Likelihood estimate test

[ktavabi]

Greetings,
Is it by any means straightforward to use Rj editor to compare likelihoods for two gamlj objects? E.g.

Code: Select all

null<-Y~ 1+x+(1|sid)
m1<-gamlj::gamljMixed(null...)
f<-Y~ 1+x+x1+(1|sid)
m2<-gamlj::gamljMixed(f...)
anova(m1, m2)

Thanks in advance,
Kam

[mcfanda@gmail.com]

sure, you can run this

Code: Select all

null<-Y~ 1+x+(1|sid)
m1<-gamlj::gamljMixed(null...)
f<-Y~ 1+x+x1+(1|sid)
m2<-gamlj::gamljMixed(f...)


anova(m1$model, m2$model)