emmeans (R) 拦截功能是否损坏?
Is the emmeans (R) intercept-only function broken?
我注意到 emmeans(在 R 中)在最新更新后不适用于仅拦截估计。
可重现的例子:
test=lm(mpg~1,mtcars)
library(emmeans)
emmeans::emmeans(test,~1)
我的两台机器(windows 和 Linux)的输出是:
> emmeans::emmeans(test,~1)
Error in `[[<-.data.frame`(`*tmp*`, ".wgt.", value = 2) :
replacement has 1 row, data has 0
这是一个已知问题,还是我搞砸了我的系统?
我相信这曾经有效。
如果包含变量, 是否有效:
test2=lm(mpg~as.factor(cyl),mtcars)
emmeans(test2,~cyl)
非常感谢您的提前帮助。
emmeans - 1.4.6
在 macOS Catalina 10.15.4
和 R 4.0
上工作正常
emmeans::emmeans(test,~1)
# 1 emmean SE df lower.CL upper.CL
# overall 20.1 1.07 31 17.9 22.3
#Confidence level used: 0.95
事实证明,我们在这里看到的问题 #197 -- and incorporated in CRAN version 1.47 -- created the issue (#206) 的修复。我想我现在把它们都修好了:
require(emmeans)
## Loading required package: emmeans
#206...
warp.lm <- lm(breaks ~ wool * tension, data = warpbreaks)
emmeans(warp.lm, "1")
## 1 emmean SE df lower.CL upper.CL
## overall 28.1 1.49 48 25.2 31.1
##
## Results are averaged over the levels of: wool, tension
## Confidence level used: 0.95
emmeans(warp.lm, "1", by = "wool")
## wool = A:
## 1 emmean SE df lower.CL upper.CL
## overall 31.0 2.11 48 26.8 35.3
##
## wool = B:
## 1 emmean SE df lower.CL upper.CL
## overall 25.3 2.11 48 21.0 29.5
##
## Results are averaged over the levels of: tension
## Confidence level used: 0.95
#197...
model <- lm(Sepal.Length ~ poly(Petal.Length,2), data = iris)
emtrends(model, ~ 1, "Petal.Length", max.degree = 2)
## degree = linear:
## 1 Petal.Length.trend SE df lower.CL upper.CL
## overall 0.4474 0.0180 147 0.4119 0.483
##
## degree = quadratic:
## 1 Petal.Length.trend SE df lower.CL upper.CL
## overall 0.0815 0.0132 147 0.0554 0.108
##
## Confidence level used: 0.95
由 reprex package (v0.3.0)
于 2020-06-01 创建
现在需要这个的用户可以通过
从github安装
remotes::install_github("rvlenth/emmeans")
我注意到 emmeans(在 R 中)在最新更新后不适用于仅拦截估计。
可重现的例子:
test=lm(mpg~1,mtcars)
library(emmeans)
emmeans::emmeans(test,~1)
我的两台机器(windows 和 Linux)的输出是:
> emmeans::emmeans(test,~1)
Error in `[[<-.data.frame`(`*tmp*`, ".wgt.", value = 2) :
replacement has 1 row, data has 0
这是一个已知问题,还是我搞砸了我的系统? 我相信这曾经有效。
如果包含变量, 是否有效:
test2=lm(mpg~as.factor(cyl),mtcars)
emmeans(test2,~cyl)
非常感谢您的提前帮助。
emmeans - 1.4.6
在 macOS Catalina 10.15.4
和 R 4.0
emmeans::emmeans(test,~1)
# 1 emmean SE df lower.CL upper.CL
# overall 20.1 1.07 31 17.9 22.3
#Confidence level used: 0.95
事实证明,我们在这里看到的问题 #197 -- and incorporated in CRAN version 1.47 -- created the issue (#206) 的修复。我想我现在把它们都修好了:
require(emmeans)
## Loading required package: emmeans
#206...
warp.lm <- lm(breaks ~ wool * tension, data = warpbreaks)
emmeans(warp.lm, "1")
## 1 emmean SE df lower.CL upper.CL
## overall 28.1 1.49 48 25.2 31.1
##
## Results are averaged over the levels of: wool, tension
## Confidence level used: 0.95
emmeans(warp.lm, "1", by = "wool")
## wool = A:
## 1 emmean SE df lower.CL upper.CL
## overall 31.0 2.11 48 26.8 35.3
##
## wool = B:
## 1 emmean SE df lower.CL upper.CL
## overall 25.3 2.11 48 21.0 29.5
##
## Results are averaged over the levels of: tension
## Confidence level used: 0.95
#197...
model <- lm(Sepal.Length ~ poly(Petal.Length,2), data = iris)
emtrends(model, ~ 1, "Petal.Length", max.degree = 2)
## degree = linear:
## 1 Petal.Length.trend SE df lower.CL upper.CL
## overall 0.4474 0.0180 147 0.4119 0.483
##
## degree = quadratic:
## 1 Petal.Length.trend SE df lower.CL upper.CL
## overall 0.0815 0.0132 147 0.0554 0.108
##
## Confidence level used: 0.95
由 reprex package (v0.3.0)
于 2020-06-01 创建现在需要这个的用户可以通过
从github安装remotes::install_github("rvlenth/emmeans")