当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R broom tidy.robustbase.glmrob 整理 a(n) glmrob 对象


Tidy 总结了有关模型组件的信息。模型组件可能是回归中的单个项、单个假设、聚类或类。 tidy 所认为的模型组件的确切含义因模型而异,但通常是不言而喻的。如果模型具有多种不同类型的组件,您将需要指定要返回哪些组件。

用法

# S3 method for glmrob
tidy(x, conf.int = FALSE, conf.level = 0.95, ...)

参数

x

robustbase::glmrob() 返回的 glmrob 对象。

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 参数时会发出警告(如果该参数将被忽略)。

细节

对于强大模型的整理器MASS包见tidy.rlm().

也可以看看

robustbase::glmrob()

其他健壮的基础整理器:augment.glmrob()augment.lmrob()glance.lmrob()tidy.lmrob()

带有列的 tibble::tibble()

conf.high

估计置信区间的上限。

conf.low

估计置信区间的下限。

estimate

回归项的估计值。

p.value

与观察到的统计量相关的两侧 p 值。

statistic

在回归项非零的假设中使用的 T-statistic 的值。

std.error

回归项的标准误差。

term

回归项的名称。

例子


if (requireNamespace("robustbase", quietly = TRUE)) {
  # load libraries for models and data
  library(robustbase)

  data(coleman)
  set.seed(0)

  m <- lmrob(Y ~ ., data = coleman)
  tidy(m)
  augment(m)
  glance(m)

  data(carrots)

  Rfit <- glmrob(cbind(success, total - success) ~ logdose + block,
    family = binomial, data = carrots, method = "Mqle",
    control = glmrobMqle.control(tcc = 1.2)
  )

  tidy(Rfit)
  augment(Rfit)
}
#> # A tibble: 24 × 5
#>    cbind(success, total - success…¹ [,""] logdose block .fitted .resid[,1]
#>                               <int> <int>   <dbl> <fct>   <dbl>      <dbl>
#>  1                               10    25    1.52 B1     -0.726      10.7 
#>  2                               16    26    1.64 B1     -0.972      17.0 
#>  3                                8    42    1.76 B1     -1.22        9.22
#>  4                                6    36    1.88 B1     -1.46        7.46
#>  5                                9    26    2    B1     -1.71       10.7 
#>  6                                9    33    2.12 B1     -1.96       11.0 
#>  7                                1    31    2.24 B1     -2.20        3.20
#>  8                                2    26    2.36 B1     -2.45        4.45
#>  9                               17    21    1.52 B2     -0.491      17.5 
#> 10                               10    30    1.64 B2     -0.737      10.7 
#> # ℹ 14 more rows
#> # ℹ abbreviated name: ¹`cbind(success, total - success)`[,"success"]
#> # ℹ 1 more variable: .resid[2] <dbl>

相关用法


注:本文由纯净天空筛选整理自大神的英文原创作品 Tidy a(n) glmrob object。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。