Displaying 1 result from an estimated 1 matches for "nvect".
Did you mean:
nvec
2007 Jun 01
2
Interaction term in lmer
...models:
Model 1: center effect (random) only
Model 2: trt (fixed) + center (random) + trt*center interaction.
Then, I want to compare these 2 models with Likelihood Ratio Test. Here are
my lmer codes that I don't feel comfortable about their correctness.
model1 <- try(lmer(cbind( yvect, nvect-yvect) ~ 1 + (1 | center),
family = binomial, niter = 25, method = "Laplace", control = list(usePQL =
FALSE) ))
model2 <- try(lmer(cbind( yvect, nvect-yvect) ~ trt*center + ( 1 | center)
,
family = binomial, niter = 25, method = "Laplace", control = list(usePQL =
FAL...