Hodges-Lehmmann调整秩和检验（非参数统计习题）

例：下表是美国三大汽车公司（A，B，C三种处理）的5种不同车型某年产品的油耗分析不同公司产品的油耗是否存在差异。

A<-c(20.3,21.2,18.2,18.6,18.5)

B<-c(25.6,24.7,19.3,19.3,20.7)

C<-c(24,23.1,20.6,19.8,21.4)

d<-rbind(A,B,C)        # 构造一个矩阵d

 x<-colMeans(d)     # 求出 矩阵d 的每列（每一个区组）平均

for (i in 1:5) {d[,i]<-d[,i]-x[i]}

d                            # 将 d 中每个元素减去该列的平均，此为调整后的观测值

 v<-as.vector(d)    # 将d转化为向量

r<-rank(v)

dr<-matrix(r,nrow=3,ncol=5)   # Hodges-Lehmmann 秩数据表