Re: [R] Question on binomial data

2009-04-21 Thread ehud cohen
David, Faraway suggests using the Hosmer Lemeshow test in the case of a binary response, and discusses the inadequacy of Wald statistics. However, I'm not sure it applies here due to the limited number of cases. Thanks, Ehud. On Wed, Apr 22, 2009 at 2:04 AM, David Winsemius wrote: > Surely Fara

Re: [R] Question on binomial data

2009-04-21 Thread David Winsemius
Surely Faraway does not suggest using the Wald statistic in preference to the deviance? Even if the distribution of deviance is not exactly chi-square, it appears generally accepted that a comparison of the difference in deviance to the chi-square statistic is better than using the ratio of

Re: [R] Question on binomial data

2009-04-21 Thread ehud cohen
I thought of testing the difference in deviance between the null model and the fitted model, assuming it is distributed as chi-sq. However, Faraway writes that if the outcome is binary, the deviance distribution is far from chisq. I've done a permutation test: N<-5000; # Towards the upper limit, a

Re: [R] Question on binomial data

2009-04-21 Thread David Winsemius
You should review your course material on interpreting general linear models. The criterion you have chosen for "significance" (looking at p values for indivdiual coefficients) is not the recommended one. Seek out the section that discusses the proper method for using deviance estimates for

[R] Question on binomial data

2009-04-21 Thread ehud cohen
Hi, We have an experiment with pass/fail outcome, and a continuous parameter which may contribute to the outcome. First, we've analyzed it by: p=c(F,T,F,F,F,T,T,T,T,T,T,T,F,T,T,T,T); w=c(53,67,59,59,53,89,72,56,65,63,62,58,59,72,61,68,63); l<-glm(p~w,family=binomial) summary(l) Which turned out