我正在尝试从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()
,然后选择?答案 0 :(得分:0)
请尝试删除 .then((result) => {
postModules.readPostImages(result);
})
.then((result) => {
console.log(result);
// return postModules.addImagesToData(rows, result);
})
中的c()
(它会发出警告,但它可以工作)。 unnest()
语法可能已更改。
unnest()