我有以下数据框架:
library(tidyverse)
df <- structure(list(genes = c(
"Ccl2", "Bar", "Wnt5a", "Ccl2", "Ccl2",
"Bar", "Ccl2", "Wnt5a", "Wnt5a", "Bar", "Ccl2", "Bar", "Wnt5a",
"Bar", "Ccl2", "Ccl2", "Ccl2", "Bar", "Ccl2", "Wnt5a", "Ccl2",
"Bar", "Bar", "Wnt5a", "Bar", "Ccl2", "Ccl2", "Ccl2", "Bar",
"Bar", "Ccl2", "Wnt5a", "Ccl2", "Wnt5a", "Bar", "Bar", "Bar",
"Ccl2", "Bar", "Bar", "Wnt5a", "Wnt5a", "Wnt5a", "Wnt5a", "Ccl2",
"Wnt5a", "Bar", "Bar", "Bar", "Bar", "Ccl2"
), gexp = c(
0,
2, 0, 0, 0, 3, 0, 0, 0, 11, 0, 1, 0, 6, 0, 0, 0, 3, 0, 0, 0,
2, 0, 0, 7, 0, 0, 0, 13, 2, 0, 0, 0, 0, 4, 2, 10, 0, 2, 3, 0,
0, 0, 0, 0, 0, 6, 0, 1, 3, 11
), Treat = c(
"Non_Treat", "Non_Treat",
"Foo", "Foo", "Non_Treat", "Foo", "Saline", "Non_Treat", "Foo",
"Saline", "Foo", "Foo", "Foo", "Foo", "Non_Treat", "Foo",
"Foo", "Saline", "Non_Treat", "Foo", "Foo", "Saline", "Non_Treat",
"Foo", "Foo", "Foo", "Foo", "Non_Treat", "Foo", "Foo",
"Foo", "Foo", "Saline", "Non_Treat", "Non_Treat", "Non_Treat",
"Foo", "Foo", "Non_Treat", "Saline", "Foo", "Foo", "Non_Treat",
"Foo", "Foo", "Foo", "Non_Treat", "Non_Treat", "Non_Treat",
"Non_Treat", "Foo"
), partition = c(
1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
), cell_name = c(
"SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP",
"SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP",
"SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP",
"SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP",
"SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP",
"SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP", "SaP"
), Qux_status = c(
"Qux-",
"Qux-", "Qux-", "Qux-", "Qux-", "Qux-", "Qux-", "Qux-",
"Qux-", "Qux-", "Qux+", "Qux-", "Qux-", "Qux-", "Qux+",
"Qux-", "Qux-", "Qux-", "Qux-", "Qux-", "Qux-", "Qux-",
"Qux-", "Qux-", "Qux-", "Qux+", "Qux+", "Qux-", "Qux-",
"Qux-", "Qux-", "Qux-", "Qux-", "Qux-", "Qux-", "Qux-",
"Qux-", "Qux-", "Qux-", "Qux-", "Qux-", "Qux-", "Qux-",
"Qux-", "Qux+", "Qux-", "Qux+", "Qux-", "Qux-", "Qux-",
"Qux-"
), sname = structure(c(
4L, 4L, 6L, 6L, 4L, 6L, 5L, 4L,
6L, 5L, 3L, 6L, 6L, 6L, 1L, 6L, 6L, 5L, 4L, 6L, 6L, 5L, 4L, 6L,
6L, 3L, 3L, 4L, 6L, 6L, 6L, 6L, 5L, 4L, 4L, 4L, 6L, 6L, 4L, 5L,
6L, 6L, 4L, 6L, 3L, 6L, 1L, 4L, 4L, 4L, 6L
), .Label = c(
"Qux+.Non_Treat",
"Qux+.Saline", "Qux+.Foo", "Qux-.Non_Treat", "Qux-.Saline",
"Qux-.Foo"
), class = "factor")), row.names = c(NA, -51L), class = c(
"tbl_df",
"tbl", "data.frame"
))
我使用ggpubr::ggbarplot来制作带有错误条的barpolot。用这个代码
ggpubr::ggbarplot(df, x = "sname", y = "gexp", add = "mean_se")
我希望错误栏会出现。但它没有
它给了我一个信息:
Warning message:
Computation failed in `stat_summary()`:
object 'mean_se_' of mode 'function' was not found
我如何解决这个问题?
发布于 2021-06-14 07:55:40
我认为这种情况不应该发生,但简单的解决办法是先加载包ggpubr
library(ggpubr)
。我不知道为什么,但它在我的机器上工作。
library(ggpubr)
ggbarplot(df, x = "sname", y = "gexp", add = "mean_se")
https://stackoverflow.com/questions/67964748
复制相似问题