我正在尝试创建此数据的 RidgePlot:
long_data <- structure(list(Cell_Type = structure(c(1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L,
6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 10L,
10L, 10L, 10L), .Label = c("CD4..T.cells", "CD4..Tcm", "CD8..Tcm",
"CD8..Tem", "Class.switched.memory.B.cells", "Monocytes", "MPP",
"naive.B.cells", "NK.cells", "Plasma.cells"), class = "factor"),
Sample = structure(c(1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L,
2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L,
1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L,
4L), .Label = c("Sample1", "Sample2", "Sample3", "Sample4"
), class = "factor"), Percentage = c(39, 35.1, 12.7405962,
4.995628, 11.4, 10.9, 8.9116355, 10.1955, 4.5, 6.1, 18.2643388,
7.005608, 5.9, 6, 4.9746425, 10.288099, 2.8, 2.8, 13.7408777,
19.657708, 18.8, 21.4, 5.7131852, 19.657708, 0.6, 0.8, 0.2308895,
18.295758, 7, 6.2, 7.210278, 3.383666, 9.9, 9.8, 16.9838566,
9.087761, 0.1, 1, 11.2297, 10.511573), outlier = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L,
1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L), .Label = c("no", "yes"
), class = "factor")), class = "data.frame", row.names = c(NA,
-40L))
列异常值是“是”和“否”两个级别的因子,但是当我创建 ggplot 时,它完全忽略了所有带有“是”的行作为异常值。
这是我得到的 ggplot 的一个例子。我似乎在这里遗漏了一些东西,因为预期的输出是一个带有“否”和“是”的图例,并且点有两种不同的颜色。
这是生成它的代码:
ggplot(long_data, aes(x=Percentage, y=reorder(Cell_Type,Percentage))) +
geom_density_ridges(scale=2, alpha=0.2)+
geom_density_ridges(alpha=0, color=NA, jittered_points=T, point_alpha=1, aes(point_color=outlier)) +
xlab("Percentage") + ylab("Cell Type")
我在这里想念什么?