0

我想知道如何通过 or 将多重插补摘要的结果导出到Ameliatex 或 html 。但是，两种类型的包都显示错误。texregstargazer

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 的多重插补结果？

附言

> 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。

4

2 回答 2

2

目前并非所有模型对象都支持在这两个包中，我猜您正在使用的模型就是这种情况。如果 stargazer 不支持您的模型对象，我不知道有解决方案。但是，texreg 是可扩展的：您可以编写一个简单的函数来从模型对象中提取您需要的内容，然后普通的 texreg 命令将与该提取函数一起使用。这在这篇关于 texreg 的文章中有清楚的解释，见第六节。

texreg 的作者还写了一个如何编写提取函数的示例来回答我在这里提出的问题，并在此处写了一篇更全面的关于如何扩展 texreg 的评论。我建议您编写与您的模型对象一起使用的提取函数，然后您可以使用 texreg。

于 2016-10-21T23:23:26.757 回答
0

texreg 现在处理 Zelig 对象并通过 Zelig 与 Amelia、Mice 和 Mi 一起工作。使用函数 to_zelig_mi 创建 Zelig 可以识别的数据对象。

于 2018-10-08T14:25:39.860 回答