data(hcrabs)
attach(hcrabs)
y <- ifelse(num.satellites>0,1,0)
bnfit <- glm(y~color+spine+width+weight,family=binomial)
rsq.partial(bnfit)
bnfitr <- glm(y~color+weight,family=binomial)
rsq.partial(bnfit,bnfitr)
quasibn <- glm(y~color+spine+width+weight,family=quasibinomial)
rsq.partial(quasibn)
quasibnr <- glm(y~color+weight,family=binomial)
rsq.partial(quasibn,quasibnr)
Run the code above in your browser using DataLab