Hodges-Lehmann modified rank sum test (nonparametric statistics exercise)

以下表格展示了美国三家主要汽车制造商（分别标记为A、B、C）针对五款不同车型在某年度所生产车辆的燃油消耗情况，旨在分析各企业产品在油耗表现上是否存在显著差异。

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-Leh