當前位置: 首頁>>代碼示例 >>用法及示例精選 >>正文

R broom tidy.epi.2by2 整理 a(n) epi.2by2 對象


Tidy 總結了有關模型組件的信息。模型組件可能是回歸中的單個項、單個假設、聚類或類。 tidy 所認為的模型組件的確切含義因模型而異，但通常是不言而喻的。如果模型具有多種不同類型的組件，您將需要指定要返回哪些組件。

用法

# S3 method for epi.2by2
tidy(x, parameters = c("moa", "stat"), ...)

參數

x

通過調用 epiR::epi.2by2() 生成的 epi.2by2 對象

parameters

返回關聯度量 ( moa ) 或測試統計量 ( stat )，默認為 moa (關聯度量)

...

附加參數。不曾用過。僅需要匹配通用簽名。注意：拚寫錯誤的參數將被吸收到 ... 中，並被忽略。如果拚寫錯誤的參數有默認值，則將使用默認值。例如，如果您傳遞 conf.lvel = 0.9 ，所有計算將使用 conf.level = 0.95 進行。這裏有兩個異常：

  • tidy() 方法在提供 exponentiate 參數時會發出警告(如果該參數將被忽略)。

  • augment() 方法在提供 newdata 參數時會發出警告(如果該參數將被忽略)。

細節

當調用 epiR::epi.2by2() 時，tibble 對於 massoc 或 massoc.detail 中包含的每個關聯度量或測試都有一列。

也可以看看

值

帶有列的 tibble::tibble()：

conf.high

估計置信區間的上限。

conf.low

估計置信區間的下限。

df

該術語在模型中使用的自由度。

p.value

與觀察到的統計量相關的兩側 p 值。

statistic

在回歸項非零的假設中使用的 T-statistic 的值。

term

回歸項的名稱。

estimate

估計關聯度

例子


# load libraries for models and data
library(epiR)
#> Package epiR 2.0.61 is loaded
#> Type help(epi.about) for summary information
#> Type browseVignettes(package = 'epiR') to learn how to use epiR for applied epidemiological analyses
#> 

# generate data
dat <- matrix(c(13, 2163, 5, 3349), nrow = 2, byrow = TRUE)

rownames(dat) <- c("DF+", "DF-")
colnames(dat) <- c("FUS+", "FUS-")

# fit model
fit <- epi.2by2(
  dat = as.table(dat), method = "cross.sectional",
  conf.level = 0.95, units = 100, outcome = "as.columns"
)

# summarize model fit with tidiers
tidy(fit, parameters = "moa")
#> # A tibble: 15 × 4
#>    term                estimate conf.low conf.high
#>    <chr>                  <dbl>    <dbl>     <dbl>
#>  1 PR.strata.wald         4.01    1.43      11.2  
#>  2 PR.strata.taylor       4.01    1.43      11.2  
#>  3 PR.strata.score        1.84    1.25       2.24 
#>  4 OR.strata.wald         4.03    1.43      11.3  
#>  5 OR.strata.cfield       4.03   NA         NA    
#>  6 OR.strata.score        4.03    1.49      10.9  
#>  7 OR.strata.mle          4.02    1.34      14.4  
#>  8 ARisk.strata.wald      0.448   0.0992     0.797
#>  9 ARisk.strata.score     0.448   0.142      0.882
#> 10 NNT.strata.wald      223.    125.      1008.   
#> 11 NNT.strata.score     223.    113.       705.   
#> 12 PARisk.strata.wald     0.176  -0.0225     0.375
#> 13 PARisk.strata.piri     0.176   0.0389     0.314
#> 14 AFRisk.strata.wald     0.750   0.301      0.911
#> 15 PAFRisk.strata.wald    0.542   0.0361     0.782
tidy(fit, parameters = "stat")
#> # A tibble: 3 × 4
#>   term               statistic    df p.value
#>   <chr>                  <dbl> <dbl>   <dbl>
#> 1 chi2.strata.uncor       8.18     1 0.00424
#> 2 chi2.strata.yates       6.85     1 0.00885
#> 3 chi2.strata.fisher     NA       NA 0.00635
源代碼：R/epiR-tidiers.R

相關用法


注：本文由純淨天空篩選整理自等大神的英文原創作品 Tidy a(n) epi.2by2 object。非經特殊聲明，原始代碼版權歸原作者所有，本譯文未經允許或授權，請勿轉載或複製。