R语言中替换掉ch[k] <- cluster.stats(iris_data, km$cluster)$ch中的cluster.stats函数

要替换掉R语言中的cluster.stats函数，您需要编写自己的函数，实现与cluster.stats相同的功能，并将其应用于ch[k]中。以下是一个示例函数，它计算k-means聚类的Calinski-Harabasz指数：

calinski_harabasz <- function(data, clusters) {
  # 计算簇内平均距离
  within_ss <- sum(sapply(unique(clusters), function(i) {
    mean(dist(data[clusters == i,]))^2
  }))
  
  # 计算簇间平均距离
  between_ss <- sum(sapply(unique(clusters), function(i) {
    n_i <- sum(clusters == i)
    mean(dist(data[clusters != i,][1:n_i,]))^2
  }))
  
  # 计算Calinski-Harabasz指数
  ch <- (between_ss / within_ss) * ((sum(table(clusters)) - length(unique(clusters))) / (length(unique(clusters)) - 1))
  
  return(ch)
}

要将此函数应用于ch[k]中，只需将cluster.stats替换为calinski_harabasz：

ch[k] <- calinski_harabasz(iris_data, km$cluster)