Показать результат Амелии, множественного вменения, texreg или stargazer

Question

Я хотел бы знать, как экспортировать результаты множественного вменения Amelia Сводка в текс или HTML по texreg или же stargazer, Однако оба типа пакетов показывают ошибку.

library(Amelia)
library(Zelig)
library(texreg)
library(stargazer)

mi <-amelia(x = d,m = 5)
summary(mi)
lmi <- zelig(y ~ x1+x2,
 data = m, model = "logit")
summary(lmi)
mi2 <-amelia(x = d,m = 5)
summary(mi2)
lmi2 <- zelig(y ~ x1+x2,
 data = m2, model = "logit")
summary(lmi2)

stargazer(lmi,lmi2, title="hogehoge", omit.stat=c("f","ser"), align=T, no.space=T)

>Error in envRefInferField(x, what, getClass(class(x)), selfEnv) :
>‘result’ is not a valid field or method name for reference class “Zelig-logit”

models <- list(lmi,lmi2)
texreg(l=models,file="texreg.tex", caption="texreg",
    digits=3, booktabs=T, dcolumn=T, center=T, use.packages=F, caption.above=T, custom.model.names=c("model1", "model2"))

>Error in (function (classes, fdef, mtable)  :
>unable to find an inherited method for function ‘extract’ for signature ‘"Zelig-logit"’

Как я могу показать результат множественного вменения R?

PS

> summary(lmi)
Model: Combined Imputations 
                    Estimate Std.Error  z value  Pr(>|z|)    
(Intercept)         -4.15186   0.41048 -10.1147 0.000e+00 ***
x1                   9.76225   0.18356  12.4336 0.000e+00 ***
x2                   1.76719   0.18031  10.9489 0.000e+00 ***
---
Signif. codes:  '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

For results from individual imputed datasets, use summary(x, subset = i:j)
Next step: Use 'setx' method

lmi2 похож на lmi.

Answer 1

Не все объекты модели в настоящее время поддерживаются ни в одном пакете, и я полагаю, что это относится к моделям, которые вы используете. Если Stargazer не поддерживает ваш модельный объект, я не знаю решения. Тем не менее, texreg расширяемый: вы можете написать простую функцию для извлечения того, что вам нужно из объекта модели, и тогда обычные команды texreg будут работать с этой функцией извлечения. Это ясно объясняется в этой статье на texreg, см. Раздел шесть.

Автор texreg также написал пример того, как написать функцию извлечения в ответ на вопрос, который я задал здесь, и написал еще более полный обзор того, как расширить texreg здесь. Я бы посоветовал вам написать функции извлечения, которые работают с объектами вашей модели, а затем вы можете использовать texreg.

Answer 2

texreg теперь обрабатывает объекты Zelig и работает с Amelia, Mice и Mi через Zelig. Используйте функцию to_zelig_mi, чтобы создать объект данных, который распознает Zelig.