ods html;
DATA peso; 
INPUT x y @@;
x1=x;
x2=x**2;
x3=x2*x;
x4=x3*x;
x5=x4*x;
x6=x5*x;
FA=x; 
CARDS;
0  134.8  0   139.7  0 147.6	0   132.3
50 161.7 50   157.7	50 150.3 	50  144.7
100 160.7 100 172.7	100 163.4	100 161.3
150 169.8 150 168.2	150 160.7	150 161.0
200 165.7 200 160.0	200 158.2	200 151.0
250 171.8 250 157.3	250 150.4	250 160.4
300 154.5 300 160.4	300 148.8	300 154.0
; 
PROC GLM DATA=peso; 
TITLE'Análise de variância'; 
CLASS x; 
MODEL y = x / SS3; 
/*VERIFICANDO A FALTA DE AJUSTE*/
PROC GLM DATA=peso PLOTS=NONE;
CLASS FA;
MODEL y= x FA/SS1;
title 'TESTE DA FALTA DE AJUSTE DO MODELO LINEAR';
RUN;
PROC GLM DATA=peso;
CLASS FA;
MODEL y= x1 x2 FA/SS1;
title 'TESTE DA FALTA DE AJUSTE DO MODELO QUADRÁTICO';
RUN;
PROC GLM DATA=peso;
CLASS FA;
MODEL y= x1 x2 x3 FA/SS1;
title 'TESTE DA FALTA DE AJUSTE DO MODELO CÚBICO';
RUN;


/*TODOS OS EFEITOS*/
PROC GLM DATA=peso;
CLASS FA;
MODEL y= x1 x2 x3 x4 x5 x6/SS1;
title 'SOMAS DE QUADRADOS';
RUN;

/* MODELO SELECIONADO: QUADRÁTICO*/

/* Com a opção "/SOLUTION" é possivel obter as 
estimativas dos parâmetros do modelo de 
regressão */
PROC GLM PLOTS=NONE;
MODEL y = x1 x2/ SS1 SOLUTION;
RUN;
proc gplot data=peso;
plot y*x;
symbol1 v=dot i=rq c=red;
run;
ods html close;