Whisker绘制比较簇之间的均值和方差

Question

给定具有N个变量的数据集的聚类：

data(iris)
df <- iris[,1:4]
fit <- kmeans(df, 3)

我想使用晶须图来显示结果，其中：

• x轴对应于变量或尺寸。
• 框表示聚类均值（水平线）和组的方差。
• 对于每个x值（变量），都有k个框，每个组一个。
• 还会显示数据点。

我明白了：

points <- df
points$cluster <- fit$cluster
points <- melt(points, id='cluster')

p <- ggplot(points, aes(x=variable, y=value)) 
p <- p + geom_boxplot(aes(fill = factor(cluster)))
p <- p + geom_jitter(aes(color=factor(cluster)))
print(p)

但抖动点应与相应的框对齐。我该怎么办？