用R可以实现。下面的程序里模拟产生了51组数据，每组5例。前50组的总体均数是相同的。只有最后一组即第51组的均数和它们不同。

用pairwise.t.test命令可以作两两比较。

times<-50
n<-5
y<-c(rnorm(times*n,0,1),rnorm(5,3,1))
grp<-rep(1:51,rep(5,51))
data<-cbind(y,grp)

oneway.test(y ~ grp, var.equal = TRUE)

pairwise.t.test(y, grp, p.adjust.method = "bonferroni")