[R-sig-ME] glm model with all zeros for one of the factor level
Juan Pablo Edwards Molina
edwardsmolina at gmail.com
Tue Dec 5 16:58:00 CET 2017
Dear List members:
I performed two independent experiments (CRD) to test if a whitefly
has preference to infect with a virus: potatos, tomatos or peppers
(target hosts, TH), wether if the virus was obtained from potato or
tomato (source hosts, SH). So I released 100 white flyes (previously
infected with the virus from one or other SH) inside cages containing
10 plants of each TH (30 total). This is how the data looks like:
exp SH TH cage tot posit
1 tom tom 1 10 4
1 tom bat 1 10 3
1 tom pep 1 10 0
1 bat tom 2 10 1
1 bat bat 2 10 2
1 bat pep 2 10 0
2 tom tom 3 10 6
2 tom bat 3 10 4
2 tom pep 3 10 0
2 bat tom 4 10 4
2 bat bat 4 10 0
2 bat pep 4 10 0
The issue I found here is that pepper was not infected at all,
however it was infected in another experiment without chance of TH
choice: i.e. I released infectious whiteflies inside cages containing
the same pepper genotyope and they present the typical virus disease
symptoms.
So, how should I consider modeling this data?
Zero-inflated negative binomial using the total plants as offset? Hurdle-model?
Should I remove the pepper level for the model?
Any help would be really helpful.
Juan Edwards
More information about the R-sig-mixed-models
mailing list