2012-11-27 80 views
1

我尝试使用stargazer程序包创建回归结果表。 不幸的是stargazer拒绝创建表,因为奇点。由于这些奇点来自多层次的重叠因素,因此不容易避免。 That's为什么我要寻找一个替代解决方案
这里的工作示例 adapted from a related threadstargazer()适用于奇点型号

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) 

因为奇异性的我得到的错误:下标越界。你有什么想法如何解决这个问题?

回答

1

此问题已得到解决。你的例子中最新版本的stargazer(从3.0开始)应该没有问题。