Procedura testu permutacyjnego - porównanie dwóch populacji
# Wprowadzenie danych z prób x1<-c(0,4); x2<-c(3,5,7) x<-c(x1,x2); N=10000 # obliczenie wartości statystyki T0=mean(x1)-mean(x2) T=c() for (i in 1:N) { xs=sample(x) T[i]=mean(xs[1:2])-mean(xs[3:5]) } # empiryczny rozkład statystyki T hist(T)