Python SciPy hierarchy.maxRstat用法及代碼示例

本文簡要介紹 python 語言中 scipy.cluster.hierarchy.maxRstat 的用法。

用法:

scipy.cluster.hierarchy.maxRstat(Z, R, i)#

返回每個非單例集群及其子集群的最大統計量。

參數 ：：

Z： array_like

層次聚類編碼為矩陣。有關詳細信息，請參閱 linkage 。

R： array_like

不一致矩陣。

i： int

用作統計量的 R 列。

返回 ：：

MR： ndarray

計算不一致矩陣的第 i 列的最大統計量R對於每個非單例集群節點。MR[j]是最大超過R[Q(j)-n, i]，其中Q(j)與以下節點相對應的所有節點 ID 的集合，包括j.

例子：

>>> from scipy.cluster.hierarchy import median, inconsistent, maxRstat
>>> from scipy.spatial.distance import pdist

給定數據集 X ，我們可以應用聚類方法來獲得鏈接矩陣 Z 。 scipy.cluster.hierarchy.inconsistent 也可用於獲取與此聚類過程相關的不一致矩陣R：

>>> X = [[0, 0], [0, 1], [1, 0],
...      [0, 4], [0, 3], [1, 4],
...      [4, 0], [3, 0], [4, 1],
...      [4, 4], [3, 4], [4, 3]]

>>> Z = median(pdist(X))
>>> R = inconsistent(Z)
>>> R
array([[1.        , 0.        , 1.        , 0.        ],
       [1.        , 0.        , 1.        , 0.        ],
       [1.        , 0.        , 1.        , 0.        ],
       [1.        , 0.        , 1.        , 0.        ],
       [1.05901699, 0.08346263, 2.        , 0.70710678],
       [1.05901699, 0.08346263, 2.        , 0.70710678],
       [1.05901699, 0.08346263, 2.        , 0.70710678],
       [1.05901699, 0.08346263, 2.        , 0.70710678],
       [1.74535599, 1.08655358, 3.        , 1.15470054],
       [1.91202266, 1.37522872, 3.        , 1.15470054],
       [3.25      , 0.25      , 3.        , 0.        ]])

scipy.cluster.hierarchy.maxRstat 可用於計算每個非單例集群及其子集群的 R 每列的最大值：

>>> maxRstat(Z, R, 0)
array([1.        , 1.        , 1.        , 1.        , 1.05901699,
       1.05901699, 1.05901699, 1.05901699, 1.74535599, 1.91202266,
       3.25      ])
>>> maxRstat(Z, R, 1)
array([0.        , 0.        , 0.        , 0.        , 0.08346263,
       0.08346263, 0.08346263, 0.08346263, 1.08655358, 1.37522872,
       1.37522872])
>>> maxRstat(Z, R, 3)
array([0.        , 0.        , 0.        , 0.        , 0.70710678,
       0.70710678, 0.70710678, 0.70710678, 1.15470054, 1.15470054,
       1.15470054])