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]]

