restart:
with(LinearAlgebra):
with(stats):H:=Matrix(5, 3, [[186, 38, 35], [227, 54, 45], [219, 78, 78], [355, 112, 140], [653, 285, 259]]);
Quantil:=0.99;## H\344ufigkeiten Merkmalkombinationen
r, s:=Dimension(H);
HY:=Vector(r):
HZ:=Vector(s):
n:=0:
for i from 1 to r do
for j from 1 to s do
HY[i]:=HY[i] + H[i, j];
HZ[j]:=HZ[j] + H[i, j];
n:=n + H[i, j];
od:
od:
## Wahrscheinlichkeit Sch\344tzwert
P:=Matrix(5, 3):
for i from 1 to r do
for j from 1 to s do
P[i, j]:=HY[i] / n * HZ[j] / n;
od:
od:## Testgr\366sse Z
Z:=0:
for i from 1 to r do
for j from 1 to s do
Z:=Z + (H[i, j] - n * P[i, j])^2 / (n * P[i, j]);
od:
od:Z:=evalf(Z);printf("Siehe Tabelle Chi-Quadrat-Verteilung: alpha=%f, n=%d", Quantil, (r - 1) * (s - 1));
Q:=20.09;
## Tabelle von Maple verwenden
Q:=statevalf[icdf, chisquare[(r - 1) * (s - 1)]](Quantil);printf("Hypothese: Y[i] und Z[i] sind stochastisch unabh\344ngig\134n");
if (Z > Q) then
printf("--> Hypothese abgelehnt");
else
printf("--> Hypothese angenommen");
end if;TTdSMApJM1JUQUJMRV9TQVZFLzk5ODI3NlgsJSlhbnl0aGluZ0c2IjYiW2dsISIlISEhIzAiJiIkIiQnPSIkRiMiJD4jIiRiJCIkYCciI1EiI2EiCiN5IiQ3IiIkJkciI04iI1hGLiIkUyIiJGYjRiYK