微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

在R中使用stargazer软件包以使用引导程序输出分位数回归

如何解决在R中使用stargazer软件包以使用引导程序输出分位数回归

生成带有自举错误的分位数回归的乳胶表时,我遇到了一些问题。我可以运行回归并使用summary()获得自举错误。但是,当我尝试使用stargazer()进行操作时,要么永久占用了数据集(因为我仍然无法查看其输出),要么打印了一些与summary()不匹配的奇怪输出当我使用RStudio的样本数据集时。

library(quantreg)
library(stargazer)


dtset <- mtcars
nreps=20
clusvar=dtset$cyl

qreg = rq(mpg ~ wt + disp,tau = 0.10,data = dtset)

set.seed(2020)
summary(qreg,se="boot",R=nreps,cluster=clusvar)

set.seed(2020)
stargazer(qreg,type = "text",rq.se="boot",cluster=clusvar)

摘要显示

Call: rq(formula = mpg ~ wt + disp,tau = 0.1,data = dtset)

tau: [1] 0.1

Coefficients:
            Value    Std. Error t value  Pr(>|t|)
(Intercept) 30.85455 26.60892    1.15956  0.25569
wt          -2.90311  8.41730   -0.34490  0.73266
disp        -0.01830  0.15225   -0.12023  0.90513

对于stargazer(),它会打印:

========================================
                 Dependent variable:    
             ---------------------------
                         mpg            
----------------------------------------
wt                    -2.903**          
                       (1.291)          
                                        
disp                  -0.018***         
                       (0.006)          
                                        
Constant              30.855***         
                       (2.939)          
                                        
----------------------------------------
Observations             32             
========================================
Note:        *p<0.1; **p<0.05; ***p<0.01

==
20
--

===============================================================
6 6 4 6 8 6 8 4 4 6 6 8 8 8 8 8 8 4 4 4 4 8 8 8 8 4 4 4 8 6 8 4
---------------------------------------------------------------

关于可能发生的事情以及如何从自举分位数回归中获得一些格式良好的表的任何想法吗?

我需要运行许多回归分析,并且非常希望使用一种自动化的方法

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。