使用此代码,我可以理解我可以绘制各个聚类。
library(fpc)
set.seed(665544)
n <- 600
x <- cbind(runif(10, 0, 10)+rnorm(n, sd=0.2), runif(10, 0, 10)+rnorm(n,sd=0.2))
ds <- dbscan(x, 0.2)
plot(x[ds$cluster %in% 1:3,])
对于上面的代码,我得到以下输出,
ds
dbscan Pts=600 MinPts=5 eps=0.2
0 1 2 3 4 5 6 7 8 9 10 11
border 28 4 4 8 5 3 3 4 3 4 6 4
seed 0 50 53 51 52 51 54 54 54 53 51 1
total 28 54 57 59 57 54 57 58 57 57 57 5
我得到了所有11个集群。
如果我想在特定群集中显示数据,我需要为此编写什么R脚本。
答案 0 :(得分:5)
为了完整性:
正如您自己发现的那样,您可以使用
绘制单个群集plot(ds$cluster = yourclusternum)
但是,您可以通过绘制所有群集,但使用不同的颜色来做得更好: 如果您还想绘制噪点,请执行
plot(x, col=ds$cluster)
如果你想隐藏噪音,你可以做到
plot(x[ds$cluster > 0,], col=ds$cluster[ds$cluster>0])