Добавил:

Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.

Вуз:

Самарский национальный исследовательский университет им. ак. С.П. Королёва (бывш. СГАУ, СамГУ)

Предмет:

[НЕСОРТИРОВАННОЕ]

Файл:

4ПМ_Б_Комп_Модел_Власова / Рекомендации студенту / Пример отчета.doc

Скачиваний:

Добавлен:

28.03.2016

Размер:

886.27 Кб

Скачать

☆

<<< < Предыдущая 1 23 / 33

Методы поиска общностей hj2 .

а) Метод Барта.

По каждому столбцу матрицы R найдем средние значения коэффициента корреляции (r_j)_ср. Если( r_j)_ср.сравнительно велико, то за общность принимается значение, которое несколько выше наибольшего в столбце коэффициента корреляции.

Если( r_j)_ср.сравнительно мало, то общность будет несколько меньше наибольшего в столбце коэффициента корреляции.

Средние (r_j)_ср.

0,195

0,068

0,328

0,334

-0,017

0,281

0,386

0,285

0,324

0,129

h_j² = 0,334

b). Метод триад.

Общности для каждого столбца вычисляются по формуле:

h_j² = r_ik*r_il/ r_{kl
,}

где r_ik и r_il коэффициенты корреляции, наибольшие в столбце.

h_j² =

0,700

h_j²=0.700

с) Метод малого центроида.

Для каждой переменной строится корреляционная матрица ( 44), включая саму переменную. В эту матрицу записывают оценки корреляции трех других переменных особенно тесно связанных с первой. По данным малой матрицы (44) рассчитываются общности по формуле

h_j² = ( r_j1)² /  r_ij,

где  r_j1- cумма элементов j-го столбца,  r_ij - сумма всех элементов матрицы.

h_j²=

0,4793

0,4746

0,6839

0,6267

0,5237

0,5455

0,6001

0,6700

0,6817

0,7807

h₁²:

0,700

0,427

0,356

0,332

0,427

0,7

-0,009

0,511

0,356

-0,009

0,7

0,08

0.332

0,511

0,08

0,7

h₂²:	0,7	0,494	-0,322	-0,471
	0,494	0,7	-0,416	-0,492
	-0,322	-0,416	0,7	0,563
	-0,471	-0,492	0,563	0,7

и т.д.

Первая итерация.

0-ой шаг. R_h Таблица 7

1	2	3	4	5	6	7	8	9	10	S_i⁽⁰⁾=	L_i⁽⁰⁾= S_i⁽⁰⁾/S_ma
0,700	0,210	0,427	0,005	0,356	0,006	0,332	0,002	0,050	-0,137	1,951	0,505
0,210	0,700	0,239	0,046	0,494	-0,208	0,242	-0,322	-0,247	-0,471	0,684	0,177
0,427	0,239	0,700	0,431	-0,009	0,398	0,511	0,287	0,324	-0,031	3,276	0,848
0,005	0,046	0,431	0,700	-0,194	0,618	0,545	0,458	0,530	0,200	3,341	0,865
0,356	0,494	-0,009	-0,194	0,700	-0,344	0,080	-0,416	-0,344	-0,492	-0,169	-0,044
0,006	-0,208	0,398	0,618	-0,344	0,700	0,370	0,494	0,519	0,254	2,808	0,727
0,332	0,242	0,511	0,545	0,080	0,370	0,700	0,406	0,500	0,175	3,861	1,000
0,002	-0,322	0,287	0,458	-0,416	0,494	0,406	0,700	0,681	0,563	2,854	0,739
0,050	-0,247	0,324	0,530	-0,344	0,519	0,500	0,681	0,700	0,526	3,240	0,839
-0,137	-0,471	-0,031	0,200	-0,492	0,254	0,175	0,563	0,526	0,700	1,288	0,334

1-ый шаг R_h² Таблица 8

S_i⁽¹⁾=

P_i⁽¹⁾=

R_h*S⁽⁰

L_i⁽¹⁾=

P_i⁽¹⁾/S_max

0,975

0,703

0,840

0,317

0,671

0,131

0,769

0,007

0,135

-0,391

4,157

0,375

0,703

1,326

0,323

-0,296

1,296

-0,654

0,182

-1,003

-0,857

-1,250

-0,230

-0,021

0,840

0,323

1,523

1,440

-0,138

1,271

1,569

1,134

1,289

0,399

9,650

0,870

0,317

-0,296

1,440

1,927

-0,891

1,849

1,697

1,832

1,949

1,129

10,953

0,987

0,671

1,296

-0,138

-0,891

1,556

-1,185

-0,316

-1,481

-1,372

-1,497

-3,357

-0,302

0,131

-0,654

1,271

1,849

-1,185

1,908

1,488

1,946

2,009

1,350

10,113

0,911

0,769

0,182

1,569

1,697

-0,316

1,488

1,806

1,477

1,645

0,726

11,043

0,995

0,007

-1,003

1,134

1,832

-1,481

1,946

1,477

2,249

2,268

1,782

10,210

0,920

0,135

-0,857

1,289

1,949

-1,372

2,009

1,645

2,268

2,318

1,714

11,098

1,000

-0,391

-1,250

0,399

1,129

-1,497

1,350

0,726

1,782

1,714

1,702

5,664

0,510

Вычисляем разности d.

d₁=L_i⁽¹⁾-L_i⁽⁰⁾
-0,1308
-0,1980
0,0211
0,1216
-0,2587
0,1841
-0,0050

Вновь возводим в квадрат матрицу R.

2-ой шаг R_h⁴ Таблица 9

Si(2)=S

Pi(2)=

Rh2*S(1)

Li(2)=

Pi(2)/Smax

3,48

3,09

4,09

2,71

2,20

1,68

4,06

0,86

1,50

-1,26

22,40

0,16

3,09

7,89

-2,38

-7,25

9,68

-8,79

-4,06

-10,85

-10,40

-10,21

-33,30

-0,24

4,09

-2,38

12,41

15,31

-6,77

14,75

14,39

15,22

16,17

9,47

92,67

0,68

2,71

-7,25

15,31

21,50

-13,28

21,74

18,55

23,44

24,25

16,59

123,55

0,91

2,20

9,68

-6,77

-13,28

13,18

-14,82

-9,33

-17,27

-17,08

-14,64

-68,13

-0,50

1,68

-8,79

14,75

21,74

-14,82

22,38

18,17

24,49

25,11

18,00

122,73

0,90

4,06

-4,06

14,39

18,55

-9,33

18,17

16,95

19,08

20,08

12,49

110,38

0,81

0,86

-10,85

15,22

23,44

-17,27

24,49

19,08

27,18

27,68

20,62

130,45

0,96

1,50

-10,40

16,17

24,25

-17,08

25,11

20,08

27,68

28,30

20,63

136,23

1,00

-1,26

-10,21

9,47

16,59

-14,64

18,00

12,49

20,62

20,63

16,75

88,43

0,65

Вновь вычисляем разноси. Нетрудно заметить, что они уже уменьшились.

d₂=L_i⁽²⁾-L_i⁽¹⁾
-0,210
-0,224
-0,189
-0,080
-0,198
-0,010
-0,185
0,038
0,000
0,139

Повторим все вычисления:

3-ий шаг R_h⁸ Таблица 10

Si(3)=S

Pi(3)=

Rh4*S(2)

Li(3)=

Si(3)/Smax

74,4

-16,2

192,7

225,9

-81,8

212,5

220,1

214,5

231,0

123,9

1396,9

0,068

-16,2

647,6

-830,5

-1312,1

1005,8

-1381,6

-1050,6

-1542,4

-1564,5

-1187,2

-7231,7

-0,351

192,7

-830,5

1464,5

2131,3

-1421,3

2185,3

1797,7

2386,8

2452,3

1741,2

12099,9

0,587

225,9

-1312,1

2131,3

3162,2

-2185,4

3263,6

2634,4

3584,1

3670,5

2651,0

17825,5

0,864

-81,8

1005,8

-1421,3

-2185,4

1605,3

-2281,8

-1779,9

-2530,0

-2576,6

-1915,8

-12161,6

-0,590

212,5

-1381,6

2185,3

3263,6

-2281,8

3375,6

2707,6

3713,7

3799,3

2759,2

18353,5

0,890

220,1

-1050,6

1797,7

2634,4

-1779,9

2707,6

2212,1

2963,1

3040,8

2172,6

14917,7

0,723

214,5

-1542,4

2386,8

3584,1

-2530,0

3713,7

2963,1

4091,9

4182,5

3051,4

20115,5

0,975

231,0

-1564,5

2452,3

3670,5

-2576,6

3799,3

3040,8

4182,5

4277,3

3112,2

20624,9

1,000

123,9

-1187,2

1741,2

2651,0

-1915,8

2759,2

2172,6

3051,4

3112,2

2296,2

14804,6

0,718

d₃=L_i⁽³⁾-L_i⁽²⁾
-0,097
-0,106
-0,094
-0,043
-0,090
-0,011
-0,087
0,018
0,000
0,069

4-ый шаг R_h¹⁶ Таблица 11

S_i⁽⁴⁾=

P_i⁽⁴⁾=

R_h⁸*S⁽³⁾

L_i⁽⁴⁾=

P_i⁽⁴⁾/S_max

3,E+05

-2,E+06

3,E+06

5,E+06

-3,E+06

5,E+06

4,E+06

5,E+06

4,E+06

3,E+07

3,E+05

-2,E+06

1,E+07

-2,E+07

-3,E+07

2,E+07

-3,E+07

-4,E+07

-3,E+07

-2,E+08

-2,E+06

3,E+06

-2,E+07

3,E+07

5,E+07

-3,E+07

5,E+07

4,E+07

5,E+07

6,E+07

4,E+07

3,E+08

3,E+06

5,E+06

-3,E+07

5,E+07

7,E+07

-5,E+07

7,E+07

6,E+07

8,E+07

6,E+07

4,E+08

5,E+06

-3,E+06

2,E+07

-3,E+07

-5,E+07

4,E+07

-5,E+07

-4,E+07

-6,E+07

-4,E+07

-3,E+08

-3,E+06

5,E+06

-3,E+07

5,E+07

7,E+07

-5,E+07

8,E+07

6,E+07

8,E+07

9,E+07

6,E+07

4,E+08

5,E+06

4,E+06

-3,E+07

4,E+07

6,E+07

-4,E+07

6,E+07

5,E+07

7,E+07

5,E+07

3,E+08

4,E+06

5,E+06

-3,E+07

5,E+07

8,E+07

-6,E+07

8,E+07

7,E+07

9,E+07

1,E+08

7,E+07

5,E+08

5,E+06

-4,E+07

6,E+07

8,E+07

-6,E+07

9,E+07

7,E+07

1,E+08

7,E+07

5,E+08

5,E+06

4,E+06

-3,E+07

4,E+07

6,E+07

-4,E+07

6,E+07

5,E+07

7,E+07

5,E+07

3,E+08

4,E+06

d₄=L_i⁽⁴⁾-L_i⁽³⁾

-0,0133

-0,0147

-0,0129

-0,0059

-0,0123

-0,0016

-0,0120

0,0024

0,0000

0,0095

Вычисляем разности:

Вычисления выполнены с точностью до d₄≤ 0,02, повторим вычисления еще раз до достижения точности.

5-ый шаг R_h³² Таблица 12

S_i⁽⁵⁾=

P_i⁽⁵⁾=

R_h¹⁶*S⁽⁴⁾

L_i⁽⁵⁾=

S_i⁽⁵⁾/S_max

2,E+14

-1,E+15

2,E+15

-2,E+15

2,E+15

3,E+15

2,E+15

1,E+16

5,E-02

-1,E+15

7,E+15

-1,E+16

-2,E+16

1,E+16

-2,E+16

-1,E+16

-2,E+16

-1,E+16

-9,E+16

-4,E-01

2,E+15

-1,E+16

2,E+16

3,E+16

-2,E+16

3,E+16

2,E+16

3,E+16

2,E+16

1,E+17

6,E-01

2,E+15

-2,E+16

3,E+16

4,E+16

-3,E+16

4,E+16

3,E+16

4,E+16

3,E+16

2,E+17

9,E-01

-2,E+15

1,E+16

-2,E+16

-3,E+16

2,E+16

-3,E+16

-2,E+16

-3,E+16

-2,E+16

-1,E+17

-6,E-01

2,E+15

-2,E+16

3,E+16

4,E+16

-3,E+16

4,E+16

3,E+16

4,E+16

5,E+16

3,E+16

2,E+17

9,E-01

2,E+15

-1,E+16

2,E+16

3,E+16

-2,E+16

3,E+16

4,E+16

3,E+16

2,E+17

7,E-01

3,E+15

-2,E+16

3,E+16

4,E+16

-3,E+16

4,E+16

5,E+16

4,E+16

2,E+17

1,E+00

3,E+15

-2,E+16

3,E+16

4,E+16

-3,E+16

5,E+16

4,E+16

5,E+16

4,E+16

2,E+17

1,E+00

2,E+15

-1,E+16

2,E+16

3,E+16

-2,E+16

3,E+16

4,E+16

3,E+16

2,E+17

7,E-01

Вычислим разности:

d₅=L_i⁽⁵⁾-L_i⁽⁴⁾

-0,0002

-0,0001

-0,0002

0,0000

-0,0002

0,0000

0,0001

Вычисления выполнены с точностью до d5≤0,0002

Вычисление нагрузок 1-го главного фактора.

Факторы	L_i⁽⁵⁾=U_i	U_i²	_R_h* L_i⁽⁵⁾	A₁ =U_i* ^U_i²)^1/2
х1	0,0542	0,0029	0,1902	130,0439
х2	-0,3655	0,1336	-1,2818	-0,2961
х3	0,5735	0,3289	2,0111	0,4646
х4	0,8582	0,7366	3,0095	0,6952
х5	-0,6022	0,3626	-2,1116	-0,4878
х6	0,8883	0,7890	3,1148	0,7195
х7	0,7111	0,5057	2,4936	0,5760
х8	0,9778	0,9561	3,4287	0,7920
х9	1,0000	1,0000	3,5066	0,8100
х10	0,7275	0,5292	2,5509	0,5892
 max{ _i }=			3,5066	Первый столбец матрицы факторных нагрузок
^{}			1,8726
Сумма (U_i²)=		5,3446
(U_i²)^1/2=		2,3118

Вычисляем воспроизведенную матрицу парных корреляций R_h⁺

0,0439

-0,2961

0,4646

0,6952

R_h⁺=

-0,4678

* 0,0439 ; -0, 2961; 0,4646; 0,6952; -0,4678; 0,7195; 0,5760; 0,7920; 0,8100; 0,5892 =

0,7195

0,5760

0,7920

0,8100

0,5892

Матрица R_h⁺

	1	2	3	4	5	6	7	8	9	10
1	0,0019	-0,0130	0,0204	0,0305	-0,0214	0,0316	0,0253	0,0348	0,0356	0,0259
2	-0,0130	0,0877	-0,1375	-0,2058	0,1444	-0,2130	-0,1705	-0,2345	-0,2398	-0,1745
3	0,0204	-0,1375	0,2158	0,3229	-0,2266	0,3342	0,2676	0,3679	0,3763	0,2737
4	0,0305	-0,2058	0,3229	0,4833	-0,3391	0,5002	0,4004	0,5506	0,5631	0,4096
5	-0,0214	0,1444	-0,2266	-0,3391	0,2379	-0,3509	-0,2809	-0,3863	-0,3951	-0,2874
6	0,0316	-0,2130	0,3342	0,5002	-0,3509	0,5177	0,4144	0,5698	0,5828	0,4240
7	0,0253	-0,1705	0,2676	0,4004	-0,2809	0,4144	0,3318	0,4562	0,4666	0,3394
8	0,0348	-0,2345	0,3679	0,5506	-0,3863	0,5698	0,4562	0,6273	0,6415	0,4667
9	0,0356	-0,2398	0,3763	0,5631	-0,3951	0,5828	0,4666	0,6415	0,6561	0,4773
10	0,0259	-0,1745	0,2737	0,4096	-0,2874	0,4240	0,3394	0,4667	0,4773	0,3472

Вычисляем остаточную матрицу R₁ по формуле R₁= R_h- R_h⁺ Матрица R₁

	1	2	3	4	5	6	7	8	9	10
1	0,6981	0,2233	0,4069	-0,0253	0,3771	-0,0256	0,3066	-0,0328	0,0148	-0,1631
2	0,2233	0,6123	0,3764	0,2522	0,3497	0,0046	0,4128	-0,0871	-0,0071	-0,2962
3	0,4069	0,3764	0,4842	0,1080	0,2173	0,0639	0,2434	-0,0814	-0,0524	-0,3049
4	-0,0253	0,2522	0,1080	0,2167	0,1451	0,1183	0,1451	-0,0921	-0,0326	-0,2096
5	0,3771	0,3497	0,2173	0,1451	0,4621	0,0070	0,3605	-0,0293	0,0511	-0,2042
6	-0,0256	0,0046	0,0639	0,1183	0,0070	0,1823	-0,0441	-0,0757	-0,0637	-0,1701
7	0,3066	0,4128	0,2434	0,1451	0,3605	-0,0441	0,3682	-0,0503	0,0331	-0,1640
8	-0,0328	-0,0871	-0,0814	-0,0921	-0,0293	-0,0757	-0,0503	0,0727	0,0398	0,0965
9	0,0148	-0,0071	-0,0524	-0,0326	0,0511	-0,0637	0,0331	0,0398	0,0439	0,0491
10	-0,1631	-0,2962	-0,3049	-0,2096	-0,2042	-0,1701	-0,1640	0,0965	0,0491	0,3528

Выполним аналогичные вычисления для остальных факторов. Всего 8 итераций.. Окончательно, решение принимает вид:

Матрица факторных нагрузок. Таблица 16.

Факторы	F1	F2	F3	F4	F5	F6	F7	F8	F9	F10
Х1	0,0439	0,6086	0,5172	0,2275	0,0651	-0,0638	-0,0614	-0,0593	0,0683	-0,1081
Х2	-0,2961	0,6755	-0,2313	-0,2395	-0,2085	-0,0003	-0,0003	-0,0003	0,0003	-0,0005
Х3	0,4646	0,5876	0,0041	0,2867	-0,2336	0,0452	0,0435	0,0420	-0,0484	0,0767
Х4	0,6952	0,2643	-0,3556	-0,0606	0,1081	-0,0695	-0,0669	-0,0646	0,0744	-0,1179
Х5	-0,4878	0,5706	0,1195	-0,2000	0,2761	0,0673	0,0647	0,0625	-0,0720	0,1140
Х6	0,7195	0,0641	-0,2658	0,2572	0,1946	-0,0305	-0,0294	-0,0284	0,0327	-0,0517
Х7	0,5760	0,5431	0,0548	-0,2517	0,0005	-0,0737	-0,0710	-0,0685	0,0789	-0,1250
Х8	0,7920	-0,1212	0,1414	-0,0862	-0,0256	0,1689	0,1625	0,1569	-0,1808	0,2863
Х8	0,8100	-0,0092	0,1076	-0,1442	0,0402	0,0836	0,0804	0,0777	-0,0895	0,1417
Х10	0,5892	-0,4423	0,2934	-0,2114	-0,0835	-0,1259	-0,1212	-0,1170	0,1348	-0,2134

Вчисление факторов и характерностей:

Ф А К Т О Р Ы											Общн.	Характ.
	F₁	F₂	F₃	F₄	F₅	F₆	F₇	F₈	F₉	F₁₀	h_j=Sa_ij²	d_j²=1-h_j
X₁	0,0439	0,6086	0,5172	0,2275	0,0651	-0,0638	-0,0614	-0,0593	0,0683	-0,1081	0,724	0,276
X₂	-0,2961	0,6755	-0,2313	-0,2395	-0,2085	-0,0003	-0,0003	-0,0003	0,0003	-0,0005	0,698	0,302
X₃	0,4646	0,5876	0,0041	0,2867	-0,2336	0,0452	0,0435	0,0420	-0,0484	0,0767	0,712	0,288
X₄	0,6952	0,2643	-0,3556	-0,0606	0,1081	-0,0695	-0,0669	-0,0646	0,0744	-0,1179	0,728	0,272
X₅	-0,4878	0,5706	0,1195	-0,2000	0,2761	0,0673	0,0647	0,0625	-0,0720	0,1140	0,725	0,275
X₆	0,7195	0,0641	-0,2658	0,2572	0,1946	-0,0305	-0,0294	-0,0284	0,0327	-0,0517	0,703	0,297
X₇	0,5760	0,5431	0,0548	-0,2517	0,0005	-0,0737	-0,0710	-0,0685	0,0789	-0,1250	0,730	0,270
X₈	0,7920	-0,1212	0,1414	-0,0862	-0,0256	0,1689	0,1625	0,1569	-0,1808	0,2863	0,864	0,136
X₉	0,8100	-0,0092	0,1076	-0,1442	0,0402	0,0836	0,0804	0,0777	-0,0895	0,1417	0,738	0,262
X₁₀	0,5892	-0,4423	0,2934	-0,2114	-0,0835	-0,1259	-0,1212	-0,1170	0,1348	-0,2134	0,788	0,212
_j= a_ij²	3,5066	2,0770	0,6531	0,4374	0,2373	0,0732	0,0678	0,0632	0,0839	0,2103	7,410	2,590

Таблица относительного Таблица относительного вклада Таблица относительного вклада факторов в суммарную общн. признаков в суммарную вклада характ. дисперсии

дисперсию дисперсию в суммарную дисперсию

фактор

Значение вклада

Вклад в %

Признак

Значение вклада

Вклад в%

Признак

Значение вклада

Вклад в%

_j= a_ij²

_j/a_ij²

h_i=a_ir²

h_i/a_ij²

1 - h_i

(1-h_i)/a_ij²

F₁

3,5066

47,32%

Х₁

0,724

9,77%

Х₁

0,276

10,67%

F₂

2,0770

28,03%

Х₂

0,698

9,42%

Х₂

0,302

11,65%

F₃

0,6531

8,81%

Х₃

0,712

9,61%

Х₃

0,288

11,13%

F₄

0,4374

5,90%

Х₄

0,728

9,82%

Х₄

0,272

10,51%

F₅

0,2373

3,20%

Х₅

0,725

9,78%

Х₅

0,275

10,62%

F₆

0,0732

0,99%

Х₆

0,703

9,48%

Х₆

0,297

11,47%

F₇

0,0638

0,86%

Х₇

0,730

9,85%

Х₇

0,270

10,42%

F₈

0,0632

0,85%

Х₈

0,864

11,66%

Х₈

0,136

5,24%

F₉

0,0839

1,13%

Х₉

0,738

9,96%

Х₉

0,262

10,13%

F₁₀

0,2103

2,84%

Х₁₀

0,788

10,64%

Х₁₀

0,212

8,17%

Проверка значимости факторов

a) Критерий илкса.

X²_Э= (n - (2*m+5)/6)*Ln(| _h⁺|).

Вычислим определитель матрицы R_h⁺для первого главного фактора.

Det ( R_h⁺)₁ = -7,374E-59

n = 8;

m = 10;

X²_Э= 5,131E+02 X²_{теорет.}= 43,8 k = m*(m-1)/2= 45 

Так как X²_Э > X²_{теорет.} , то первый главный фактор значимый.

Вычислим определитель матрицы R_h⁺для второго главного фактора.

Det ( R_h⁺)₂=	-1,107E-64
n = 8;
m = 10;

X²_Э= 5,645E+02 X²_{теорет.}= 43,8 k = m*(m-1)/2= 45 

Так как X²_Э >X²_{теорет.} , то второй главный фактор значимый.

Вычислим определитель матрицы R_h⁺для третьего главного фактора.

Det ( R_h⁺)₃=	-1,1E-160
n = 8;
m = 10;

X²_Э= 1,412E+03 X²_{теорет.}= 43,8 k = m*(m-1)/2= 45 

Так как X²_Э > X²_{теорет.} , то третий главный фактор значимый.

b) Критерий Лоули.

Запишем матрицу факторных нагрузок с тремя факторами

		Матрица А

	F₁	F₂	F₃
X₁	0,0439	0,6086	0,5172
X₂	-0,2961	0,6755	-0,2313
X₃	0,4646	0,5876	0,0041
X₄	0,6952	0,2643	-0,3556
X₅	-0,4878	0,5706	0,1195
X₆	0,7195	0,0641	-0,2658
X₇	0,5760	0,5431	0,0548
X₈	0,7920	-0,1212	0,1414
X₉	0,8100	-0,0092	0,1076
X₁₀	0,5892	-0,4423	0,2934

		Матрица А^Т

	X₁	X₂	X₃	X₄	X₅	X₆	X₇	X₈	X₉	X₁₀
F₁	0,0439	-0,2961	0,4646	0,6952	-0,4878	0,7195	0,5760	0,7920	0,8100	0,5892
F₂	0,6086	0,6755	0,5876	0,2643	0,5706	0,0641	0,5431	-0,1212	-0,0092	-0,4423
F₃	0,5172	-0,2313	0,0041	-0,3556	0,1195	-0,2658	0,0548	0,1414	0,1076	0,2934

		Матрица R_h= AА*^Т

	X₁	X₂	X₃	X₄	X₅	X₆	X₇	X₈	X₉	X₁₀
X₁	0,6399	0,2785	0,3802	0,0075	0,3877	-0,0669	0,3842	0,0341	0,0856	-0,0915
X₂	0,2785	0,5975	0,2585	0,0550	0,5022	-0,1083	0,1837	-0,3491	-0,2709	-0,5411
X₃	0,3802	0,2585	0,5611	0,4768	0,1092	0,3708	0,5869	0,2973	0,3713	0,0150
X₄	0,0075	0,0550	0,4768	0,6796	-0,2307	0,6116	0,5245	0,4683	0,5224	0,1884
X₅	0,3877	0,5022	0,1092	-0,2307	0,5778	-0,3462	0,0355	-0,4386	-0,3875	-0,5047
X₆	-0,0669	-0,1083	0,3708	0,6116	-0,3462	0,5924	0,4346	0,5245	0,5536	0,3176
X₇	0,3842	0,1837	0,5869	0,5245	0,0355	0,4346	0,6297	0,3981	0,4675	0,1153
X₈	0,0341	-0,3491	0,2973	0,4683	-0,4386	0,5245	0,3981	0,6619	0,6579	0,5618
X₉	0,0856	-0,2709	0,3713	0,5224	-0,3875	0,5536	0,4675	0,6579	0,6678	0,5129
X₁₀	-0,0915	-0,5411	0,0150	0,1884	-0,5047	0,3176	0,1153	0,5618	0,5129	0,6289

Вычислим определители матриц R_h⁺и R_h
.

Опр. ( R_h⁺) = 2,00E-113 X²_Э=(n-1)*Ln( |R_h⁺| / |R_h| )

Опр. ( R_h) = 2,742E-52

n = 8; число объектов k -число степеней свободы.

m = 10; число признаков уровень значимости. r - количество факторов.

X²_Э=9,854E+02 X²_{теорет.}= 28,9; k =((m-r)2-m-r)/2=18 

Так как X²_Э >X²_{теорет.} , то количество выбранных общих факторов не достаточно.

Вращение пространства общих факторов

Выберем три первых главных фактора. Вклад трех главных факторов в общую дисперсию составляет

Выполним вращение пространтва факторов на 60 градусов. Тогда матрицы примут вид :

f = 60^o

cos f

sin f

0,5

0,866025

T₁₂ =

- sin f

cos f

0 =

-0,8660

0,5

cos f

sin f

0,5

0,86603

T₁₃ =

0 =

- sin f

cos f

-0,86603

0,5

T₂₃ =

cos f

sin f =

0,5

0,86603

- sin f

cos f

-0,86603

0,5

Найдем произведение матриц T = T₁₂ *T₁₃ *T₂₃

0,2500

0,8660

0,4330

T₁₂*T₁₃ =

-0,4330

0,5000

-0,7500

-0,8660

0,0000

0,5000

0,2500

0,0580

0,9665

T₁₂ *

T₁₃ *T₂₃=

-0,4330

0,8995

0,0580

-0,8660

-0,4330

0,2500

	Результирующая матрица имеет вид

		Матрица W

	W₁	W₂	W₃
X₁	-0,7005	0,3261	0,2071
X₂	-0,1662	0,6906	-0,3048
X₃	-0,1419	0,5538	0,4841
X₄	0,3673	0,4321	0,5983
X₅	-0,4725	0,4332	-0,4085
X₆	0,3823	0,2144	0,6327
X₇	-0,1387	0,4982	0,6019
X₈	0,1281	-0,1243	0,7938
X₉	0,1133	-0,0079	0,8092
X₁₀	0,0847	-0,4907	0,6172

<<< < Предыдущая 1 23 / 33

Соседние файлы в папке Рекомендации студенту

#
28.03.201630.72 Кб11Методические рекомендации студенту.doc
#
28.03.2016179.71 Кб13Образец отчета по ДС_4ПМ.doc
#
28.03.2016886.27 Кб20Пример отчета.doc