非参数统计__秩相关分析和秩回归.ppt

第七章 秩相关分析和秩回归;相关系数的度量;7.1 Spearman秩相关系数及检验;检验;例7.1;解答; 相关系数及检验 ;这样的样本共有 个数对，用 表示协同的数对的数目， 表示不协同的数对数目。则 系数定义为： 其中 ，易知;另一种转换形式: 将X的数据由小到大排序, 由于协同性考虑Y的秩, 记为: d1,d2,…,dn, 计算;例7.2;d1,d2,……, d10;1 0; x-c(75,95,85,70,76,68,60,66,80,88) y-c(2.62,2.91,2.94,2.11,2.17,1.98,2.04,2.2,2.65,2.69) cor.test(x,y,meth=kendall);练习: 双胞胎儿童间的智力相关程度分析。 某幼儿园对9对双胞胎的智力进行测验，并按百分制打分。现将资料列示如表 ： ; ;计算Kendall秩相关系数 即双胞胎儿童间的智力相关程度为0.722;多变量Kendall协同系数检验 ;每列的秩和为:;实际检验时，可以查零分布表，在n固定， 时： 拒绝域：{Wc};当样本中有结点时，采用修正的Kendall协和系数 ;例7.3; x1-c(41,43,39.5,38,40.5,41,40,38.5,44,39) x2-c(55.7,56.3,54.5,54.2,55.1,55.4,54.5,54.2,56.9,54.5) x3-c(8.6,9.2,8,5.6,6.8,8,8.6,7.4,9.8,7.4) y1-rank(x1) y2-rank(x2) y3-rank(x3) Rh-y1+y2+y3 SSR-sum(Rh*Rh)-(sum(Rh))^2/10 SSR [1] 657.5 Wc-12*SSR/(9*(10^3-10)-3*(5*(2^3-2)+(3^3-3))) ka-3*(10-1)*Wc qchisq(0.95,9) [1] 16.91898 %查表值 ka [1] 24.35185 %计算值 (拒绝H0, 三个因素一致相关);Kappa一致性检验 ;Kappa一致性检验 ;Kappa一致性检验 ;Kappa一致性检验 ;Kappa一致性检验 ;例 ;解答 ;解答 ;一元线性回归;多元线性回归;多元线性回归系数估计;例;X1=c(-0.05, 0.25,0.60,0, 0.25,0.20, 0.15,0.05,-0.15, 0.15, 0.20, 0.10,0.40,0.45,0.35,0.30, 0.50,0.50, 0.40,-0.05, -0.05,-0.10,0.20,0.10,0.50,0.60,-0.05,0, 0.05, 0.55) X2=c( 5.50,6.75,7.25,5.50,7.00,6.50,6.75,5.25,5.25,6.00, 6.50,6.25,7.00,6.90,6.80,6.80,7.10,7.00,6.80,6.50, 6.25,6.00,6.50,7.00,6.80,6.80,6.50,5.75,5.80,6.80) Y=c( 7.38,8.51,9.52,7.50,9.33,8.28,8.75,7.87,7.10,8.00, 7.89,8.15,9.10,8.86,8.90,8.87,9.26,9.00,8.75,7.95, 7.65,7.27,8.00,8.50,8.75,9.21,8.27,7.67,7.93,9.26) lm.sol-lm(Y~X1+X2) summary(lm.sol);Theil和中位数回归系数估计 ;1）Brown-Mood方法：;例7.5 ;分组程序 ;例7.5 ;例7.5 ;一元线性回归;一元线性回归;1）Theil方法：;Theil方法 ;三种方法的效果图形 ;课后习题7.5;1）BM方法;2）Theil方法;关于α和β的检验问题 ;关于α和β的检验问题 ;对于斜率的检验用如下统计量：;关于α和β的检验问题