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


R broom tidy.anova 整理 a(n) anova 對象


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

用法

# S3 method for anova
tidy(x, ...)

參數

x

anova 對象,例如由 stats::anova()car::Anova()car::leveneTest()car::linearHypothesis() 創建的對象。

...

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

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

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

細節

ANOVA 表的 term 列可以帶有前導或尾隨空格,此整理方法會修剪這些空格。

有關 car::leveneTest() 輸出整理的文檔,請參閱 tidy.leveneTest()

帶有列的 tibble::tibble()

df

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

meansq

平均平方和。等於平方和除以自由度。

p.value

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

statistic

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

sumsq

該術語解釋的平方和。

term

回歸項的名稱。

例子


# fit models
a <- lm(mpg ~ wt + qsec + disp, mtcars)
b <- lm(mpg ~ wt + qsec, mtcars)

mod <- anova(a, b)

# summarize model fit with tidiers
tidy(mod)
#> # A tibble: 2 × 7
#>   term                  df.residual   rss    df    sumsq statistic p.value
#>   <chr>                       <dbl> <dbl> <dbl>    <dbl>     <dbl>   <dbl>
#> 1 mpg ~ wt + qsec + di…          28  195.    NA NA       NA         NA    
#> 2 mpg ~ wt + qsec                29  195.    -1 -0.00102  0.000147   0.990
glance(mod)
#> # A tibble: 1 × 2
#>   deviance df.residual
#>      <dbl>       <dbl>
#> 1     195.          29

# car::linearHypothesis() example
library(car)
mod_lht <- linearHypothesis(a, "wt - disp")
tidy(mod_lht)
#> # A tibble: 1 × 10
#>   term   null.value estimate std.error statistic p.value df.residual   rss
#>   <chr>       <dbl>    <dbl>     <dbl>     <dbl>   <dbl>       <dbl> <dbl>
#> 1 wt - …          0    -5.03      1.23      16.6 3.39e-4          28  195.
#> # ℹ 2 more variables: df <dbl>, sumsq <dbl>
glance(mod_lht)
#> # A tibble: 1 × 2
#>   deviance df.residual
#>      <dbl>       <dbl>
#> 1     195.          28

相關用法


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