我尝试使用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)由于奇异点,我得到了错误:下标越界。你有什么办法解决这个问题吗?
发布于 2013-01-31 03:28:35
此问题已修复。对于您的示例,较新版本的stargazer (从3.0开始)应该没有问题。
https://stackoverflow.com/questions/13593433
复制相似问题