I completly changed my answer as a result of a long conversation with Daniel. I'll try to provide some background information so interested readers can understand my answer.
As I understand the question, Daniel is trying to assess the effect of no.Green on the probability of subjects choosing red in an experiment. no.Green is the number of balls (1, 2, 3 or 4). And the experiment was conducted under several conditions, nameley: under no.Red equals to 5, 7 and 9; and under condition in which message is blue and message is red. So, we have a total of 4 * 3 * 2 = 24 conditions (4 conditions from no.Green, 3 conditions from no.Red and 2 conditions from message blue or red).
One single regression with all interactions terms is quite complex to interpret. However, his main task is fairly simple, namely: to show that no.green has an effect on the probability of choosing red. So, my sugestion is to run a separece regression for message == blue and message == red condition, and also a separate regression for each no.red condition. Moreover, I'll simplify thing by assuming that no.Green is continous (it seems to be possible treat it as continous, or at least an interval variable). In R, for the message == blue case, just do this:
fit.1 = glm(DecisionasReceiver ~ no.GREEN, family=binomial, data=subset(lue, messagereceived=="blue" & no.RED==5) )
fit.2 = glm(DecisionasReceiver ~ no.GREEN, family=binomial, data=subset(lue, messagereceived=="blue" & no.RED==7) )
fit.3 = glm(DecisionasReceiver ~ no.GREEN, family=binomial, data=subset(lue, messagereceived=="blue" & no.RED==9) )
Now, in order to properly asses the effect of no.green, you have to take in consideration the uncertainty of estimates. Looking to standard errors, you will see that no.green is significant. However, looking only for standard errors does not allow you to proper understadn the range of uncertainty. Say, for instance, that you are interested to know how less likely to choose red are subjects (under condition no.red ==5) with no.green == 2 compared with subjects with no.green == 1. To answer this type of question, it's better, I think to look at predicted probability, but taking into consideration the uncertainty on estimates. To do this, i'll use the "sim" function, of arm package.
n.sims = 1000
sim.1 = sim(fit.1, n.sims)
with(subset(lue, messagereceived=="blue" & no.RED==5), plot(no.GREEN,jitter(DecisionasReceiver, .1),
ylab="Probability of Choosing Red", xlab="Number of Green",
main="Effect of Green under no.Red equals 5"))
for (s in 1:100)
curve(invlogit(coef(sim.1)[s,1] +coef(sim.1)[s,2]*x), col="gray", xlim=c(1,4), add=T)
The result is a graphic with 100 logistic curves. Each curve represents a possible effect of no.green on the probability of choosing red. From the graphic, we see which is the most likely range of predicted probability for each value of no.green.
I hope it helps.