Kod 3.6 - NMSwBE

# Kod 3.6 Symulacyjne wyznaczenie rozkładu statystyki T= a

x=c(1,1,1,1,0,0,0,0)

y=c(0,0,0,0,1,1,1,1)

A=table(x,y)

T0= A[1,1]

T=c()

N=10000

for (i in 1:N) {

As=table(x,sample(y))

T=c(T,As[1,1])

}

sum(T>=T0)/N

plot(table(T),type='h',ylab='n')