[R] Quasi-binomial GLM and model selection

[Heather Kharouba]

Hi,

I'm using a GLM with a quasi binomial error distribution and I would like
to do a model selection method similar to step(AIC) to carry out a
restricted search for the "best" model. I would like to know which of my 5
predictor variables would be included in the "best" model if I start with
a 'full' model (fullbinom in this case). However, AIC can't be calculated
on a quasi-binomial model. Is there another similar method I can use?

Here's my R code so far:

(AUC ranges from 0 to 1)

fail<-1-AUC;  succfail<-cbind(AUC, fail)
fullbinom<-glm(succfail~variables+log(area)+resolution+model+taxa,family=quasibinomial(link="logit"),data=auc);

Thanks in advance,
Heather

-- 
Heather Kharouba
PhD candidate
University of British Columbia