上述代码j未定义

非常抱歉，可能是因为我在前面没有定义一些必要的库和函数导致代码无法运行。下面是完整的代码：

import numpy as np
import matplotlib.pyplot as plt
from sklearn.datasets import load_iris
from scipy.spatial.distance import pdist, squareform

# 定义 DIANA 聚类函数
def diana(X, k):
    dists = squareform(pdist(X))
    np.fill_diagonal(dists, np.inf)
    labels = np.zeros(X.shape[0])
    labels[:k] = np.arange(k)
    for i in range(k, X.shape[0]):
        idx = int(np.argmin(np.mean(dists[np.ix_(labels == j, labels == j)], axis=0)))
        labels[i] = labels[labels == idx][0]
        labels[labels == idx] = i
    return labels.astype(int)

# 加载鸢尾花数据集
iris = load_iris()
X = iris.data
y = iris.target

# 使用 DIANA 算法进行聚类
k = 3
labels = diana(X, k)

# 可视化聚类结果
colors = ['r', 'g', 'b']
for i in range(k):
    plt.scatter(X[labels == i, 0], X[labels == i, 1], c=colors[i], label='Cluster {}'.format(i+1))
plt.title('DIANA Clustering of Iris Dataset')
plt.xlabel('Sepal Length')
plt.ylabel('Sepal Width')
plt.legend()
plt.show()

希望这次可以成功运行。