Reunião 2/11/2024 Alvaro PCA e Cluster - Não PCA IMC buzcar outro

 

Reunião 2/11/2024

data imc_dat;

 

input cat $ imc corr kcal;

 

cards;

AT 81.9 12400 217.4

PR 75.8 7950 8.1

SE 102.3 10800 11.7

SEM 92.2 11800 66.4

;

 

proc print;

 

run;

 

title "PCA - Biplot - Com todas as Obsevações";

 

title "Melhor fazer com Medias ou Medianas ou Trimedias";

 

proc prinqual data=imc_dat plots=(MDPref)

 

                 /* project onto Prin1 and Prin2 */

 

              ; /* use COV scaling */

 

   transform identity(imc corr kcal);  /* identity transform */

 

   id cat;

 

   ods select MDPrefPlot;

 

run;

 

 

Robust ANOVA

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;

 

/* input cat $ imc corr kcal;

proc anova;

   class cat ;

   model imc corr kcal = cat;

   means tukey / lines;

run; */

Title "Robust ANOVA ou NPANOVA ou Kuskal Wallis (um fator)";

 

proc npar1way data=imc_dat wilcoxon dscf;

 

 class cat;

 

 var imc corr kcal;

 

run;

 

data  pessoas;

 

input cat $ imc corr kcal;

 

cards;

 

AT 20.5 54.4 3159

 

PR 25.3 2.7 2650

 

SE 25.6 2.9 2700

 

SEM 23.1 16.6 2959

 

;

 

proc cluster data=pessoas outtree = arvore method = average;

 

var imc corr;

 

id cat;

 

run;

 

PROC TREE DATA = arvore;

 

RUN;