非参数统计秩相关分析和秩回归演示.pptVIP

非参数统计秩相关分析和秩回归;相关系数的度量;7.1Spearman秩相关系数及检验;检验;例7.1;解答;相关系数及检验;这样的样本共有个数对，用表示协同的数对的数目，表示不协同的数对数目。则系数定义为：

其中，易知;另一种转换形式:

将X的数据由小到大排序,由于协同性考虑Y的秩,记为:d1,d2,…,dn,计算;例7.2;d1,d2,……,d10;10;当前13页，总共62页。;练习:双胞胎儿童间的智力相关程度分析。

某幼儿园对9对双胞胎的智力进行测验，并按百分制打分。现将资料列示如表：

;

;计算Kendall秩相关系数

即双胞胎儿童间的智力相关程度为0.722;多变量Kendall协同系数检验;每列的秩和为:;实际检验时，可以查零分布表，在n固定，时：

拒绝域：{Wc};当样本中有结点时，采用修正的Kendall协和系数

;例7.3;当前22页，总共62页。;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一致性检验;例;解答;解答;一元线性回归;当前33页，总共62页。;当前34页，总共62页。;当前35页，总共62页。;当前36页，总共62页。;多元线性回归;当前38页，总共62页。;多元线性回归系数估计;例;当前41页，总共62页。;当前42页，总共62页。;当前43页，总共62页。;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方法;当前59页，总共62页。;关于α和β的检验问题;关于α和β的检验问题;谢谢大家！