/* 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;
