多案例 Cook 距离

必需输入

• 单列响应值
• 多列预测变量值。

可选输入

运行宏

用于运行宏的语法会略有区别，具体取决于您使用的版本。

下面的示例使用的样本数据是 Rousseeuw 和 Leroy (1987) 中的“Modified Data on Wood Specific Gravity”（木材特定比重的修改数据）数据集，其中包括二十个案例和五个预测变量。五个选定案例子集的计算结果与 Seaver、Triantis 和 Reeves (1999) 中给出的结果一致。

假设响应值 Y，即特定比重，在 C1 中，五个预测变量的值 X1-X5 在列 2-6 中。五个子集案例已选定。

%MULTDIST C1-C6;
SUB1 5;
SUB2 8 19;
SUB3 6 8 19;
SUB4 4 8 19;
SUB5 4 6 8 19。

单击运行。

输出

下面是该宏生成的结果。

Multiple Case Cook's Distance

Model Information
------------------------
Response:     Y

Predictors:   X1 , X2 , X3 , X4 , X5                                            

Parameters:    6
 
Threshold value:    1.00
------------------------
 
*** Cook's Distance for Case Pairs ***
 
     Cases        Cook's Distance

     7 , 11             1.03

 
*** Cook's Distance for a Subset ***

     Cases:  5   Cook's Distance:  0.06                                              


     Cases:   8 , 19   Cook's Distance:  0.33                                        


     Cases:   6 ,  8 , 19   Cook's Distance:  1.99                                   


     Cases:   4 ,  8 , 19   Cook's Distance:  0.49                                   


     Cases:   4 ,  6 ,  8 , 19   Cook's Distance:  53.93 

更多信息

数据集大小

对于案例对和案例三要素，计算 Cook 距离的数据集大小限制分别为 60 和 30。计算案例子集的数据集大小限制为 500。您可以在此宏内更改案例对和三要素限制。要更改这些限制，请转到宏代码中标记为“MSE 检查，三要素，无配对”部分，然后将 30 和 60 更改为您所需的大小。请注意，计算所需的时间会随着数据集大小的增加而增加，尤其是计算所有三要素时更是如此。

逆矩阵不存在

如果分析混合模型，必须指定非常数子命令。如果没有指定非常数子命令，将会看到一则错误消息，指示 XTX 矩阵的逆矩阵不存在。通常，如果有任何预测变量（几乎）完全相关，则会看到此错误消息。

缺失值

该宏通过删除包含缺失数据的行来处理缺失数据。这会在输出和图形中显示。

参考资料

Rousseeuw, P. J. 和 Leroy, A. M. (1987)，《Robust Regression & Outlier Detection》（可靠回归与异常值检测），John Wiley & Sons, Inc.

Seaver, B.、Triantis, K. 和 Reeves, C. (1999)，The Identification of Influential Subsets in Regression Using a Fuzzy Clustering Strategy（在回归中使用模糊聚类策略识别有影响的子集）,《Technometrics》（技术计量学），第 41 期，第 340 到 351 页。