data ms;
input espe espac blo resp;
datalines;
1 1 1 10.69
1 1 2 10.05
1 1 3 10.31
1 2 1 12.95
1 2 2 10.99
1 2 3 12.37
2 1 1 9.68
2 1 2 10.58
2 1 3 10.44
2 2 1 11.78
2 2 2 12.28
2 2 3 12.82
3 1 1 9.59
3 1 2 9.78
3 1 3 10.17
3 2 1 11.46
3 2 2 10.89
3 2 3 10.97
4 1 1 9.4
4 1 2 7.5
4 1 3 8.66
4 2 1 10.02
4 2 2 10.45
4 2 3 10.71
;
proc print;
title 'Dados para verificação';
run;
proc glm data=ms PLOTS(ONLY)=(DIAGNOSTICS);
title 'Análise de variância e resíduos';
class espe espac blo;
model resp = espe espac blo espe*espac/ss3;
output out=residuos PREDICTED=pred RESIDUAL=res_ord STUDENT=res_stud;
lsmeans espe*espac/adjust=tukey pdiff slice=espe slice=espac;
run;
proc univariate data=residuos normal;
title 'Verificação da normalidade (teste Shapiro-Wilk)';
var res_stud;
run;
proc glm data=ms PLOTS=NONE;
title 'Análise de variância sem discriminar os fatores e Teste de Tukey para médias de tratamentos';
class trat;
model resp = trat /ss3;
means trat/Tukey;
run;
proc glm data=ms PLOTS=NONE;
title 'Análise de variância discriminando os fatores';
class espe espac blo;
model resp = espe espac blo espe*espac/ss3;
run;
/*daqui em diante o programa apresenta erro em sauser.letras*/
proc glm data=ms PLOTS(ONLY)=INTPLOT(LIMITS);
title 'Desdobramento Calagem dentro de cada Irrigação';
class Irrig Calag;
model prod = Irrig Calag Irrig*Calag/ss3;
lsmeans Irrig*Calag/slice=Irrig adjust=tukey PDIFF=all;
store sasuser.letras;
run;
proc PLM restore=sasuser.letras PLOTS=NONE;
lsmeans Irrig*Calag / lines adjust=tukey;
slice Irrig*Calag / sliceby=Irrig lines adjust=tukey;
run;
proc glm data=ms PLOTS(ONLY)=INTPLOT(LIMITS);
title 'Desdobramento Irrigação dentro de cada nível de Calagem';
class Calag Irrig;
model prod = Calag Irrig Calag*Irrig/ss3;
lsmeans Calag*Irrig/slice=Calag adjust=tukey PDIFF=all;
store sasuser.letras;
run;
proc PLM restore=sasuser.letras PLOTS=NONE;
lsmeans Calag*Irrig/lines adjust=tukey;
slice Calag*Irrig/ sliceby=Calag lines adjust=tukey;
run;