1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
| #include <iostream>
#include <ctime>
using namespace std;
int main() {
setlocale(LC_ALL , "");
cout << "Критерий Манна-Уитни.";
//int x[6] = { 15, 23, 25, 26, 28, 29 };
//int y[8] = { 12, 14, 18, 20, 22, 24, 27, 30 };
/*7 8 9 10 11 12 13 14 15 16 17 18*/
int u[18][12] = { { 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7 }, //В данной таблице представлены критические значения
/*4*/ { 3, 4, 4, 5, 6, 7, 8, 9, 10, 11, 11, 12 }, //критерия U-Манна-Уитни для уровня значимости 0,05
/*5*/ { 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 17, 18 },
/*6*/ { 6, 8, 10, 11, 13, 14, 16, 17, 19, 21, 22, 24 },
/*7*/ { 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30 },
/*8*/ { 10, 13, 15, 17, 19, 22, 24, 26, 29, 31, 34, 36 },
/*9*/ { 12, 15, 17, 20, 23, 26, 28, 30, 34, 37, 39, 42 },
/*10*/ { 14, 17, 20, 23, 26, 29, 33, 36, 39, 42, 45, 48 },
/*11*/ { 16, 19, 23, 26, 30, 33, 37, 40, 44, 48, 51, 55 },
/*12*/ { 18, 22, 26, 29, 33, 37, 41, 45, 49, 53, 57, 61 },
/*13*/ { 20, 24, 28, 33, 37, 41, 45, 50, 54, 59, 63, 67 },
/*14*/ { 22, 26, 31, 36, 40, 45, 50, 55, 59, 64, 67, 74 },
/*15*/ { 24, 29, 34, 39, 44, 49, 54, 59, 64, 70, 75, 80 },
/*16*/ { 26, 31, 37, 42, 47, 53, 59, 64, 70, 75, 81, 86 },
/*17*/ { 28, 34, 39, 45, 51, 57, 63, 67, 75, 81, 87, 93 },
/*18*/ { 30, 36, 42, 48, 55, 61, 67, 74, 80, 86, 93, 99 },
/*19*/ { 32, 38, 45, 52, 58, 65, 72, 78, 85, 92, 99, 106 },
/*20*/ { 34, 41, 48, 55, 62, 69, 76, 83, 90, 98, 105, 112 }
};
int n1, n2;
cout << "\nВведите количество значений первой выборки n1: ";
cin >> n1;
cout << "\nВведите количество значений второй выборки n2: ";
cin >> n2;
int* x = new int[n1];
srand(time(0));
for (int i(0); i < n1; i++)
x[i] = rand() % 1 + 99;
int* y = new int[n1];
srand(time(0));
for (int i(0); i < n2; i++)
y[i] = rand() % 1 + 99;
int *ind1 = new int[n1];
int *ind2 = new int[n2];
int m(0), sum1(0), sum2(0), nx;
if (n1 > 2 && n2 > 6) {
int** xy = new int*[2];
for (int i(0); i < 2; i++)
xy[i] = new int[n1 + n2];
for (int i(0); i < n1 + n2; i++)
xy[0][i] = i + 1;
for (int i(0); i < n2; i++)
xy[1][i] = y[i];
for (int i(n2); i < n1 + n2; i++)
xy[1][i] = x[i - n2];
for (int i(1); i < n1 + n2; i++)
for (int j(i); j > 0 && xy[1][j - 1] > xy[1][j]; j--)
swap(xy[1][j], xy[1][j - 1]);
for (int i(0); i < n1 + n2; i++)
cout << xy[0][i] << " ";
for (int i(0); i < n1 + n2; i++)
for (int j(0); j < n1; j++)
if (xy[1][i] == x[j]) {
m++;
ind1[m - 1] = i + 1;
sum1 += ind1[m - 1];
}
m = 0;
for (int i(0); i < n1 + n2; i++)
for (int j(0); j < n2; j++)
if (xy[1][i] == y[j]) {
m++;
ind2[m - 1] = i + 1;
sum2 += ind2[m - 1];
}
if (sum2 > sum1) {
sum1 = sum2;
nx = n2;
}
else
nx = n1;
int U = n1*n2 + nx*(nx + 1) / 2 - sum1;
cout << "\nЗначение U-критерия Манна-Уитни: " << U;
cout << "\nКритическое значение критерия для данных n1 и n2: " << u[n1 - 3][n2 - 8];
if (U > u[n1 - 3][n2 - 8])
cout << "\nНулевая гипотеза об однородности двух выборок принимается.";
else
cout << "\nНулевая гипотеза об однородности двух выборок отвергается.";
for (int j(0); j < n1 + n2; j++)
delete[] xy[j];
delete[] xy;
}
else
cout << "\nНеверно заданы значения!";
delete[] x;
delete[] y;
system("pause");
return 0;
} |