非层次聚类方法

4683 字
23 分钟
非层次聚类方法

5. 非层次聚类方法(Non-hierarchical clustering methods)#

非层次聚类方法从以下二者之一出发:

  • 将对象初始划分为若干组,或
  • 一组初始种子点,它们将构成各簇的核心。

一旦指定了要形成的簇数,对象便会被分配到各簇中。

Important

六簇解并不是由七簇解中的两个簇简单组合而成,而是通过直接寻找最优的六簇解得到的。

5.1 K-means 方法#

5.1.1 步骤#

  1. 将对象初始划分为 KK 个簇。
    • 随机分配或预先指定。
  2. 计算各簇的质心。
  3. 逐一遍历各对象,将每个对象分配到质心(均值)距其最近的簇。
  4. 重复第 2、3 步,直到不再发生重新分配为止。

对象到簇的最终分配在一定程度上取决于初始划分或种子点的初始选取。

5.2 K-means 示例#

对四个对象测量两个变量 X1X_1 和 X2X_2:

对象x1x_1x2x_2
A53
B-11
C1-2
D-3-2

目标是将这些对象划分为

K=2K=2

个簇。

第 1 步:初始划分#

任意地将对象划分为:

  • 簇 1:(AB)(AB)
  • 簇 2:(CD)(CD)

初始质心为:

簇成员xˉ1\bar{x}_1xˉ2\bar{x}_2
1(AB)(AB)22
2(CD)(CD)-1-2
质心

质心(centroid)就是簇内所有成员在每个变量上的均值。

第 2 步:到初始质心(centroids)的距离#

对于对象 A:

d2(A,(AB))=(5−2)2+(3−2)2=10,d^2(A,(AB)) = (5-2)^2+(3-2)^2 = 10,d2(A,(CD))=(5+1)2+(3+2)2=61.d^2(A,(CD)) = (5+1)^2+(3+2)^2 = 61.

A 距簇 1 更近,因此不重新分配。

对于对象 B:

d2(B,(AB))=(−1−2)2+(1−2)2=10,d^2(B,(AB)) = (-1-2)^2+(1-2)^2 = 10,d2(B,(CD))=(−1+1)2+(1+2)2=9.d^2(B,(CD)) = (-1+1)^2+(1+2)^2 = 9.

B 距簇 2 更近,因此将 B 重新分配到簇 2。

所有对象到两个质心的平方距离为:

簇质心ABCD
(AB)(AB)10101741
(CD)(CD)61944

C 和 D 距簇 2 更近,因此它们留在簇 2 中。

更新后的质心#

簇变为:

  • 簇 1:A
  • 簇 2:(BCD)(BCD)

它们的质心为:

簇成员xˉ1\bar{x}_1xˉ2\bar{x}_2
1A53
2(BCD)(BCD)-1-1

新的平方距离为:

簇质心ABCD
A0404189
(BCD)(BCD)52455

不再发生进一步的重新分配。

Note

最终的 K=2K=2 个簇为:

  • 簇 1:A
  • 簇 2:(BCD)(BCD)

6. 簇的数目#

选择簇的数目需要在以下二者之间权衡:

  • 中层管理者,他们可能偏好数量较少、范围较宽的战略性细分市场(KK 小),以及
  • 销售人员,他们可能偏好数量较多、界定狭窄的细分市场(KK 大)。

讲义中讨论的准则包括:

  • ESSESS
  • 簇间距离
  • 方差比准则(Variance-ratio criterion)
  • 轮廓系数(Silhouette coefficient)
  • 簇的相对大小

6.1 方差比准则(Variance-ratio criterion)#

方差比准则为

VRC=BSSESS×N−kk−1.VRC = \frac{BSS}{ESS} \times \frac{N-k}{k-1}.
  • NN :对象总数,整个数据集里有多少个观测对象。
    • 以上面的例子为例,例子里有 A、B、C、D 四个对象,所以 N=4N=4.
  • kk :簇的个数,即 K-means 里设定的 KK.
  • BSSBSS :组间平方和(Between-group Sum of Squares),度量各簇之间分得有多开。
  • ESSESS :组内平方和(Error/Within-cluster Sum of Squares),度量每个簇内部有多紧凑。

其中,

BSS=∑i=1mni∑j=1p(xˉij−xˉj)2,ESS=∑所有对象∑j=1p(xj−xˉ所在簇j)2\begin{aligned} BSS &= \sum_{i=1}^{m} n_i \sum_{j=1}^{p} \left( \bar{x}_{ij}-\bar{x}_j \right)^2, \\ ESS &= \sum_{\text{所有对象}}\sum_{j=1}^p(x_j-\bar x_{\text{所在簇} j})^2 \end{aligned}

式中:

  • xˉij\bar{x}_{ij} 是第 ii 个簇中第 jj 个变量的均值,
  • xˉj\bar{x}_j 是第 jj 个变量的总均值。
  • nin_i :第 ii 个簇的成员个数

6.1.1 整个 VRC 公式的含义#

VRC=BSSESS⏟组间/组内 变异比×N−kk−1⏟自由度修正VRC = \underbrace{\frac{BSS}{ESS}}_{\text{组间/组内 变异比}} \times \underbrace{\frac{N-k}{k-1}}_{\text{自由度修正}}
  • BSSESS\frac{BSS}{ESS}:就是 ANOVA 里 F 统计量的核心——组间差异相对组内差异越大,说明簇分得越 “干净利落”。
  • N−kk−1\frac{N-k}{k-1}:自由度修正。k−1k-1 是 BSSBSS 的自由度(kk 个质心相对总均值的变异),N−kN-k 是 ESSESS 的自由度(NN 个对象相对各自簇质心的变异)。乘上它才构成真正的 F 型统计量,使不同 kk 之间可比。
变异

所谓的“变异”,指的就是“参差不齐的程度”:一组数围绕某个参照点散得有多开。统计里度量它的标准做法是:

变异=∑(每个值−参照点)2\text{变异}=\sum(\text{每个值}-\text{参照点})^2

为什么平方?因为单纯算差值的话,正负会抵消(比如 +3 和 -3 加起来是 0,但明明散得很开)。平方之后全部为正,离参照点越远,惩罚越重。

所以,我们上面说“kk 个质心相对总均值的变异”,翻译成人话就是,各个簇的中心点,离全体数据的中心点有多远(平方后加起来)。

自由度

自由度的含义是:在约束条件下,还能”自由变动”的数值个数。

我们可以用一个简单的例子来理解。假设我们有 3 个数,均值是 10,也就是说,总和固定为 30.

  • 第一个数你可以随便取,比如 7 → 自由
  • 第二个数也可以随便取,比如 2 → 自由
  • 第三个数呢?没得选了,必须是  → 不自由

所以”3 个数 + 均值已知”这个约束下,只有 2 个数是自由的,自由度 = 3 - 1 = 2。

回到 VRC 的两个自由度,BSSBSS 的自由度是 k−1k-1,我们可以想一想, BSSBSS 用到了 kk 个质心,但这些质心是围绕总均值来算偏差的。关键事实是:所有质心(按簇大小加权)的平均恰好等于总均值——这是一个约束。和上面的例子一样,kk  个质心里只有 k−1k-1  个能”自由跑”,最后一个被总均值锁死了。

而对于 ESSESS,它的自由度是 N−kN-k:ESSESS 用到了全部  个对象,但每个对象是和自己簇的质心比较。每个簇内部都有一个约束:簇内成员的偏差之和必须为零(因为质心就是它们的均值)。一共有 kk 个簇,就有 kk 个约束,所以自由度 =N−k=N-k.

那么,为什么 VRC 要除以自由度?

这是因为 BSSBSS 和 ESSESS 都是”加总”出来的量,加总的项越多,数值天然就越大,直接比不公平。类比一下:

  • 3 个质心产生的 BSSBSS 天然比 2 个质心的大——不是聚类更好,只是”分母项”多了;
  • 100 个对象产生的 ESSESS 天然比 4 个对象的大。

除以自由度相当于算”平均每一份自由变动贡献了多少变异”,把它变成一个”单位化”的公平比值。这正是 F 统计量的标准构造:

F=BSS/(k−1)ESS/(N−k)=VRCF=\frac{BSS/(k-1)}{ESS/(N-k)}=VRC

可以看到,VRC 只是把分子分母调换写法,本质一模一样。也正因为做了这个校正,我们才能拿k=2k=2 的 VRC 和 k=5k=5 的 VRC 直接比较大小——这正是 6.4 节那张表里”比较不同簇数的 VRC”做法合法的原因。

💡 记忆:BSS 看簇与簇之间(越大越好),ESS 看簇内部(越小越好),VRC 是两者之比再做自由度校正。

6.2 轮廓系数(Silhouette coefficient)#

对于对象 ii,

s(i)={bi−aimax⁡(ai,bi),ni>1,0,ni=1.s(i) = \begin{cases} \dfrac{b_i-a_i}{\max(a_i,b_i)}, & n_i>1,\\[6pt] 0, & n_i=1. \end{cases}

式中:

  • aia_i = 对象 ii 到同一簇中所有其他对象的平均距离,“我和我的队友有多近”;
  • bib_i = 对象 ii 到最近邻簇中所有对象的平均距离,“我和隔壁有多近”。

我们可以设想两种情形:

情形 1:ni>1n_i > 1(簇里不止自己一个人)#

s(i)=bi−aimax⁡(ai,bi)s(i) = \frac{b_i - a_i}{\max(a_i, b_i)}

分子 bi−aib_i - a_i 是”离隔壁的距离”减去”离队友的距离”:

  • bi≫aib_i \gg a_i:队友很近、隔壁很远 → 分对了!分子为正且大
  • bi≈aib_i \approx a_i:两边差不多远 → 我站在边界上,分到哪边都说得过去
  • bi<aib_i < a_i:隔壁反而比队友更近 → 我可能被分错簇了!分子为负

分母为什么除以 max⁡(ai,bi)\max(a_i, b_i)? 这是为了把结果归一化到 [−1,1][-1, 1] 区间。可以验证两种极端:

  • 完美聚类:ai→0a_i \to 0(紧贴队友),此时 s(i)=bi−0bi=1s(i) = \dfrac{b_i - 0}{b_i} = 1
  • 完全分错:bi→0b_i \to 0,此时 s(i)=0−aiai=−1s(i) = \dfrac{0 - a_i}{a_i} = -1

这样不同对象、不同数据集的分数就可以直接比较了。

情形 2:ni=1n_i = 1(簇里只有自己)#

如果簇里只有一个对象,aia_i 根本没法算(没有”队友”),所以约定 s(i)=0s(i) = 0——表示”无法评价”,既不表扬也不批评。


平均轮廓分数满足

−1≤sˉ≤1.-1 \leq \bar{s} \leq 1.

解释:

  • 取值接近 −1-1:被错误分类,
  • 取值接近 11:聚类良好。

也就是说,轮廓系数回答的问题就是:对象 ii 待在现在这个簇里,合不合适?

6.3 示例:VRC 与轮廓系数#

使用上一示例中最终的 K=2K=2 解:

上图的 1,2,2,21,2,2,2 表示 A 属于簇 1,BCD 属于簇 2.

总均值为

xˉ1=0.5,xˉ2=0.\bar{x}_1=0.5, \qquad \bar{x}_2=0.

则

BSS=(5−0.5)2+(3−0)2+3(−1−0.5)2+3(−1)2=39,BSS = (5-0.5)^2 + (3-0)^2 + 3(-1-0.5)^2 + 3(-1)^2 = 39,

上述四项的结构是:

(5−0.5)2⏟簇1, X1+(3−0)2⏟簇1, X2+3(−1−0.5)2⏟簇2, X1, 乘 n2+3(−1)2⏟簇2, X2, 乘 n2\underbrace{(5-0.5)^2}_{\text{簇1, }X_1} + \underbrace{(3-0)^2}_{\text{簇1, }X_2} + \underbrace{3(-1-0.5)^2}_{\text{簇2, }X_1\text{, 乘 }n_2} + \underbrace{3(-1)^2}_{\text{簇2, }X_2\text{, 乘 }n_2}

且

ESS=(5−5)2+(3−3)2+(−1+1)2+(1+1)2+(1+1)2+(−2+1)2+(−3+1)2+(−2+1)2=14.ESS = (5-5)^2 + (3-3)^2 + (-1+1)^2 + (1+1)^2 + (1+1)^2 + (-2+1)^2 + (-3+1)^2 + (-2+1)^2 = 14.

因此,

VRC=3914×4−22−1=5.5714.VRC = \frac{39}{14} \times \frac{4-2}{2-1} = 5.5714.

四个对象之间的欧氏距离为:

ABCD
A06.32466.40319.4340
B6.324603.60563.6056
C6.40313.605604.0000
D9.43403.60564.00000

轮廓系数的计算如下:

iiaia_ibib_imax⁡(ai,bi)\max(a_i,b_i)s(i)s(i)
A07.38727.38720
B3.60566.32466.32460.4299
C3.80286.40316.40310.4061
D3.80289.43409.43400.5969

从而,

sˉ=0+0.4299+0.4061+0.59694=0.3582.\bar{s} = \frac{ 0+0.4299+0.4061+0.5969 }{4} = 0.3582.

6.4 示例 1:顾客支出细分#

一项针对 24 名顾客的调查记录了其在九个类别上的支出百分比:

变量类别
X1X_1食品
X2X_2交通
X3X_3生活
X4X_4通讯
X5X_5娱乐
X6X_6其他
X7X_7衣着
X8X_8教育
X9X_9储蓄

目标是通过聚类分析进行细分。

6.4.1 方法#

共有九个连续属性。

  • 采用欧氏距离平方作为相异性度量。
  • 对观测值进行标准化。
  • 所考虑的层次凝聚方法:
    • 单连接法,
    • Ward 法。

6.4.2 标准化示例#

原始数据为

对于观测值 1,

zfood=3.3−15.9110.73=−1.176,z_{\text{food}} = \frac{3.3-15.91}{10.73} = -1.176,

且

ztransport=0.9−1.271.00=−0.365.z_{\text{transport}} = \frac{0.9-1.27}{1.00} = -0.365.
标准化公式
z=原始值−该变量的均值该变量的标准差z=\frac{\text{原始值}-\text{该变量的均值}}{\text{该变量的标准差}}

对于 z-score 的含义,以上面的 -1.176 为例,它读作“这位顾客的食品支出比平均水平低 1.176 个标准差。”

对于观测值 1 和 2,标准化空间中的欧氏距离平方为

(−1.176+0.626)2+(−0.365+1.163)2+⋯+(0.362−0.273)2=4.74.(-1.176+0.626)^2 + (-0.365+1.163)^2 + \cdots + (0.362-0.273)^2 = 4.74.

可以看出:

  • 观测值 2 与 16 属于距离最近的一对,
  • 观测值 18 与 22 也属于距离最近的一对。

6.4.3 单连接法#

凝聚过程开始:

  • 第 1 阶段:观测值 2 与 16 形成一个簇,距离 =2.2562=2.2562。
  • 第 2 阶段:观测值 18 与 22 形成一个簇,距离 =3.3621=3.3621。
  • 第 3 阶段:观测值 1 与 3 形成一个簇,距离 =3.4443=3.4443。
  • 第 4 阶段:簇 (1,3)(1,3) 与簇 (2,16)(2,16) 合并,距离 =3.6693=3.6693。
  • ……

簇数准则为:

簇的数目2345
VRC2.10481.82854.37513.6540
轮廓系数0.19830.06940.16310.0595

讲义中的观察:

  • VRC 最大:4 簇解。
  • 轮廓系数最大:2 簇解。
  • 距离图显示从 4 簇到 3 簇、或从 2 簇到 1 簇处出现跳跃。
  • 讲义建议取 4 个簇。

这里简单介绍一下怎么读上面的图。

首先是 height 的含义:height = 该阶段合并两个簇时的距离(相异度),其实就是树状图的纵轴。凝聚聚类每一步都把”距离最近的两个簇”合并成一个,height 就记录这次合并花了多大的”代价”:

  • height 小 → 合并的两个簇本来就很像,合并顺理成章
  • height 大 → 被硬凑在一起的两个簇其实差异很大,合并很勉强

其次就是图的读法。读法是从右往左看:阶段越靠后,剩下的簇越少。


假定分为四个簇:

  • C1=(1,2,3,4,5,6,7,8,9,10,11,12,13,14,16,17)C_1=(1,2,3,4,5,6,7,8,9,10,11,12,13,14,16,17)
  • C2=(15)C_2=(15)
  • C3=(18,19,20,21,22,23)C_3=(18,19,20,21,22,23)
  • C4=(24)C_4=(24)

6.4.4 Ward 法#

凝聚过程包括:

  • 第 1 阶段:观测值 2 与 16 形成一个簇,
d(2,16)=2.26.d(2,16)=2.26.
  • 第 2 阶段:观测值 18 与 22 形成一个簇,
d(18,22)=3.36.d(18,22)=3.36.
  • ……
  • 第 5 阶段:簇 (1,3)(1,3) 与观测值 8 合并成一个簇。

讲义中展示的 Ward 更新公式为

d(13),8=(N8+N1)d18+(N8+N3)d38−N8d13N1+N3+N8.d_{(13),8} = \frac{ (N_8+N_1)d_{18} + (N_8+N_3)d_{38} - N_8d_{13} }{ N_1+N_3+N_8 }.

代入所列数值,

d(13),8=(1+1)(4.62)+(1+1)(4.52)−1(3.44)1+1+1=4.95.d_{(13),8} = \frac{ (1+1)(4.62) + (1+1)(4.52) - 1(3.44) }{ 1+1+1 } = 4.95.

簇数准则为:

簇的数目2345
VRC7.86986.95387.28536.7826
轮廓系数0.25520.22430.21680.2133

讲义指出:

  • VRC 最大:2 簇解。
  • 轮廓系数最大:2 簇解。
  • 距离图显示从 4 簇到 3 簇出现相对较大的增幅。

为了与单连接法解对比,讲义还给出了一种四簇的 Ward 解:

  • C1=(1,2,3,4,8,9,11,13,16)C_1=(1,2,3,4,8,9,11,13,16)
  • C2=(5,12,14,18,22,23,24)C_2=(5,12,14,18,22,23,24)
  • C3=(6,7,10,15,17)C_3=(6,7,10,15,17)
  • C4=(19,20,21)C_4=(19,20,21)
Important

单连接法与 Ward 法的解并不相同。讲义指出,单连接法倾向于形成较大的簇。

6.4.5 四个 Ward 簇的解释#

讲义将标准化均值轮廓总结如下:

  • C1: 在”其他、衣着、教育”方面支出占比较大;在”食品、交通”方面支出占比较小。
  • C2: 在”食品”方面支出占比较大;在”生活、通讯”方面支出占比较小。
  • C3: 在”娱乐”方面支出占比较大;在”教育、储蓄”方面支出占比较小。
  • C4: 在”交通、生活、通讯、储蓄”方面支出占比较大;在”其他、衣着”方面支出占比较小。

6.5 示例 2:顾客面部特征#

一个下午内共有 12 名顾客进入一家商店。

目标是基于五个二元面部特征对顾客进行细分:

  • 性别
  • 眼镜
  • 胡须
  • 微笑
  • 帽子

6.5.1 方法#

共有五个二元属性。

采用简单失配法(Simple mismatching):

  • 对”是”与”否”两种取值均统计失配数,
  • 并将该计数作为相异性度量。

采用完全连接法的层次凝聚方法。

邻近性度量为五个属性中失配特征的个数。

例如:

  • 观测值 1 与 2 有两处失配,故距离 =2=2,
  • 观测值 1 与 3 有一处失配,故距离 =1=1。

对于完全连接法,簇间距离由分属两簇、彼此相距最远的两个观测值确定。

6.5.2 凝聚过程#

讲义重点提到:

  • 第 1 阶段:观测值 1 与 7 形成一个簇;失配数 =0=0。
  • 之后:观测值 3 并入簇 (1,7)(1,7);失配距离 =1=1。

6.5.3 簇数目的选择#

簇的数目3456
VRC13.0316.7517.7915.55

VRC 在五簇时最大,但距离图显示从三簇到两簇出现大幅跳跃。

因此讲义建议采用三簇解:

  • C1=(1,3,7,9)C_1=(1,3,7,9)
  • C2=(2,4,5,8,10,12)C_2=(2,4,5,8,10,12)
  • C3=(6,11)C_3=(6,11)

6.5.4 簇轮廓#

簇% 男性% 眼镜% 胡须% 微笑% 帽子
C110010025750
C25000500
C31000100050

讲义给出的解释:

  • C1: 男性、戴眼镜、不戴帽子。
  • C2: 不戴眼镜、无胡须、不戴帽子。
  • C3: 男性、不戴眼镜、有胡须、不微笑。

6.6 示例 3:美国公用事业公司#

收集了 22 家美国公用事业公司 1975 年的数据。

变量如下:

变量说明
X1X_1固定费用保障比率(收益/债务)
X2X_2资本回报率
X3X_3单位装机容量成本(每 KW)
X4X_4年负荷率
X5X_51974 至 1975 年高峰 KWH 需求增长
X6X_6销售量(每年 KWH 用电量)
X7X_7核电百分比
X8X_8燃料总成本(美分/KWH)

目标是通过聚类分析对这些公司进行定位。

6.6.1 方法#

共有八个连续属性。

  • 采用欧氏距离平方作为相异性度量。
  • 对观测值进行标准化。
  • 采用 K-means 方法。

6.6.2 簇的数目#

kk23456
VRC5.60936.09134.58135.51134.9967
轮廓系数0.18530.22000.11690.19470.1343

讲义接下来在 5 个簇的假定下继续分析:

6.6.3 簇轮廓#

  • 簇 1

    • 7 家公司。
    • 在资本回报率 (X2)(X_2) 上最高。
    • 在成本 (X3)(X_3) 上最低。
  • 簇 2

    • 6 家公司。
    • 在核电 (X7)(X_7) 上最高。
  • 簇 3

    • 1 家公司。
    • 在固定费用保障比率 (X1)(X_1) 与燃料成本 (X8)(X_8) 上最高。
    • 在资本回报率 (X2)(X_2)、年负荷率 (X4)(X_4)、高峰 KWH 需求增长 (X5)(X_5) 与销售量 (X6)(X_6) 上最低。
  • 簇 4

    • 5 家公司。
    • 在年负荷率 (X4)(X_4) 与燃料成本 (X8)(X_8) 上最高。
    • 在固定费用保障比率 (X1)(X_1)、资本回报率 (X2)(X_2) 与核电 (X7)(X_7) 上最低。
  • 簇 5

    • 3 家公司。
    • 在成本 (X3)(X_3)、高峰 KWH 需求增长 (X5)(X_5) 与销售量 (X6)(X_6) 上最高。
    • 在固定费用保障比率 (X1)(X_1)、资本回报率 (X2)(X_2)、核电 (X7)(X_7) 与燃料成本 (X8)(X_8) 上最低。

6.6.4 与 Ward 法的比较#

讲义指出,对于这一五簇解,K-means 与 Ward 法的划分结果一致。

采用讲义中给出的 Ward 簇标签:

C1:南方公司群体#

  • Oklahoma Gas and Electric Company
  • Texas Utilities Company
  • Arizona Public Service
  • The Southern Company
  • Central Louisiana Electric Co.
  • Kentucky Utilities Co.
  • Florida Power & Light Co.

C2:东部+中部公司群体#

  • Madison Gas & Electric Co.
  • Northern States Power Co.
  • Commonwealth Edison Co.
  • Wisconsin Electric Power Co.
  • Boston Edison Co.
  • Virginia Electric & Power Co.

C3:独立公司#

  • Consolidated Edison Co. (New York)

C4:东西部近海岸公司#

  • New England Electric Co.
  • United Illuminating Co.
  • Pacific Gas & Electric Co.
  • Hawaiian Electric Co.
  • San Diego Gas & Electric Co.

C5:西北部公司群体#

  • Idaho Power Co.
  • Puget Sound Power & Light Co.
  • Nevada Power Co.

6.6.5 地理解释#

讲义的结论是:公用事业公司往往按地理位置聚成群体或展开竞争。

  • C4:东部与西部,但靠近海岸。
  • C2:东部+中部。
  • C1:南部。
  • C5:西北部。
  • Consolidated Edison 独自成簇。

讲义给出的可能原因:

  • 位置相近的公司往往在电厂中使用相似的燃料,并面对共同的市场,
  • 因此这些公司在发电类型、成本、增长率及相关变量上应相对同质。

7. 与聚类的关联分析#

7.1 单因素方差分析模型(1-way ANOVA model)#

单因素方差分析研究以下两者之间的关系:

  • 一个定量因变量(quantitative dependent variable),以及
  • 一个定性自变量(qualitative independent variable)。

设:

  • rr = 因子的水平数(Number of levels of factor),
  • nTn_T = 案例总数(Total number of cases)。

假设为

H0:μ1=μ2=⋯=μr,H_0: \mu_1 = \mu_2 = \cdots = \mu_r,H1:not all μi are equal.H_1: \text{not all }\mu_i\text{ are equal}.

当下式成立时拒绝 H0H_0:

F>Fα,r−1,nT−rF > F_{\alpha,r-1,n_T-r}

显著性水平为 α\alpha。

7.2 χ2\chi^2 检验#

χ2\chi^2 检验研究两个定性变量之间的关系。

设:

  • rr = 变量 A 的水平数,
  • cc = 变量 B 的水平数。

假设为

H0:There is no association between the two variables,H_0: \text{There is no association between the two variables},H1:There is an association between the two variables.H_1: \text{There is an association between the two variables}.

当下式成立时拒绝 H0H_0:

χ2>χα,(r−1)(c−1)2\chi^2 > \chi^2_{\alpha,(r-1)(c-1)}

显著性水平为 α\alpha。

7.3 示例:与顾客簇的关联分析#

讲义提出以下问题:

  1. 各聚类变量在不同簇之间是否存在差异?
  2. 聚类结果与顾客类型之间是否存在关联?其中:
    • shop =1=1:百货商店,
    • shop =2=2:超市,
    • shop =3=3:本地杂货店?

7.3.1 单因素方差分析#

对标准化后的聚类变量进行单因素方差分析。

临界值为

F0.05,3,20=3.098.F_{0.05,3,20} = 3.098.
变量FFp-value
食品18.7<0.001<0.001
交通7.40.0016
生活4.20.0186
通讯4.590.0133
娱乐5.30.0075
其他8.46<0.001<0.001
衣着19.9<0.001<0.001
教育13.1<0.001<0.001
储蓄1.960.1528

[!result] 在 5% 显著性水平下,除储蓄外的所有聚类变量在四个簇之间均存在显著差异。

7.3.2 χ2\chi^2 检验#

带列百分比的二维频数表如下:

顾客类型簇 1簇 2簇 3簇 4
百货商店6 (66.67%)1 (14.29%)2 (40.00%)0 (0%)
超市3 (33.33%)2 (28.57%)3 (60.00%)0 (0%)
本地杂货店0 (0%)4 (57.14%)0 (0%)3 (100%)

检验统计量为

χ2=17.4>12.59=χ0.05,62.\chi^2 = 17.4 > 12.59 = \chi^2_{0.05,6}.

并且,

p-value=0.0076<0.05.p\text{-value} = 0.0076 < 0.05.

因此,商店类型与聚类结果之间存在关联。

讲义对各簇的刻画如下:

  • 簇 1:百货商店
  • 簇 2:本地杂货店
  • 簇 3:超市
  • 簇 4:本地杂货店

文章分享

如果这篇文章对你有帮助,欢迎分享给更多人!

非层次聚类方法
https://hac.he.cn/posts/stat-04/
作者
HAC
发布于
2026-09-29
许可协议
CC BY-NC-SA 4.0

评论区

Profile Image of the Author
HAC
观之非易,行且克难
Greetings
欢迎来到我的博客!这里主要分享我的学习笔记与兴趣爱好。
音乐
封面

音乐

暂未播放

0:000:00
暂无歌词
分类
标签
站点统计
文章
46
分类
8
标签
19
总字数
116,126
运行时长
0 天
最后活动
0 天前
文章目录