- •Обозначения
- •Классификация фактов (случаев)
- •Функции классификации
- •Расстояние Махалонобиса
- •Пошаговый дискриминантный анализ
- •Пошаговый анализ с включением
- •Пошаговый анализ с исключением
- •Разработка статистической модели данных
- •Напишем макрос, который заполнит таблицу с векторами
- •Подвергнем данные дискриминантному анализу
- •Теоретический расчет ошибок классификации для модели
Разработка статистической модели данных
Модель представляет собой набор многомерных векторов , имеющих заданные вектора математических ожиданий и заданные ковариационные матрицы (одинаковые по классам). Компоненты векторов имеют нормальное распределение. Количество классов M= 2;3
Вариант №10 :
-
Размерность = 2
-
Объем выборки = 250
-
-
-
σ1 = 1.0 σ2 = 1.2
-
М = 2
-
Напишем макрос, который заполнит таблицу с векторами
Option Base 1
Sub Main
Dim n As Long
Dim i As Long
Dim m As Long
Dim St As Double
n = 250 'число наблюдений
m = 3 'число переменных
Dim s As New Spreadsheet
s.Header.Value = "Разработка статистической модели данных"
s.SetSize(n, m)
s.VariableName(1) = "Вектор 1"
s.VariableName(2) = "Вектор 2"
s.VariableName(3) = "Номер класса"
For i=1 To 100
s.Value(i, 1) = VNormal(Rnd(1), 1.5, 1.0)
s.Value(i, 2) = VNormal(Rnd(1), 1.0, 1.0)
s.Value(i, 3) = 1
Next i
For i = 101 To n
s.Value(i, 1) = VNormal(Rnd(1), 3.0, 1.2)
s.Value(i, 2) = VNormal(Rnd(1), 2.7, 1.2)
s.Value(i, 3) = 2
Next i
s.Visible = True
End Sub
Результат выполнения макроса (Один из результатов реализации)
|
Вектор 1 |
Вектор 2 |
Номер класса |
1 |
-0,150339054 |
-0,622227117 |
1 |
2 |
1,87506845 |
0,741934406 |
1 |
3 |
1,53765748 |
0,788950051 |
1 |
4 |
1,22441489 |
0,30858198 |
1 |
5 |
3,16478303 |
-0,203254201 |
1 |
6 |
0,149660006 |
0,259164289 |
1 |
7 |
2,53789502 |
0,900317028 |
1 |
8 |
1,24126184 |
0,116137968 |
1 |
9 |
3,02660053 |
1,67308541 |
1 |
10 |
0,971518277 |
1,29290362 |
1 |
11 |
0,470650132 |
1,79399462 |
1 |
12 |
1,33101628 |
0,232741505 |
1 |
13 |
1,91156835 |
1,39066945 |
1 |
14 |
0,572141902 |
1,58183723 |
1 |
15 |
2,7487247 |
1,82476327 |
1 |
16 |
-0,0949873617 |
0,981848025 |
1 |
17 |
-0,663212298 |
1,78791701 |
1 |
18 |
0,515379784 |
1,23464946 |
1 |
19 |
3,52067594 |
0,217110849 |
1 |
20 |
1,01363384 |
1,11369761 |
1 |
21 |
0,360868882 |
2,39373817 |
1 |
22 |
1,59753365 |
0,416170537 |
1 |
23 |
1,14014364 |
-1,32070685 |
1 |
24 |
1,63883867 |
1,85747626 |
1 |
25 |
1,94804129 |
0,772435936 |
1 |
26 |
1,50728203 |
0,645480827 |
1 |
27 |
1,99795086 |
0,47076958 |
1 |
28 |
0,0902818661 |
0,246969169 |
1 |
29 |
1,27983228 |
1,02907286 |
1 |
30 |
0,783814568 |
1,62105383 |
1 |
31 |
1,83803705 |
0,337541667 |
1 |
32 |
2,4267437 |
-0,138715106 |
1 |
33 |
0,718730948 |
0,665417966 |
1 |
34 |
1,54420909 |
1,54278558 |
1 |
35 |
2,02663518 |
0,860161009 |
1 |
36 |
-0,00319889246 |
0,607466284 |
1 |
37 |
1,1262229 |
0,696518577 |
1 |
38 |
2,29239108 |
0,739156613 |
1 |
39 |
0,413573446 |
-0,814077154 |
1 |
40 |
1,2676794 |
1,29145122 |
1 |
41 |
0,00829016592 |
0,25781585 |
1 |
42 |
3,40624041 |
2,54387721 |
1 |
43 |
2,13478293 |
2,15447339 |
1 |
44 |
2,40497064 |
2,08754914 |
1 |
45 |
1,66837993 |
1,60395899 |
1 |
46 |
2,16224823 |
0,107077806 |
1 |
47 |
0,908672887 |
1,66257009 |
1 |
48 |
-0,912178683 |
2,03105223 |
1 |
49 |
2,81182879 |
0,882504238 |
1 |
50 |
0,767826476 |
0,917573256 |
1 |
51 |
-0,633930523 |
0,682783081 |
1 |
52 |
1,84795646 |
1,47206785 |
1 |
53 |
1,98022805 |
1,90242359 |
1 |
54 |
3,39618616 |
-0,322999379 |
1 |
55 |
1,14029802 |
1,57699543 |
1 |
56 |
0,446273796 |
-0,545562693 |
1 |
57 |
0,375189681 |
0,589001887 |
1 |
58 |
1,74238805 |
1,36626652 |
1 |
59 |
2,58069728 |
-0,3465018 |
1 |
60 |
2,07134282 |
-0,0644481042 |
1 |
61 |
1,02880295 |
3,57112046 |
1 |
62 |
2,21638945 |
-1,27347471 |
1 |
63 |
2,71825552 |
1,87911756 |
1 |
64 |
0,540030691 |
-0,832514597 |
1 |
65 |
-0,0189660818 |
2,62190018 |
1 |
66 |
1,65832799 |
0,617062411 |
1 |
67 |
1,48130983 |
-0,555447794 |
1 |
68 |
2,20359075 |
0,348581938 |
1 |
69 |
1,53393988 |
1,5901975 |
1 |
70 |
1,32630065 |
2,0938983 |
1 |
71 |
0,537544584 |
1,28999986 |
1 |
72 |
0,953324093 |
0,864142511 |
1 |
73 |
1,77489738 |
2,03044713 |
1 |
74 |
2,39661646 |
0,375947563 |
1 |
75 |
2,71607319 |
-1,45876895 |
1 |
76 |
1,85603676 |
0,604248446 |
1 |
77 |
1,5683509 |
1,46593167 |
1 |
78 |
2,74805179 |
1,5144042 |
1 |
79 |
0,731502883 |
0,0718289501 |
1 |
80 |
2,1358069 |
2,79722097 |
1 |
81 |
0,570800959 |
-0,789400621 |
1 |
82 |
2,53005371 |
1,6114006 |
1 |
83 |
0,865644606 |
0,124472746 |
1 |
84 |
3,95156389 |
2,42480003 |
1 |
85 |
2,24892491 |
0,899908248 |
1 |
86 |
1,63328906 |
0,78086037 |
1 |
87 |
2,8174537 |
1,01989661 |
1 |
88 |
1,76379529 |
0,992755496 |
1 |
89 |
2,77019692 |
1,82729679 |
1 |
90 |
0,54500199 |
0,601071908 |
1 |
91 |
1,00187678 |
-0,807230494 |
1 |
92 |
0,925931223 |
2,4315203 |
1 |
93 |
2,26335779 |
1,63976599 |
1 |
94 |
0,306189608 |
1,0309593 |
1 |
95 |
3,38520932 |
1,21314095 |
1 |
96 |
1,23120334 |
0,120725589 |
1 |
97 |
-0,587943729 |
0,593088648 |
1 |
98 |
0,825809691 |
0,992194992 |
1 |
99 |
2,91856369 |
2,75445206 |
1 |
100 |
2,15198608 |
3,31280703 |
1 |
101 |
0,921620698 |
2,51730537 |
2 |
102 |
5,25056148 |
2,74610085 |
2 |
103 |
3,06074939 |
4,27534224 |
2 |
104 |
4,54643793 |
5,25989484 |
2 |
105 |
2,55156796 |
4,93073704 |
2 |
106 |
3,23502987 |
1,13631444 |
2 |
107 |
4,28389669 |
2,57794358 |
2 |
108 |
1,24123169 |
2,67338437 |
2 |
109 |
3,66599773 |
4,4215258 |
2 |
110 |
3,4008117 |
2,96835736 |
2 |
111 |
3,63917457 |
1,76729935 |
2 |
112 |
1,9299533 |
1,17907276 |
2 |
113 |
3,85431903 |
2,31554985 |
2 |
114 |
3,15328278 |
3,01569061 |
2 |
115 |
4,52440885 |
2,16124173 |
2 |
116 |
2,95339832 |
2,25708907 |
2 |
117 |
2,50792811 |
3,44012782 |
2 |
118 |
1,1087236 |
2,51896115 |
2 |
119 |
0,852300736 |
2,8236525 |
2 |
120 |
1,29147757 |
4,2081061 |
2 |
121 |
2,17505164 |
3,09913014 |
2 |
122 |
3,21950699 |
3,65032983 |
2 |
123 |
2,95750696 |
3,77657365 |
2 |
124 |
4,83081449 |
2,72846394 |
2 |
125 |
4,24087611 |
2,68742529 |
2 |
126 |
0,868765395 |
2,71798013 |
2 |
127 |
2,67847785 |
4,05953506 |
2 |
128 |
2,15643678 |
5,30148297 |
2 |
129 |
5,1804419 |
3,65313313 |
2 |
130 |
2,00280231 |
1,11213979 |
2 |
131 |
2,30043793 |
1,44342341 |
2 |
132 |
2,71033185 |
1,4533575 |
2 |
133 |
2,97023854 |
1,25445665 |
2 |
134 |
2,45322711 |
2,68012175 |
2 |
135 |
3,88398079 |
0,528093877 |
2 |
136 |
2,26153486 |
1,16342409 |
2 |
137 |
4,21936373 |
4,3842643 |
2 |
138 |
4,49876312 |
2,05166625 |
2 |
139 |
4,18587804 |
3,5965504 |
2 |
140 |
2,81735457 |
3,77672186 |
2 |
141 |
4,26510996 |
2,96839447 |
2 |
142 |
1,67755961 |
0,947529234 |
2 |
143 |
3,58166375 |
2,16587904 |
2 |
144 |
2,56670641 |
2,83325165 |
2 |
145 |
3,97630012 |
2,7850313 |
2 |
146 |
3,44251468 |
3,21791179 |
2 |
147 |
4,14755574 |
2,13874185 |
2 |
148 |
0,189899551 |
3,18081875 |
2 |
149 |
1,75038335 |
2,47867577 |
2 |
150 |
3,43854676 |
3,8472104 |
2 |
151 |
5,14706181 |
3,52243115 |
2 |
152 |
3,03507529 |
2,80395722 |
2 |
153 |
1,98528036 |
2,24186146 |
2 |
154 |
2,53803749 |
2,49437459 |
2 |
155 |
3,42421636 |
4,03935267 |
2 |
156 |
1,80691942 |
1,43525874 |
2 |
157 |
4,36863309 |
1,85696642 |
2 |
158 |
1,45113802 |
4,86716957 |
2 |
159 |
3,43217742 |
3,62238413 |
2 |
160 |
2,38341066 |
4,16349869 |
2 |
161 |
3,28206064 |
3,34516458 |
2 |
162 |
4,48183837 |
5,16883985 |
2 |
163 |
3,88173011 |
2,3117759 |
2 |
164 |
3,62773256 |
4,13968357 |
2 |
165 |
2,85177498 |
2,94906138 |
2 |
166 |
3,79761793 |
1,9755539 |
2 |
167 |
3,68532775 |
3,95648845 |
2 |
168 |
3,14226422 |
0,989412939 |
2 |
169 |
1,73224895 |
3,62234598 |
2 |
170 |
2,44422419 |
2,26949372 |
2 |
171 |
3,9818284 |
1,19767225 |
2 |
172 |
3,58369121 |
1,74598152 |
2 |
173 |
4,57298248 |
3,46182778 |
2 |
174 |
2,29378768 |
2,94038864 |
2 |
175 |
2,65730303 |
4,4191122 |
2 |
176 |
2,32724389 |
3,49559291 |
2 |
177 |
3,56591641 |
0,138576175 |
2 |
178 |
5,4536175 |
2,05485141 |
2 |
179 |
2,51619952 |
3,48452454 |
2 |
180 |
3,89800034 |
2,55245571 |
2 |
181 |
1,58582997 |
0,68701053 |
2 |
182 |
2,8664371 |
2,93776788 |
2 |
183 |
3,20112221 |
2,47244585 |
2 |
184 |
3,3447498 |
0,722568614 |
2 |
185 |
2,59264731 |
3,55848601 |
2 |
186 |
3,50231735 |
0,38777725 |
2 |
187 |
2,72032569 |
1,14319053 |
2 |
188 |
1,78529608 |
4,83236237 |
2 |
189 |
3,12370028 |
1,46605891 |
2 |
190 |
4,2264267 |
1,72601906 |
2 |
191 |
2,84668455 |
2,88845558 |
2 |
192 |
5,12225409 |
0,574533848 |
2 |
193 |
5,03936733 |
1,67629485 |
2 |
194 |
5,23178119 |
1,89350044 |
2 |
195 |
2,13318278 |
0,901133078 |
2 |
196 |
2,39574873 |
1,32053501 |
2 |
197 |
2,36821835 |
2,92956088 |
2 |
198 |
5,28149617 |
2,58814012 |
2 |
199 |
4,34114091 |
2,18518079 |
2 |
200 |
3,1524441 |
2,95547501 |
2 |
201 |
3,49970217 |
-0,218792001 |
2 |
202 |
2,79528485 |
1,04637728 |
2 |
203 |
1,37351915 |
1,74245726 |
2 |
204 |
0,721826972 |
0,0459476254 |
2 |
205 |
2,58580372 |
4,01046889 |
2 |
206 |
2,60804381 |
4,91280567 |
2 |
207 |
3,5425961 |
2,86337197 |
2 |
208 |
3,29726032 |
4,63765288 |
2 |
209 |
3,20649743 |
1,75920604 |
2 |
210 |
2,89301259 |
3,04941901 |
2 |
211 |
3,59866968 |
2,37903056 |
2 |
212 |
4,88198305 |
1,15466895 |
2 |
213 |
1,08496674 |
4,97157684 |
2 |
214 |
1,16020051 |
2,64105407 |
2 |
215 |
3,38932385 |
3,33931226 |
2 |
216 |
3,57925387 |
3,48592884 |
2 |
217 |
2,26992689 |
1,63162986 |
2 |
218 |
2,29309301 |
4,00569384 |
2 |
219 |
2,22244466 |
0,897335711 |
2 |
220 |
3,55492442 |
3,71227889 |
2 |
221 |
3,3636163 |
1,9035168 |
2 |
222 |
3,78564703 |
2,49300945 |
2 |
223 |
3,28730048 |
3,33131082 |
2 |
224 |
3,96971987 |
2,01673246 |
2 |
225 |
3,94020885 |
2,04818141 |
2 |
226 |
2,16671952 |
2,69641988 |
2 |
227 |
4,07871494 |
0,624979999 |
2 |
228 |
2,97537048 |
3,82129398 |
2 |
229 |
3,45555726 |
2,30758026 |
2 |
230 |
2,54709619 |
2,16334117 |
2 |
231 |
0,793236568 |
2,28078795 |
2 |
232 |
2,23553049 |
2,01467342 |
2 |
233 |
1,33951715 |
4,44166968 |
2 |
234 |
5,0931045 |
3,01695479 |
2 |
235 |
4,41863315 |
4,08422495 |
2 |
236 |
4,61697845 |
1,50656905 |
2 |
237 |
2,87885801 |
2,49113321 |
2 |
238 |
2,74404918 |
1,12815437 |
2 |
239 |
3,54411261 |
3,07736147 |
2 |
240 |
3,3550049 |
3,68395511 |
2 |
241 |
2,62232485 |
3,53521763 |
2 |
242 |
3,13759002 |
2,98945313 |
2 |
243 |
1,87801025 |
2,50498641 |
2 |
244 |
2,84515371 |
1,75464622 |
2 |
245 |
4,24258439 |
2,54476082 |
2 |
246 |
4,28576655 |
5,1630757 |
2 |
247 |
3,89285755 |
1,57402782 |
2 |
248 |
2,98899519 |
2,66440347 |
2 |
249 |
3,66563604 |
2,60526013 |
2 |
250 |
2,51282119 |
3,10431443 |
2 |