问题描述
我正在尝试从Clause Willke的数据可视化指南(https://clauswilke.com/dataviz/visualizing-uncertainty.html#fig:cocoa-data-vs-CI)中复制图形,他编写的代码对我来说很令人窒息。
如果时间太长,我事先表示歉意,但尝试解释我的想法。 (此外,您可以在https://github.com/clauswilke/dviz.supp/blob/master/data/cacao.rda上找到cacao
数据。)
完整代码块:
cacao %>%
filter(location == "Canada") -> cacao_single
fit <- lm(rating ~ 1,data = cacao_single)
CI_df <- data.frame(type = c(0.8,0.95,0.99)) %>%
mutate(df = map(type,~tidy(emmeans(fit,~ 1,options = list(level = .x)
)))) %>%
unnest(c()) %>%
select(type,estimate,std.error,conf.low,conf.high) %>%
mutate(type = paste0(signif(100*type,2),"% confidence interval"))
Error: Can't subset columns that don't exist.
x Column `estimate` doesn't exist.
emmeans()函数通常返回beta估计值,SE,df以及置信区间的上限和下限,如:
> emmeans(fit,~1)
1 emmean SE df lower.CL upper.CL
overall 3.32 0.0379 124 3.25 3.4
Confidence level used: 0.95
当我运行上面的管道语句时,我得到了一个小标题列表,但是列表元素小标题没有附加配置项。这是select()
找不到estimate
的补充(错误消息似乎表明)。因此,
- 为什么小标题列表的
type
列中的CI与每个置信度都没有关联?应该吗? - 如何获得
select()
,然后选择?
解决方法
请尝试删除 .then((result) => {
postModules.readPostImages(result);
})
.then((result) => {
console.log(result);
// return postModules.addImagesToData(rows,result);
})
中的c()
(它会发出警告,但它可以工作)。 unnest()
语法可能已更改。
unnest()