[R] Correct nested design for GLM
Hi, I am currently running the following negative binomial GLM: glm89.nb <- glm.nb(AvGUD ~ Year*Trt*Micro + (0 + Micro/Trt/Year)) where Year has 3 levels, Trt has 2 levels, and Micro has 3 levels. >From what I have read the above model has a 3 way interaction (Year*Trt*Micro), and Micro is nest within Trt and Trt is nested with Year (0 + Micro/Trt/Year). I was hoping someone could confirm that this is actually what the model is running? Thank you very much for your time, Rebecca [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] Zero inflated negative binomial
Hi all, I am running the following model: > glm89.nb <- glm.nb(AvGUD ~ Year*Trt*Micro) where Year has 3 levels, Trt has 2 levels and Micro has 3 levels. However when I run it has a zero inflated negative binomial (as I have lots of zeros) I get the below error message: > Zinb <- zeroinfl(AvGUD ~ Year*Trt*Micro |1, data = AvGUD89, dist = "negbin") Error in optim(fn = loglikfun, gr = gradfun, par = c(start$count, start$zero, : non-finite value supplied by optim For what I have read I think the problem is that for Year level 3 there is no Trt 1 and for Year level 1 there is no Micro level 3. I cannot find a solution to this problem, is there any way I can solve this so I can run the zero inflated model? Thank you for your time Rebecca [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] LMM and unequal variances
Hi all, I am running a LMM and have two fixed factors: Trt (with 2 levels) and Microhabitat (with 3 levels), and a randon factor: Group nested within Trt. I have unequal variances between the two Trt levels. Does this matter in a LMM? Thanks for your time Rebecca [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
