# Kod 3.14 Test permutacyjny dla porównania dwóch wartości oczekiwanych
x1=rnorm(5,10,1)
x2=rnorm(7,11,1)
n1=length(x1)
n2=length(x2)
x=c(x1,x2)
N=1000
T0=mean(x1)-mean(x2)
T=c()
for (i in 1:N)
{
xs=sample(x)
T=c(T,mean(xs[1:n1])-mean(xs[(n1+1):(n1+n2)]))
}
hist(T)
sum(abs(T)>=abs(T0))/N
