多个变量之间的关系描述
全糖奶茶屋
2022年08月31日 23:22

        我们在研究变量之间的关系的时候,如果能从一个变量的信息中得到另一个变量的信息,那么这两个变量之间就是相关的。在建模题目中, 经常要我们分析某变量与某变量之间的关系, 我们就需要分以下各方面进行分析.

1. 散点图

        研究两个变量之间关系的最简单方法是散点图(scatter plot),但好的散点图的绘制并不简单。

        下面绘制BRFSS调查参与者的体重与身高关系的散点图。但这个绘制结果数据都成列聚集,这是因为身高数据四舍五入到相邻的英寸,转换为厘米后,再次四舍五入。在这个转换过程中,丢失了一些信息。

 BRFSS调查参与者的体重与身高关系散点图,未抖动(左),抖动(右)

        即使经过了抖动处理,散点图也不是展示数据的最佳方法。图中有很多重叠的点,遮盖了密集部分的数据,使离群值显得特别突出。这种效果称为饱和(saturation)。

2. 描述关系特征

        散点图能让我们对变量关系有个大体了解,其他可视化方法则可以让我们更深入了解变量关系的本质。一种方法是对一个变量进行分区,绘制另一个变量的百分位数。

        计算分区的过程如下:

  • 去除指定列含有nan值的数据行,对数据进行分区

  • 遍历每个分区的数据,计算其身高均值和体重CDF

  • 绘制身高对应的体重百分位数

        下图展示了绘制结果,在140-200厘米,变量关系几乎是线性的。140-200厘米这个范围涵盖了超过99%的数据。

一组高度分区的体重百分位数

3. 相关性

        相关性(correlation)是一个统计量,用于量化两个变量之间关系的强弱。

        度量相关性的困难在于需要比较的变量通常使用不同的单位,即使变量使用相同的单位,也可能来自不同的分布。通常有两个常见的解决方法:

  • 将每个值转换为标准分数(standard score),即其偏离均值的标准差数。这种转换会产生“Pearson乘积矩相关系数”。

  • 将每个值都转换为秩,即其在所有值的排列表中的索引。这种转换会产生“Spearman秩相关系数”。

4. Pearson相关性

        协方差在一些计算中非常有用,但其含义很难解释,因此很少作为摘要统计量。解决这个问题的方法之一是将偏差除以标准差,得到标准分数,然后计算标准分数的乘积:

p_i%3D%5Cfrac%7Bx_i-%5Cbar%7Bx%7D%7D%7BS_X%7D%20%5Cfrac%7By_i-%5Cbar%7By%7D%7D%7BS_Y%7D

其中, S_X和S_Y%20分布是X和Y的标准差。

        这些乘积的均值为:

%5Crho%3D%5Cfrac%7B1%7D%7Bn%7D%20%5Csum%20p_i

或者,可以通过分解 S_X 和SyS_Y  改写标准差:

%5Crho%3D%5Cfrac%7BCov(X%2CY)%7D%7BS_X%20S_Y%7D

        这个公式以统计学家Karl Pearson的名字命名,称为Pearson相关性。Pearson相关性取值介于-1~+1之间(包含端点)。

        %5Crho 的大小表明了相关性的强弱程度,如果 %5Crho 为1或-1,两个变量完全相关。

5. 非线性关系

        Pearson相关性只是度量了线性关系,如果变量之间存在非线性关系,那么\rhoρ对变量相关性强弱的估计就可能是错误的。

        下图摘取自Wiki,展示了数据集的散点图和相关系数。

各种相关性的数据集示例

        第一行展示了不同线性相关性的数据集,第二行展示了具有不同斜度的完全相关,第三行展示了变量非线性的相关性。

6. Spearman秩相关

        如果变量之间的关系是线性的,而且变量大致符合正态分布,那么Pearson相关性能很好地说明相关性的强弱,但离群值会影响Pearson相关性的稳健性。Spearman秩相关能缓解离群值以及偏斜分布的影响,也可用于描述变量的相关性。要计算Spearman相关性,必须计算每个值的秩,即该值在排序样本中的索引,然后计算这些秩的Pearson相关性。

7. 相关性和因果关系

        如果变量A和变量B相关,那么有3种可能:A导致B,B导致A,或其他因素导致A和B。这些解释称为“因果关系”。但相关性并不意味着因果关系。

        统计关系,无论多么强、多么富有启示性,都不能确立因果关系。因果关系的思想必须来自于统计学之外,来源于一些理论或者其他方面

——Kendall & Stuart(1961)

        

        证明因果关系的方式有:

  • 时间: 如果A在B之前发生,那么A可能导致B,而B不可能导致A。

  • 随机性: 如果将一个大型样本随机分为两组,计算任意变量的均值,那么两组结果的差别应该很小。

        正是这些想法催生了随机对照试验。在随机对照试验中,试验对象被随机分配到两个(或多个)组:试验组和对照组,试验组接受某些干预,对照组不接受干预。