【问题标题】:Computation failed in `stat_signif()`: data are essentially constant`stat_signif()` 中的计算失败:数据本质上是恒定的
【发布时间】:2021-03-29 10:36:33
【问题描述】:

我正在尝试制作箱线图,我想查看组比较,但是由于某种原因它没有显示

library(ggbeeswarm)
library(ggpubr)
library(ggplot2)
ggplot(data = data2, aes(wcc_group, 2750.3)) +
geom_boxplot() +
geom_quasirandom(aes(colour = as.factor(X28_day_mortality) , shape = sex, fill = 
as.factor(X28_day_mortality)), groupOnX = FALSE) +
labs(x = "WCC Group", y = "Apolipoprotein A-I", title = "Expression of Apolipoprotein A-I in WCC 
Groups", fill = "Mortality", shape = "Sex", colour = "Mortality") +  stat_compare_means(method = 
"t.test", comparisons = list(c("Low", "Medium"), c("Low", "High"), c("Medium", "High") ), aes(label = 
..p.signif..)) +
 scale_colour_manual(values =c("red", "black")) +
 scale_fill_manual(values =c("red", "black")) +
 scale_shape_manual(values =c(24,22))

enter image description here

这个错误一直出现的原因是什么?

【问题讨论】:

  • 我无法重现您遇到的错误,因为找不到您在代码中使用的data2。您能否提供data2 或提供一个代表您的问题的最小的、可重现的示例?

标签: r ggplot2 comparison boxplot p-value


【解决方案1】:

当您在示例中指定 y 值时,您将其设为 2750.3,这是有原因的吗?

这个奇异的 y 值是为每个组绘制黑线的原因。它不是箱线图,因为只有一个值 - 不能有 Q1、中位数、Q3 等。如果您将 y 值作为该变量的值(载脂蛋白 AI?),那么箱线图将是能够显示分数的分布。

t 检验不起作用,并且由于类似原因引发错误 - 计算 t 检验值没有方差,因为 y 的唯一值是 2750.3。

ggplot(data = data2, aes(wcc_group, [enter associated protein var here])) +
geom_boxplot() +
geom_quasirandom(aes(colour = as.factor(X28_day_mortality) , shape = sex, fill = 
as.factor(X28_day_mortality)), groupOnX = FALSE) +
labs(x = "WCC Group", y = "Apolipoprotein A-I", title = "Expression of Apolipoprotein A-I in WCC 
Groups", fill = "Mortality", shape = "Sex", colour = "Mortality") +  stat_compare_means(method = 
"t.test", comparisons = list(c("Low", "Medium"), c("Low", "High"), c("Medium", "High") ), aes(label = 
..p.signif..)) +
 scale_colour_manual(values =c("red", "black")) +
 scale_fill_manual(values =c("red", "black")) +
 scale_shape_manual(values =c(24,22))

【讨论】:

    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2022-08-16
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多