Ex. III Obrigatorio - Winzor. - MANOVA - Cluster - PCA - Rede Neural

data rastreab;

input Registro Regiao $ Rb Se Cs Br;

datalines;

1 Anazonia 73.99 100.00 96.97 91.??

2 Anazonia 78.02 90.28 93.94 80.00

3 Anazonia 70.70 81.71 94.50 85.71

4 Anazonia 77.29 58.65 100.00 68.57

5 Anazonia 75.00 85.42 93.94 57.14

6 Cerrado 82.05 24.55 78.79 94.29

7 Cerrado 80.22 29.32 81.82 60.00

8 Cerrado 87.18 30.64 96.97 71.43

9 Cerrado 88.28 24.88 100.00 97.14

10 Cerrado 84.43 27.35 89.39 77.14

11 Sul_1 100.00 4.12 81.82 76.40

12 Sul_1 93.40 4.10 79.00 77.14

13 Sul_1 92.31 3.79 78.79 71.43

14 Sul_1 96.70 4.28 81.60 71.43

15 Sul_1 93.68 4.82 81.82 72.86

16 Sul_2 95.97 3.20 81.30 65.71

17 Sul_2 88.64 3.46 81.82 97.14

18 Sul_2 93.04 3.13 80.30 100.00

19 Sul_2 89.74 3.79 81.21 57.14

20 Sul_2 91.85 4.28 80.53 60.00

;

proc print; run;

/* input Registro Regiao $ Rb Rb Se Cs Br; */;

proc anova;

  class Regiao;

  model Rb Se Cs Br = Regiao;

  run;

  

Programa com Dados Winsorizados

data rastreab;

input Registro Regiao $ Rb Se Cs Br;

datalines;

1 Amazonia 73.99 100.00 96.97 91.59

2 Amazonia 78.02 90.28 93.94 80.00

3 Amazonia 70.70 81.71 94.50 85.71

4 Amazonia 77.29 81.8 100.00 68.57

5 Amazonia 75.00 85.42 93.94 57.14

6 Cerrado 82.05 24.55 78.79 94.29

7 Cerrado 80.22 29.32 81.82 60.00

8 Cerrado 87.18 30.64 96.97 71.43

9 Cerrado 88.28 24.88 100.00 97.14

10 Cerrado 84.43 27.35 89.39 77.14

11 Sul_1 100.00 4.12 81.82 76.40

12 Sul_1 93.40 4.10 79.00 77.14

13 Sul_1 92.31 3.79 78.79 71.43

14 Sul_1 96.70 4.28 81.60 71.43

15 Sul_1 93.68 4.82 81.82 72.86

16 Sul_2 95.97 3.20 81.30 65.71

17 Sul_2 88.64 3.46 81.82 97.14

18 Sul_2 93.04 3.13 80.30 100.00

19 Sul_2 89.74 3.79 81.21 57.14

20 Sul_2 91.85 4.28 80.53 60.00

;

proc print; run;

/* input Registro Regiao $ Rb Rb Se Cs Br; */;

proc anova;

  class Regiao;

  model Rb Se Cs Br = Regiao;

  run;

Dados Winsorizados

1	Anazonia	73,99	100	96,97	91,59
2	Anazonia	78,02	90,28	93,94	80
3	Anazonia	70,7	81,71	94,5	85,71
4	Anazonia	77,29	81,9	100	68,57
5	Anazonia	75	85,42	93,94	57,14
6	Cerrado	82,05	24,55	78,79	94,29
7	Cerrado	80,22	29,32	81,82	60
8	Cerrado	87,18	30,64	96,97	71,43
9	Cerrado	88,28	24,88	100	97,14
10	Cerrado	84,43	27,35	89,39	77,14
11	Sul_1	100	4,12	81,82	76,4
12	Sul_1	93,4	4,1	79	77,14
13	Sul_1	92,31	3,79	78,79	71,43
14	Sul_1	96,7	4,28	81,6	71,43
15	Sul_1	93,68	4,82	81,82	72,86
16	Sul_2	95,97	3,2	81,3	65,71
17	Sul_2	88,64	3,46	81,82	97,14
18	Sul_2	93,04	3,13	80,3	100
19	Sul_2	89,74	3,79	81,21	57,14
20	Sul_2	91,85	4,28	80,53	60

Programa SAS para IAI Não Supervisionada para Agrupamentos e Distancias Multivariados

data cluster;

input Regioes $ Rb Se Cs Br;

datalines;

Anazonia 75 87.862 95.87 76.602

Cerrado 84.432 27.348 89.394 80

Sul_1 95.218 4.222 80.606 73.852

Sul_2 91.848 3.572 81.032 75.998

;

proc print; run;

/* IAI ou ML Não Supervisioando para Agrupamentos e Distancias 

   Multivariados - Cluster Analysis */;

/* input Regioes $ Rb Se Cs Br; */;  

proc cluster outtree = arvore method = average;

var Rb Se Cs Br;

id Regioes;

run;

PROC TREE DATA = arvore;

RUN;  

Programa MANOVA

data rastr_soja;

input Registro Regiao $ Rb Se Cs Br;

datalines;

1 Anazonia 73.99 100 96.97 91.59

2 Anazonia 78.02 90.28 93.94 80

3 Anazonia 70.7 81.71 94.5 85.71

4 Anazonia 77.29 81.9 100 68.57

5 Anazonia 75 85.42 93.94 57.14

6 Cerrado 82.05 24.55 78.79 94.29

7 Cerrado 80.22 29.32 81.82 60

8 Cerrado 87.18 30.64 96.97 71.43

9 Cerrado 88.28 24.88 100 97.14

10 Cerrado 84.43 27.35 89.39 77.14

11 Sul_1 100 4.12 81.82 76.4

12 Sul_1 93.4 4.1 79 77.14

13 Sul_1 92.31 3.79 78.79 71.43

14 Sul_1 96.7 4.28 81.6 71.43

15 Sul_1 93.68 4.82 81.82 72.86

16 Sul_2 95.97 3.2 81.3 65.71

17 Sul_2 88.64 3.46 81.82 97.14

18 Sul_2 93.04 3.13 80.3 100

19 Sul_2 89.74 3.79 81.21 57.14

20 Sul_2 91.85 4.28 80.53 60

;

proc print; run;

/* input Registro Regiao $ Rb Se Cs Br; */

/* Ordem Alfabetica dos Niveis do Fator Regiao: [Amazonia Cerrado Sul_1 Sul_2] */

proc glm;

 class Regiao;

 model  Rb Se Cs Br = Regiao;

 contrast "Sul_1 Vs Sul_2" Regiao  0 0 1 -1;

 contrast " Cerrado Vs Sul_1 e Sul_2" Regiao 0 2 -1 -1;

 contrast " Amazonia Regia Vs Cerrado, Sul_1 e Sul_2" Regiao 3 -1 -1 -1; 

 manova h=_all_ / printe printh;

  contrast "Sul_1 Vs Sul_2" Regiao  0 0 1 -1;

  contrast " Cerrado Vs Sul_1 e Sul_2" Regiao 0 2 -1 -1;

  contrast " Amazonia Regia Vs Cerrado, Sul_1 e Sul_2" Regiao 3 -1 -1 -1;

run;

Arquivo para Download - Rede Neural

Clicar Aqui !!!