# Kod 3.3. Zastosowanie testu pustych cel dla weryfikacji hipotezy o jednostajnym dwuwymiarowym rozkładzie badanej zmiennej
n=50
x=iris[1:n,1]
y=iris[1:n,2]
m=matrix(NA,4,3)
for (i in 1:4)
for (j in 1:3) {
m[i,j]=sum(x>=(i/2+3.5) & x<(i/2+4) & y>=(j+1) & y<(j+2))
}
Kn=sum(m==0)
# Wyznaczenie p-wartości
Ks=c()
N=1000
for (sym in 1:N) {
xs=runif(n,4,8)
ys=runif(n,2,5)
ms=matrix(NA,4,3)
for (i in 1:4)
for (j in 1:3) {
ms[i,j]=sum(xs>=(i+3) & xs<(i+4) & ys>=(j+1) & ys<(j+2))
}
Ks=c(Ks,sum(ms==0))
}
p.val=sum(Ks>=Kn)/n
p.val