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

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


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

用法

# S3 method for systemfit
tidy(x, conf.int = TRUE, conf.level = 0.95, ...)

參數

x

通過調用 systemfit::systemfit() 生成的 systemfit 對象。

conf.int

邏輯指示是否在整理的輸出中包含置信區間。默認為 FALSE 。

conf.level

用於置信區間的置信水平(如果 conf.int = TRUE )。必須嚴格大於 0 且小於 1。默認為 0.95，對應於 95% 的置信區間。

...

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

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

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

細節

這個整潔的方法適用於類 systemfit 的任何模型對象。默認返回六列的小標題。

也可以看看

值

帶有列的 tibble::tibble()：

conf.high

估計置信區間的上限。

conf.low

估計置信區間的下限。

estimate

回歸項的估計值。

p.value

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

std.error

回歸項的標準誤差。

term

回歸項的名稱。

例子


set.seed(27)

# load libraries for models and data
library(systemfit)
#> 
#> Please cite the 'systemfit' package as:
#> Arne Henningsen and Jeff D. Hamann (2007). systemfit: A Package for Estimating Systems of Simultaneous Equations in R. Journal of Statistical Software 23(4), 1-40. http://www.jstatsoft.org/v23/i04/.
#> 
#> If you have questions, suggestions, or comments regarding the 'systemfit' package, please use a forum or 'tracker' at systemfit's R-Forge site:
#> https://r-forge.r-project.org/projects/systemfit/

# generate data
df <- data.frame(
  X = rnorm(100),
  Y = rnorm(100),
  Z = rnorm(100),
  W = rnorm(100)
)

# fit model
fit <- systemfit(formula = list(Y ~ Z, W ~ X), data = df, method = "SUR")

# summarize model fit with tidiers
tidy(fit)
#> # A tibble: 4 × 7
#>   term            estimate std.error statistic p.value conf.low conf.high
#>   <chr>              <dbl>     <dbl>     <dbl>   <dbl>    <dbl>     <dbl>
#> 1 eq1_(Intercept)   0.109     0.0981     1.11    0.269  -0.0857    0.304 
#> 2 eq1_Z            -0.0808    0.0934    -0.865   0.389  -0.266     0.105 
#> 3 eq2_(Intercept)  -0.0495    0.110     -0.449   0.655  -0.269     0.170 
#> 4 eq2_X            -0.133     0.103     -1.30    0.198  -0.337     0.0707
tidy(fit, conf.int = TRUE)
#> # A tibble: 4 × 7
#>   term            estimate std.error statistic p.value conf.low conf.high
#>   <chr>              <dbl>     <dbl>     <dbl>   <dbl>    <dbl>     <dbl>
#> 1 eq1_(Intercept)   0.109     0.0981     1.11    0.269  -0.0857    0.304 
#> 2 eq1_Z            -0.0808    0.0934    -0.865   0.389  -0.266     0.105 
#> 3 eq2_(Intercept)  -0.0495    0.110     -0.449   0.655  -0.269     0.170 
#> 4 eq2_X            -0.133     0.103     -1.30    0.198  -0.337     0.0707

相關用法


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