Question

希望有人能指出正确的文档或提示此查询。

鉴于以下内容（dist =“weibull”是默认值）....

library(survival)
data(lung)
ff <- survreg(Surv(time, status) ~ ph.ecog + age + sex, lung)
summary(ff)

汇总调用中的系数如下所示列出，其中包含（幸存）'scale'参数作为输出的一部分......

Call:
survreg(formula = Surv(time, status) ~ ph.ecog + age + sex, data = lung)
           Value Std.     Error     z        p
(Intercept)  6.27344    0.45358 13.83 1.66e-43
ph.ecog     -0.33964    0.08348 -4.07 4.73e-05
age         -0.00748    0.00676 -1.11 2.69e-01
sex          0.40109    0.12373  3.24 1.19e-03
Log(scale)  -0.31319    0.06135 -5.11 3.30e-07

Scale= 0.731  

Weibull distribution
<...snip...>

然后如果我运行以下内容......

library(stargazer)
stargazer(ff, type="text", keep.stat="all")

我得到（非常整洁！）的观星者输出，但没有比例参数...

====================================
             Dependent variable:    
         ---------------------------
                    time
------------------------------------
ph.ecog           -0.340***         
                   (0.083)          

age                -0.007           
                   (0.007)          

sex               0.401***          
                   (0.124)          

Constant          6.273***          
                   (0.454)             

====================================
====================================
Note:    *p<0.1; **p<0.05; ***p<0.01

stargazer pdf（https://cran.r-project.org/web/packages/stargazer/stargazer.pdf）第22页建议支持scale参数，所以我可能误解了意图，或者在某个地方跳过了页面。

Answer 1

与包裹作者保持联系（谢谢，马雷克！）并得到了使用add.lines的提示，以下似乎是一个良好的开端。最后一个空字符串是强制换行，在我挖掘ff对象时，我提取了分配类型，以便在列标题中使用。

stargazer(ff, 
          type ="text", 
          keep.stat="all",
          column.labels = ff$dist,
          add.lines=list(c("Log(scale)", signif(log(ff$scale),3)),
                         c("Scale", signif(ff$scale,3)),
                         c("")
                         )
          )

这会将以下内容添加到上面提到的默认系数中......

--------------------------------------
Log(scale)           -0.313           
Scale                 0.731

将Survreg'scale'参数提取到Stargazer输出表

1 个答案: