PCA - Biplot - Com todas as Obsevações

 

PCA - Biplot - Com todas as Obsevações

Peguei o programa da Internet

Palavras chave: pca biplot sas

Achei esse site

https://blogs.sas.com/content/iml/2019/11/13/create-biplots-sas.html

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;

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;

data imc_dat;

input cat $ imc corr kcal;

cards;

AT 20.65 53.55 3150

SEM 23.1 16.45 2950

SE 25.8 2.55 2650

PR 25.4 2.1 2650

;

proc print;

run;

title "PCA - Biplot - Com Medianas";

/*data imc_dat;

input cat $ imc corr kcal; */

proc prinqual data=imc_dat plots=(MDPref) 

                 /* project onto Prin1 and Prin2 */

              ; /* use COV scaling */

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

   id cat;

   ods select MDPrefPlot;

run;