Таблица вывода итоговых данных латекса или HTML для объектов регрессии vglm (VGAM)

Question

Я пытаюсь получить латексный или html-вывод результатов регрессии модели VGAM (в приведенном ниже примере это обобщенный порядковый логит). Но пакеты, которые я знаю для этой цели, не работают с объектом vglm.

Здесь вы можете увидеть небольшой пример игрушки с сообщениями об ошибках, которые я получаю:

library(VGAM)
n <- 1000
x <- rnorm(n)
y <- ordered( rbinom(n, 3, prob=.5) )

ologit <- vglm(y ~ x,
            family =  cumulative(parallel = F , reverse = TRUE), 
            model=T)

library(stargazer)
stargazer(ologit)

Ошибка в объектах [[i]]$zelig.call: оператор $ не определен для этого класса S4

library(texreg)
htmlreg(ologit)

Ошибка в (функция (классы, fdef, mtable): невозможно найти унаследованный метод для функции "extract" для подписи "vglm" '

library(memisc)
mtable(ologit)

Ошибка в UseMethod("getSummary"): нет применимого метода для getSummary, примененного к объекту класса "c('vglm', 'vlm', 'vlmsmall')"

Answer 1

У меня просто была такая же проблема. Моя первая работа - запустить регрессию OLogit с polr функция MASS пакет. Получающиеся объекты легко визуализируются / суммируются обычными пакетами (я рекомендую sjplot "s tab_model функция для вывода таблицы!)

Второй вариант - создать свою собственную таблицу, которую вы затем превращаете в аккуратный HTML-объект с помощью stargazer,

Для этого вам нужно знать, что объекты s4 не могут быть поднаборными так же, как обычные объекты ( http://adv-r.had.co.nz/Subsetting.html). Наиболее прямолинейным решением является поднабор объекта, то есть извлечение соответствующих аспектов с помощью @ вместо $ условное обозначение:

sumobject <- summaryvglm(yourvglmobject)
stargazer(sumpbject@coef3, type="html", out = "RegDoc.doc")

Немного громоздко, но это помогло мне. Надеюсь это поможет!