《相关和回归》PPT课件.ppt

,第七章 相关和回归,,第七章 相关和回归,Spearman秩相关检验,1,2,3,5,Kendall秩相关系数检验,Theil 非参数回归和几种稳健回归,,一、适用范围 Spearman 秩相关系数是最早、最著名评秩统计量,主要用于研究两变量间的相关程度及其显著性检验,其资料要求两变量都至少是以定序尺度测量的。 二、理论依据和方法 1理论依据： Spearman秩相关系数，用rs代表。是对容量为n的xi和yi的秩(i=1,2,n)进行相关性测量。如两变量x与y完全正相关，则应有xi=yi；如完全负相关，应有x1=yn，x2=yn-1，xn=y1。Spearman秩相关系数通过di=xiyi研究总的偏离程度。,第一节 Spearman秩相关,,在计算相关系数时，直接研究di是不合适的，因为正的di与负的di相互抵消，因此采用di2，当di越大时，di2也越大。 rs的计算公式为： (7.1) 当di2为0时，rs =1，可认为两个变量完全正相关。rs所量度的是两个等级之间的联系强度，rs处于+1和-1之间。,第一节 Spearman秩相关,,2.显著性检验 假设组Ho：X和Y相互独立(X和Y正或负相关) H1：X与Y相互不独立(X与Y负或正相关) 检验rs的显著性,在小样本情况下，即n从4到30时，可查附表13来检验。该表列出了在H0成立时相伴概率分别为 a0.05和a0.01的rs值。这是一个单尾表适合用于检验单侧假设。即当rs大于或等于表中临界值时，拒绝Ho。,第一节 Spearman秩相关,,在大样本情况下，即n10时，在零假设成立时得到的rs的显著性可用统计量t来检验： t统计量近似服从df=n-2的T 分布。 如n很大时，即n30,还可用统计量Z来检验。 近似服从正态分布。 (7.2) 3.耦合修正。 两变量的秩相等即耦合。这时用它们的平均秩作耦合项的秩。,第一节 Spearman秩相关,,在耦合现象出现的比例不大时，可以忽略它们对rs的影响；但当其比例较大时， rs用下式来修正： (7.3) 式中, Tx=(t3-t)/12,t等于x变量中同一个秩的耦合数;Ty=(t3-t)/12，t代表y变量中耦合的观察数。,第一节 Spearman秩相关,,三、检验步骤 1.据题意，作正确假设； 2.将变量X和Y的观察值分别从1到n评秩，如观察值相同，用平均秩代替。 3.将两样本配对成(xi，yi),xi，yi分别代表两变量的秩。 4.定出di=xi-yi；算出di2。 5.如无耦合现象或比例较小时，用公式(7.1)算出rs，如耦合现象比例较大，则用(7.3)公式计算rs 。 6.小样本4n30时,查附表13;大样本时,N30用公式(7.2) 。,第一节 Spearman秩相关,,四、例7.1(小样本举例)学习时间长短与学生考试成绩间是否有关。调查某大学10个学生每周学习的时间与期末平均成绩的资料如表7-1所示 解：假设H0：学习时间X与平均成绩等级Y之间是相互独立的；H1：学习时间X与平均成绩等级Y之间是正相关。根据(9.22)式计算得到: 取a=0.05的显著性水平，样本容量n10，查附表13临界值rs(n，a)=0.5515,第一节 Spearman秩相关,,表7-1 大学生的学习时间与期末成绩调查表,第一节 Spearman秩相关,,因为rs=0.9460.5515（临界值）故拒绝Ho假设,而接受H1假设，即学生的学习时间与学生的平均成绩等级之间存在着正相关关系。 这里变量x和y中都分别有一对存在耦合即t2，所占比例不大，可以不必进行修正。但为说明方法，在此作一修正。 用公式(7.3), 首先要算出 代入公式得： 修正后比修正前略有减小。,第一节 Spearman秩相关,,例7.2(大样本举例)从一个大企业的生产线上抽出15名雇员组成一个随机样本，然后让这些雇员的同事和管理人员根据他们对工作的兴趣及合作精神分别对他们进行排序，结果列于表7-2，企业当局想知道同事和管理人员的看法是否正相关。 解：假设Ho：同伴对雇员的等级评定与管理人员对雇员的等级评定间相互独立；H1：同伴对雇员的等级评定与管理人员对雇员的等级评定正相关 由(7.1)式算得rs=1-(6148)/(15 (152-1)=0.7357 查附表13得r15，0.050.446rs。因此我们以0.05的显著性水平否定零假设,作出“两种等级正相关”的结论,第一节 Spearman秩相关,,这里n=1510。我们也可用T统计量来检验。 查T分布表当df13，a=0.05时，临界值ta/2=2.160，因t=3.92ta/2，所以我们可以在a=005的显著性水平上拒绝Ho，从而得出结论：同伴和管理人员对这15个雇员的评价正相关。,第一节 Spearman秩相关,,一、适用范围 肯德尔（Kendall）秩相关系数（读为Tao）是测度两组秩间的相关程度。 二、理论依据和方法 在容量为n的成对样本(xi，yi)中，如果某对观察值都比另一对观察值大，则称这两对观察值为一致对(协同)。如(1.3,2.2)和(1.6,2.7)就是一致对(用Nc表示在n 2个可能对中一致对的数目)；反之如果某对观察值与另一对观察值以相反的方向变化(即一个减少而另一个增加)，我们称之为非一致对，并用Nd表示在n 2个可能对中非一致对(不协同)的数目。,第二节 Kendall相关检验,,如果两对观察值相等，就不能算是一致对，也不能算作非一致对。 肯德尔秩相关系数的计算公式为： (7.4) 如所有观察值对都是一致对，则1；如所有观察值对都是非一致对，则1 在计算值时,可以将所有观察对(xi，yi)按变量值xi的递增次序排列,将与其对应的yi与后面的每一个yi比较,有多少对yi是按递增次序排列的,就有多少个一致对；有多少对yi不是递增次序排列,就有多少个非一致对。,第二节 Kendall相关检验,,例9.21: 双胞胎儿童间的智力相关程度分析。 某幼儿园对9对双胞胎的智力进行测验，并按百分制打分。现将资料列示如表7-3： 表7-3 解：根据计算要求，首先将资料（xi，yi）按xi的递增次序排列，并将计算结果列如表7-4所示：,第二节 Kendall相关检验,,表74,第二节 Kendall相关检验,,计算Kendall秩相关系数 即双胞胎儿童间的智力相关程度为0.722 三、显著性检验 求出Nc和Nd后,用统计量T(T= Nc-Nd)检验下面的假设： Ho：X与Y是相互独立的;H1：X与Y不是相互独立的。（正相关或负相关） 再根据样本大小n和显著水平a，可查附表14T的临界值如果TWn(1-a)就接受正相关的备择假设；如果n(a)就接受负相关的备择假设。,第二节 Kendall相关检验,,上例中，T=Nc-Nd=31-5=26。 又因n9，a=0.05，查附表14，得出双侧检验的两个临界值为：W0.97518，W0。02518 而T26W097518，因此以0.05的显著水平拒绝Ho，接受X与Y之间存在正相关的备择假设。 四、耦合处理 在观察对(Xi，Yi)中有两个或两个以上的数值相等时，用平均秩来代替。这种情况只影响(7.4)式中的分母 (7.5) 其中Tx=Ty=1/2 t(t-1),t为观察对中有相同值X(Y)的数目.,第二节 Kendall相关检验,,五、大样本检验 当N10时，遵从如下的正态分布： 均值U=0,标准差： 所以正态变量： 上例中N12于10，可用正态近似法计算 =0.822 查正态分布表，当Z3.72时大于1.96，所以拒绝Ho假设。而接受学生所花的学习时间与取得的学习成绩之间存在相关的结论。,第二节 Kendall相关检验,,六、功效讨论 斯氏秩相关系数rs与肯氏秩相关系数之间对比，两者采用两种不同的方法计算相关系数，所以即使用同样的资料，也会得出不同的数值。在大多数情况下，斯氏的rs的绝对值比肯氏的的绝对值为大（不是所有的情况均如此）。但我们绝不能根据它们计算的数值不同，而得出前者的相关程度密切，