Cochran检验和Durbin不完全区组分析法.pptVIP

第五章多总体统计推断

Cochran检验

针对定性数据表的检验方法。例如，考察三名推销员的推销能力，顾客对推销员进行打分，满意——1分，不满意——0分。

顾客

123456789101112

甲

乙

丙

111111001110

010100010000

101001010000

问推销员的推销效果是否一样？

Cochran检验

检验

则

检验

需要估计：

检验

CochranQ检验统计量：

Q近似服从分布，当Q值偏大的时候，考虑拒绝零假设。

Durbin不完全区组分析

原理：

可能存在处理非常多，但是每个区组中允许的样本量有限的时候，每一个区组中不可能包含所有的处理，比如重要的均衡不完全区组BIB设计。Durbin检验便是针对这种问题。例如，

不同城市保险公司的绩效的BIB设计（万份）

(b=4)

城市（区组）k=4

ⅠⅡⅢⅣ

A

B

C

D

342859

303645

364448

405460

构造统计量：

当D值较大的时候，可以考虑拒绝零假设，认为处理之间存在差异。在零假设成立时，大样本情况下，D近似服从分布。打结的时候，只要长度不大，对结果影响不太大。

D等同于Friedman统计量。(t=k,r=b)

k=4,b=4,t=3,r=3.

解答

R函数文件

Cochran检验

x1=c(0,0,0,1,0,0,0,0,0,1)

x2=c(1,1,0,1,0,1,0,0,1,1)

x3=c(1,1,1,1,1,1,1,1,1,0)

x=matrix(c(x1,x2,x3),nrow=10,ncol=3)

ni=apply(x,1,sum)

nj=apply(x,2,sum)

k=ncol(x)

Q=(k-1)*((k*sum(nj^2))-(sum(nj))^2)/(k*sum(ni)-sum(ni^2))

pvalue.x=pchisq(Q,k-1,lower.tail=F);pvalue.x;Q

谢谢大家！

ThankYou！