图例未显示K-Means聚类

Question

我使用K-means聚类聚类了一些值。我绘制了三个图表，我想要第三个图表的图例。我使用lengend()函数绘制图例，但图例中没有出现图例（见图），并且没有报告错误。

library(cluster) 
library(MASS)

par(mfrow=c(1,3))
par(oma=c(5,0,3,0))

dat <- read.table("representative_for_r.txt", header=TRUE)


data.p <- as.matrix(dat$Average_RMSD)

# Determine number of clusters
wss <- (nrow(data.p)-1)*sum(apply(data.p,2,var))
for (i in 2:15) wss[i] <- sum(kmeans(data.p,centers=i)$withinss)
plot(1:15, wss, type="b", xlab="Number of Clusters",ylab="Within groups sum of squares")



# K-Means Cluster Analysis
fit <- kmeans(data.p, 6) # 6 cluster solution

# get cluster means 
aggregate(data.p,by=list(fit$cluster),FUN=mean)

# append cluster assignment
data.p <- data.frame(data.p, fit$cluster)


# PLot Clusters
clusplot(data.p, fit$cluster, color=TRUE, shade=TRUE, labels=2, lines=0)

data.p <- data.frame(dat$PDB, data.p)
#print(data.p)

plot(data.p[,2],col=data.p$fit.cluster) # takes the RMSD column of data.p(by indexing) then colours the points defined by the clustering
par(xpd=NA)
legend(0,0, c("Group 1","Group 2", "Group 3", "Group 4", "Group 5", "Group 6", "Group 7"), pch=1, col=1:7)

write.matrix(data.p, file = "kmeans_output.txt", sep = " ")

输出如下：

我尝试使用以下内容更改X和Y值：

200,200
-200,200
200,-200
-200,-200

Answer 1

尝试：

legend("bottomleft", legend = paste("Group", 1:7), pch=1, col=1:7)

注意y轴上的比例。传递给legend()的第二个坐标必须位于图的范围内。 0以及-200和200都超出了图3的限制范围。

如果您不想猜测这些内容，则使用字符来定义放置更加直观。有关可提供的字符位置的详细信息，请参阅?legend。