search for: model_11

...mating mu as the fitted model like so: model_1 <-glm.nb(y~x1+x2+x3, data = datafilename) mu_1 <- fitted(model_1) theta_1 <- theta.ml(y, mu_1, length(data), limit = 10, eps = .Machine$double.eps^0.25, trace = FALSE) Then, I conduct the lmer, using the estimated theta: model_11 <-lmer(y~x1+x2+(1|x3), family = negative.binomial(theta = theta_1, link = “log”), method = “Laplace”) First, I wondered if this sounds like a reasonable method to accomplish my goals. Second, I wondered if the theta I use for reduced models (nested within model_11) should be estima...