{
        attach.targets()
	kyph.step$null.deviance <- 83.2344746889
	T
}
all.equal(kyph.step, step(kyph.fit1,  ~ . + Age^2 + Age:Start, trace = F))
# families
all.equal(binomial(), bin.fam)
#offsets
testfit(kyph.off1)