Greetings, I have a question that I'd like to get input on. I have a classic toxicology study where I artificially fertilized and exposed embryos to a chemical and counted defects. In addition, I kept track of male-female pairs that I used to artificially fertilize and generate embryos with. I need to use logistic regression to model the response, but also check that the genetics of the pairings did or did not have an effect on the response. My data looks a bit like this:
response matrix chemical concentration Genetic Matrix Present Absent Male Female 2 152 0.13 a 1 6 121 1 a 2 21 92 2 b 3 24 89 5 b 4 0 141 10 c 5 5 95 15 c 6 R code: DA<-cbind(Present, Absent) glm<-(DA ~ chemical concentration) If I do glm<-(DA ~ chemical concentration + Male + Female, I get every possible combination, but I only want specific pairs. So, I am thinking about doing: MF<-cbind(Male, Female) glm<-(DA ~ chemical concentration + MF) Is this correct? Any help on how to model this would be greatly appreciated! Thanks in advance! [[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.