在R中,我想创建一个也显示所有数据点的箱形图。您可以在许多帖子和网站上找到此信息,但它们似乎都以“抖动”或“随机”样式显示数据点。这是将ToothGrowth数据集与R中的ggplot2结合使用的示例代码。
library(datasets)
data(ToothGrowth)
ToothGrowth$dose <- as.factor(ToothGrowth$dose)
library(ggplot2)
ggplot(ToothGrowth, aes(x=dose, y=len)) +
geom_boxplot(notch = TRUE) +
geom_jitter(position=position_jitter(0.2))
但是,我希望数据点从左下角的最低点到右上角的最高点排序。请参见此链接中的示例: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3966983/figure/F1/(可免费访问)。具体来说,我指的是图1a的顶部(“纯度”)。
有人可能有建议吗?我将不胜感激。谢谢。
答案 0 :(得分:2)
我不知道这是您要追求的,但是也许您可以从下面的代码中激发自己的灵感。
set.seed(1234)
n <- 20
x <- rnorm(n)
boxplot(x)
points(seq(0.75, 1.25, length.out = n), sort(x))
df1 <- sapply(1:4, function(i) rnorm(n, mean = i))
df1 <- as.data.frame(df1)
df1 <- reshape2::melt(df1)
boxplot(value ~ variable, df1)
sp <- split(df1, df1$variable)
for(i in 1:4){
points(seq(i - 0.25, i + 0.25, length.out = n), sort(sp[[i]]$value))
}
编辑。
ggplot2
解决方案使用类似的技巧来定义点的x轴坐标。 “ strange” 唯一的事情是依靠R的因子从1开始的连续整数的内部表示形式。请注意,这必须被视为一种技巧,但作为可靠的方法,我不认为相信它将永远改变。
library(ggplot2)
library(tidyverse)
df1 %>%
group_by(variable) %>%
arrange(value) %>%
mutate(xcoord = seq(-0.25, 0.25, length.out = n())) %>%
ggplot(aes(x = variable, y = value, group = variable)) +
geom_boxplot() +
geom_point(aes(x = xcoord + as.integer(variable)))