Exercício III - Obrigatório

 Exercício III - Obrigatório

Download do Relatorio do Exercício III

Clicar Aqui !!!

Colocar os dois ultims digitos do RG no lugar dos sinais de interrogação.

Regiao	Rb	Se	Cs	Br
Anazonia	20,2	60,7	32??	32
Anazonia	21,3	54,8	3100	28
Anazonia	19,3	49,6	2800	30
Anazonia	21,1	35,6	3300	24
Anazonia	20,475	51,85	3100	20
Cerrado	22,4	14,9	2600	33
Cerrado	21,9	17,8	2700	21
Cerrado	23,8	18,6	3200	25
Cerrado	24,1	15,1	3300	34
Cerrado	23,05	16,6	2950	27
Sul_1	27,3	2,5	2700	33
Sul_1	23,4	4,3	2300	27
Sul_1	25,2	2,3	2600	25
Sul_1	26,4	2,6	3200	25
Sul_1	25,575	2,925	2700	22
Sul_2	26,2	4,1	2600	23
Sul_2	24,2	2,1	2700	34
Sul_2	25,4	1,9	2650	35
Sul_2	24,5	2,3	2680	20
Sul_2	25,075	2,6	2657,5	21

MANOVA - Cluster e PCA - Random Forest (Machine Learning Supervisionado para Clasificação)

 MANOVA - Cluster e PCA - Random Forest (Machine Learning Supervisionado para Clasificação)

 

MANOVA 

Testes de comparaçoes multiplas de medias de Tukey é mais consevador que o teste de Duncan (usado na genetica).

Seria melhor que o Tukey o testes Teste de Dunn (Dunn's Test) ou comparaçoes multiplas de K Wallis que sao baseados mas medianas nao nas medias.

data imc_dat;

input cat $ imc corr kcal;

cards;

AT 20.2 60.7 3200

AT 21.3 54.8 3100

AT 19.3 49.6 2800

AT 21.1 52.3 3300

SEM 22.4 14.9 2600

SEM 21.9 17.8 2700

SEM 23.8 18.6 3200

SEM 24.1 15.1 3300

SE  27.3 2.5 2700

SE  23.4 4.3 2300

SE  25.2 2.3 2600

SE  26.4 2.6 3200

PR 26.2 4.1 2600

PR 24.2 2.1 2700

PR 25.4 1.9 2650

;

proc print;

run;

/* Programa de ANOVA fornecido pela IA - Gabriel Maciel */;

proc glm;

 class cat;

 model imc corr kcal  = cat;

 contrast " Atl e Semiat Vs Seden e Prof"  cat 1 -1 -1 1;

 contrast " Professor Vs Sedentario" cat  0 1 -1 0;

 contrast " Atleta Vs Semiatleta" cat 1 0 0 -1;

 contrast "Semiatleta Vs Sedentario e Professor" cat 0 1 1 -2;

 manova h=_all_ / printe printh;

 contrast " Atl e Semiat Vs Seden e Prof" 1 -1 -1 1;

 contrast " Professor Vs Sedentario" cat 0 1 -1 0;

 contrast " Atleta Vs Semiatleta" cat 1 0 0 -1;

  contrast "Semiatleta Vs Sedentario e Professor" cat 0 1 1 -2;

run;

/* 
contrast " Atl e Semiat Vs Seden e Prof"  cat 1 -1 -1 1;
Pega a ordem alfabetica dos niveis dos fatores para os coeficientes
   dos contrastes 
Se colocar o contraste depois do comando MANOVA faz os contrastes univariados 
tambem

*/

Correlação Residual MANOVA:

Partial Correlation Coefficients from the Error SSCP Matrix / Prob > |r|
DF = 11	imc	corr	kcal
imc	1.000000	0.017729 0.9564	0.697770 0.0116
corr	0.017729 0.9564	1.000000	0.132496 0.6814
kcal	0.697770 0.0116	0.132496 0.6814

Outro Exemplo de Contrastes, Ordem Alfabética:





- Análise Multivariada Exemplos: 

Cluster An. para Dados de MANOVA

data QV;

input categ $ IMC Corr_Caminh Kcal;

datalines;

AT 20.475 54.35 3100

PR 25.26666667 2.7 2650

SE 25.575 2.925 2700

SEM 23.05 16.6 2950

;

proc print; run;

/* input categ $ IMC Corr_Caminh Kcal; */;

proc cluster outtree = arvore method = average;

var IMC Corr_Caminh Kcal;

id categ;

run;

PROC TREE DATA = arvore;

RUN;

Cluster com Probabilidades de MANOVA

PCA - Biplot

Comparar resultados ANOVA com Kruskal-Wallis.

data imc_dat;

input cat $ imc corr kcal;

cards;

AT    20.2  60.7  3200

AT    21.3  54.8  3100

AT    19.3  49.6  2800

AT    21.1  52.3  3300

SEM   22.4  14.9  2600

SEM   21.9  17.8  2700

SEM   23.8  18.6  3200

SEM   24.1  15.1  3300

SE    27.3  2.5   2700

SE    23.4  4.3   2300

SE    25.2  2.3   2600

SE    26.4  2.6   3200

PR    26.2  4.1   2600

PR    24.2  2.1   2700

PR         25.4         1.9        2650

;

proc print;

run;

proc npar1way data=imc_dat wilcoxon dscf;

 class cat;

 var imc corr kcal = cat;

run;

Tem um erro nesse programa, onde esta? Compare com os slides apresentados antes do programa SAS (Dois slides de Kruskal Wallis com trechos de programa SAS).

 var imc corr kcal = cat; ==> ANOVA, GLM, MANOVA
 var imc corr kcal;          ==> Kruskal Wallis

Resultados de Kruskal Wallis

Arquivo Word para Download:

Download_K_Wallis

Resultados ANOVA dados brutos, transformados e da Estat. Robusta. Tese doutorado do Gabriel.