1

表 data1 有大约 35 万个观测值。我想估计下面的 6 个模型并使用 stargazer 将结果输出到乳胶中。我应该指定 y1 是一个二进制变量，我正在处理 100 家公司。恐怕我没有任何数据可以发布。

这是我的代码。问题是每个估计都存储在 RAM 中，并且没有足够的内存来运行 Stargazer。

我有两个问题？

  1. 有没有办法将“glm”对象存储在磁盘上，然后用 stargazer 调用它们？

  2. stargazer 是否需要整个“glm”对象来输出乳胶代码。？

    l0 <- glm(y1~ x1 + log(x2)+ x3+ factor(x5)+ factor(firms) ,data=data1, family=binomial(link=logit), model=FALSE)

    l1 <- glm(y1~ x1 + log(x2)+ x3+
            factor(x5)+  x4+factor(firms) ,data = data1,family=binomial(link=logit), model=FALSE)

l2 <- glm(y1~x1 +log(x2)+  x3+ 
            factor(x5)+x4+ factor(x6) + factor(firms), data=data1,family=binomial(link=logit)
          , model=FALSE)

l3 <- glm(y1~x1 +log(x2) + x3+
            factor(x5) + x4+factor(x6)   + x7+ factor(firms)
         ,data = data1,family=binomial(link=logit), model=FALSE)

l4 <- glm(y1~x1 + log(x2) + x3+ 
            factor(x5) +x4+ factor(x6)+ x7+installments + 
            factor(firms) ,data = data1,family=binomial(link=logit), model=FALSE)

l5 <- glm(y1~x1 + log(x2) + x3+ 
            factor(x5) +x4+ factor(x6)+ x7 +installments + x8+
            factor(firms) ,data = data1,family=binomial(link=logit), model=FALSE)


stargazer(l0,l1,l2,l3,l4,title="Regression Results with Fixed Effects",     align=TRUE,apply.coef=or
      ,out = "path.tex", covariate.labels=covlabel,omit="firms",
      omit.labels="Firms", omit.yes.no=c("Yes","No"))
4

2 回答 2

0

根据stargazer manual，它不支持您在撰写本文时提到的软件包。至于其他选项，我会先尝试使用 stock 进行分析lm()。如果它不会导致与 RAM 相关的问题，您将可以在估算速度和格式化表格的难易程度之间做出选择。

您也可以选择联系开发人员，要求将其添加bigglm到speedglm功能愿望清单中。

于 2014-03-22T12:54:06.593 回答
0

面对这个问题时，我的方法是首先使用lmtest包将 *lm 对象转换为“coeftest”类。有关更多详细信息，请参阅我在此处对相关问题的回答。

于 2016-10-13T19:41:32.637 回答