我希望在R中对基于对数似然比的glm()获得的对数线性模型执行一些逐步回归(可能是前向和后向)。
理想情况下,我正在寻找一个类似于step()的函数,它不使用AIC而是使用对数似然比。
model<-glm(Frequency~ x+s+d+i, data=table_ll, poisson(link=log))
model_sat<-glm(Frequency~ x*s*d*i, data=table_ll, poisson(link=log))
step(model, scope=list(lower=model, upper=model_sat), direction="forward")
谢谢! ç