【发布时间】:2020-08-23 12:17:08
【问题描述】:
以下代码来自 help("tidy.speedglm")。它没有得到正确的 p 值。 broom 的先前版本运行良好,但 broom (0.7.0) 的当前版本无法运行。我想知道这是一个错误还是什么?谢谢。
library(speedglm)
library(broom)
clotting <- data.frame(
u = c(5, 10, 15, 20, 30, 40, 60, 80, 100),
lot1 = c(118, 58, 42, 35, 27, 25, 21, 19, 18)
)
fit <- speedglm(lot1 ~ log(u), data = clotting, family = Gamma(log))
tidy(fit)
#> # A tibble: 2 x 5
#> term estimate std.error statistic p.value
#> <chr> <dbl> <dbl> <dbl> <dbl>
#> 1 (Intercept) 5.50 0.190 28.9 NA
#> 2 log(u) -0.602 0.0553 -10.9 NA
由reprex package (v0.3.0) 于 2020 年 8 月 23 日创建
【问题讨论】:
-
看起来这是
broom和speedglm的当前迭代不能很好地协同工作,显然是因为speedglm将p 值存储为字符而不是数字。一些讨论可以在github.com/tidymodels/broom/pull/881 和github.com/MarcoEnea/speedglm/issues/2 上找到) -
@duckmayr 谢谢。我的包依赖于
broom和speedglm。我会等待speedglm的cran更新。