# Kod 3.5. Testowanie hipotezy o dwuwymiarowej normalności
library(MASS)
library(ellipse)
liczba_cel=n=10
Sigma <- matrix(c(10,3,3,8),2,2)
z=mvrnorm(n=10, rep(0, 2), Sigma)
zsr=c(mean(z[,1]),mean(z[,2]))
S=cov(z)
tablica=matrix(0,nrow=10)
for (i in 1:liczba_cel) {
for (j in 1:n) {
if ((z[j,])%*%solve(Sigma)%*%(z[j,]) > qf((i-1)/10,2,8)*2*9/8) tablica[i]=tablica[i]+1
}
}
tablica
##Wyznaczenie liczby pustych cel
lpc=0
for (i in 2:10) {
if (tablica[i]==tablica[i-1]) lpc=lpc+1
}
lpc
Ka=4.9 # wartość krytyczna odczytana z tablic
if (lpc>4.9) print('Hipoteza H0 odrzucona') else print('Brak podstaw')