Question

我尝试使用stargazer包创建回归结果表。不幸的是stargazer由于奇点而拒绝创建表格。由于这些奇点来自多个层次的重叠因素，因此不容易避免。这就是我寻找替代解决方案的原因这是一个有效的例子 adapted from a related thread：

library(sandwich)
library(stargazer)
set.seed(101)
dat<-data.frame(one=c(sample(1000:1239)),
            two=c(sample(200:439)),
            three=c(sample(600:839)),
            Jan=c(rep(1,20),rep(0,220)),
            Feb=c(rep(0,20),rep(1,20),rep(0,200)),
            Mar=c(rep(0,40),rep(1,20),rep(0,180)),
            Apr=c(rep(0,60),rep(1,20),rep(0,160)),
            May=c(rep(0,80),rep(1,20),rep(0,140)),
            Jun=c(rep(0,100),rep(1,20),rep(0,120)),
            Jul=c(rep(0,120),rep(1,20),rep(0,100)),
            Aug=c(rep(0,140),rep(1,20),rep(0,80)),
            Sep=c(rep(0,160),rep(1,20),rep(0,60)),
            Oct=c(rep(0,180),rep(1,20),rep(0,40)),
            Nov=c(rep(0,200),rep(1,20),rep(0,20)),
            Dec=c(rep(0,220),rep(1,20))) 
model <- lm(one ~ two + three + Jan + Feb + Mar + Apr + May + Jun + Jul + Aug + Sep + Oct + Nov +    Dec, data=dat)
summary(model)
stargazer(model)

由于奇点，我得到错误：下标越界。你对如何解决这个问题有任何想法吗？

Answer 1

此问题已得到修复。最新版本的stargazer（来自3.0）应该对你的例子没有任何问题。

stargazer（）用于具有奇点的模型

1 个答案: