r - vglm オブジェクトの段階的回帰: step4() エラー

Question

vglm オブジェクトに対して段階的な回帰を実行したいと考えています。それらはS4なので、このような方法step()は機能しません。step4()とstep4vglm()inを見つけましたVGAM。これは と同等のように見えますstep()が、S4オブジェクト用です。ただし、同じエラーが発生し続けましたError in Deviance.categorical.data.vgam(mu = mu, y = y, w = w, residuals = residuals, : arguments 'y' and 'mu' must have at least 2 columns。この問題を理解して修正するのに苦労しています。誰が何が悪いのか知っていますか？または、ステップワイズ回帰を実行するために変換S4してS3使用する方法はありますか? step()役立つ場合は、reprex を添付しました。

##libraries
library(VGAM)
#> Loading required package: stats4
#> Loading required package: splines

# import data
dff<-list()
dff[[1]] <- data.table::data.table(
                   sbq_1 = c(1L,2L,1L,2L,2L,
                             1L,2L,2L,4L,2L,1L,4L,2L,3L,1L,4L,2L,1L,
                             3L,1L,4L,1L,3L,4L,2L,1L,2L,2L,2L,2L,
                             1L,4L,3L,4L,1L,2L,3L,4L,1L,1L,3L,3L,1L,
                             1L,2L,1L,1L,1L,1L),
                     age = c(53L,32L,47L,58L,
                             34L,29L,57L,25L,38L,29L,51L,50L,37L,33L,
                             41L,36L,49L,38L,32L,51L,37L,36L,37L,39L,
                             55L,55L,45L,25L,24L,27L,26L,57L,47L,30L,
                             37L,39L,38L,34L,38L,56L,39L,43L,49L,39L,48L,
                             44L,63L,61L,46L),
                  gender = c(1L,1L,2L,1L,1L,
                             1L,1L,2L,2L,1L,1L,1L,1L,1L,1L,2L,1L,1L,
                             2L,1L,1L,1L,1L,1L,2L,1L,1L,2L,1L,1L,
                             1L,1L,2L,1L,2L,1L,1L,2L,1L,1L,1L,2L,1L,
                             1L,2L,1L,1L,2L,1L),
               ethnicity = c(2L,2L,2L,2L,2L,
                             2L,2L,2L,2L,2L,2L,1L,1L,2L,2L,2L,2L,2L,
                             2L,2L,2L,2L,2L,2L,2L,2L,2L,2L,2L,2L,
                             2L,2L,2L,2L,2L,2L,2L,2L,2L,2L,2L,2L,2L,
                             2L,2L,2L,2L,2L,2L)
            )

ppom.1 <- vglm(sbq_1 ~ ., family=cumulative(parallel=F~ age), maxit=50, data=dff[[1]])
ppom.2 <- step4(ppom.1)
#> Start:  AIC=135.61
#> sbq_1 ~ age + gender + ethnicity
#> Error in Deviance.categorical.data.vgam(mu = mu, y = y, w = w, residuals = residuals, : arguments 'y' and 'mu' must have at least 2 columns
summary(ppom.2)
#> Error in h(simpleError(msg, call)): error in evaluating the argument 'object' in selecting a method for function 'summary': object 'ppom.2' not found

^{reprex パッケージ(v2.0.1)により 2021-10-26 に作成}