簡體   English   中英

Modelsummary:估計和統計的不同格式

[英]Modelsummary: Different formats for estimates and statistics

首先感謝modelsummary package的創建者——非常有用!

我對統計和估計的不同fmt有疑問? 這是一個代表:

    url <- 'https://vincentarelbundock.github.io/Rdatasets/csv/HistData/Guerry.csv'
dat <- read.csv(url)
models <- list(
  "OLS 1"     = lm(Donations ~ Literacy + Clergy, data = dat),
  "Poisson 1" = glm(Donations ~ Literacy + Commerce, family = poisson, data = dat),
  "OLS 2"     = lm(Crime_pers ~ Literacy + Clergy, data = dat),
  "Poisson 2" = glm(Crime_pers ~ Literacy + Commerce, family = poisson, data = dat),
  "OLS 3"     = lm(Crime_prop ~ Literacy + Clergy, data = dat)
)
modelsummary(models)
modelsummary(models, output = "flextable", estimate=glue_col("{estimate}{stars}"),
             statistic = 'statistic', stars = c('*' = .1, '**' = .05, '***'=0.01))

如何使系數下的 t 統計量具有 2 個有效數字和估計值 3? 所以第一個系數統計對將是:7948.667 (2078.27)

我仔細閱讀了文檔,但沒有找到答案。 先感謝您!

從 0.9.4 版開始,沒有直接的方法可以僅使用fmt參數或glue字符串來實現此目的。

然而,利用modelsummary網站上描述tidy_customglance_custom機制來對您的估計和統計數據進行任何后期處理是相當容易的。 這為用戶提供了自定義 output 格式的無限可能。

例如,

library(modelsummary)
library(broom)

tidy_custom.lm <- function(x) {
    out <- broom::tidy(x)
    out$statistic <- sprintf("%.2f", out$statistic)
    out$estimate <- sprintf("%.3f", out$estimate)
    return(out)
}

mod <- lm(mpg ~ hp + drat, mtcars)

modelsummary(mod, statistic = "statistic", stars = TRUE)
Model 1
(截距) 10.790*
(2.12)
生命值 -0.052***
(-5.57)
德拉特 4.698***
(3.94)
數量。 32
R2 0.741
R2 調整。 0.723
AIC 169.5
BIC 175.4
日志。 -80.752
F 41.522

注: ^^ + p < 0.1,* p < 0.05,** p < 0.01,*** p < 0.001

在 Github 上打開了一個問題。 如果您想為對話做出貢獻,請隨時在此處發帖。

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM