Heterocedasticidade Regular e Transformação Box Cox (1/4/2020)

/* Exemplo Heterocedasticidade Regular

 e Transformação Box Coxs*/;

 /* Pequeno 2000 outlier */;

Title "Hetero. Regular";

Data Lucrativ;

Input Porte $  Lucro_L;

cards;

Pequeno 2000

Pequeno 198

Pequeno 300

Pequeno 140

Pequeno 201

Pequeno 350

Pequeno 150

Pequeno 199

Medio 800

Medio 410

Medio 270

Medio 800

Medio 390

Medio 250

Medio 430

Medio 700

Grande 600

Grande 1900

Grande 550

Grande 900

Grande 1950

Grande 900

Grande 700

Grande 2020

;

run;

proc print;

run;

proc glm;

class Porte;

Model Lucro_L = Porte;

Means Porte / Tukey Lines;

Run;

TITLE “DIC - Pressuposições do modelo”;

TITLE “Análise de variância e obtenção dos resíduos”;

PROC GLM DATA=Lucrativ;

  Class Porte;

  Model Lucro_L = Porte/SS3;

  OUTPUT OUT=residuos PREDICTED=y_est RESIDUAL=res_ord STUDENT=res_Stud;

RUN;

proc print data=residuos;

 run;

TITLE “Verificação da Normalidade dos resíduos, Normal-plot e Teste de Shapiro-Wilk”;

PROC UNIVARIATE DATA=residuos PLOT NORMAL;

VAR res_Stud;

RUN;

PROC UNIVARIATE DATA=residuos PLOT; VAR res_Stud;

RUN;

TITLE 'Análise gráfica dos resíduos e

 Teste de homogeneidade de variâncias de

 Brown-Forsythe';

 PROC PLOT DATA=residuos;

 PLOT res_Stud*trat="*"/vpos=20; RUN;

 PROC GLM DATA=Lucrativ;

   CLASS Porte;

   MODEL Lucro_L = Porte/SS3;

   MEANS Porte/HOVTEST=bf;

RUN;

TITLE 'Análise gráfica dos resíduos e Transformação de Box-Cox';

PROC PLOT DATA=residuos; PLOT res_Stud*y_est="*"/vpos=20;

RUN; ODS GRAPHICS ON;

PROC TRANSREG DATA=milho;

MODEL BOXCOX(Lucro_L/CONVENIENT LAMBDA=-2 to 2 BY 0.05) = CLASS(trat);

RUN; ODS GRAPHICS OFF;

*/

Programa Rodando

TITLE “DIC - Pressuposições do modelo”*/ ;

PROC SORT DATA=Lucrativ; BY Porte;

RUN;

TITLE “Análise de variância e obtenção dos resíduos”;

PROC GLM DATA=Lucrativ;

CLASS Porte;

MODEL Model Lucro_L = Porte/SS3;

OUTPUT OUT=residuos PREDICTED=y_est RESIDUAL=res_ord STUDENT=res_Stud;

RUN;

proc print data=residuos;

 run;

TITLE “Verificação da Normalidade dos resíduos,

       Normal-plot e Teste de Shapiro-Wilk”;

PROC UNIVARIATE DATA=residuos PLOT NORMAL;

VAR res_Stud;

RUN;

PROC UNIVARIATE DATA=residuos PLOT; VAR res_Stud;

RUN;

TITLE 'Análise gráfica dos resíduos e Teste de homogeneidade de

      variâncias de Brown-Forsythe';

PROC PLOT DATA=residuos; PLOT res_Stud*trat="*"/vpos=20;

RUN;

PROC GLM DATA=Lucrativ;

  class Porte;

  Model Lucro_L = Porte/SS3;

  MEANS Porte/HOVTEST=bf;

 RUN;

TITLE 'Análise gráfica dos resíduos e Transformação de Box-Cox';

PROC PLOT DATA=residuos;

PLOT res_Stud*y_est="*"/vpos=20;

 RUN;

 ODS GRAPHICS ON;

 PROC TRANSREG DATA=Lucrativ;

 MODEL BOXCOX(Lucro_L/CONVENIENT LAMBDA=-2 to 2 BY 0.05) = CLASS(Porte);

 RUN;

 ODS GRAPHICS OFF;