r - 尤度比に基づく対数線形モデルの順方向ステップワイズ回帰

Question

LogLikelihood 比に基づいて、glm() を使用して取得した対数線形モデルに対して、R で段階的な回帰 (おそらく前方と後方の両方) を実行したいと考えています。

理想的には、AIC ではなく対数尤度比を使用する step() に似た関数を探しています。

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")

ありがとう！c