options nodate nonumber;
data hierarq1;
input machos femeas descend PESO;
cards;
1      1      1      1550
1      1      2      1500
1      1      3      1550
1      1      4      1600
1      2      1      1700
1      2      2      1600
1      2      3      1300
1      2      4      1500
1      3      1      1600
1      3      2      1500
1      3      3      1550
1      3      4      1650
1      4      1      1350
1      4      2      1700
1      4      3      1500
1      4      4      2050
2      1      1      1550
2      1      2      1650
2      1      3      1500
2      1      4      1550
2      2      1      1400
2      2      2      1400
2      2      3      1350
2      2      4      1300
2      3      1      1550
2      3      2      1450
2      3      3      1650
2      3      4      1650
2      4      1      1300
2      4      2      1350
2      4      3      1500
2      4      4      1500
3      1      1      1650
3      1      2      1600
3      1      3      1600
3      1      4      1550
3      2      1      1300
3      2      2      1600
3      2      3      1250
3      2      4      1500
3      3      1      1250
3      3      2      1300
3      3      3      1350
3      3      4      1400
3      4      1      1700
3      4      2      1500
3      4      3      1700
3      4      4      1350
4      1      1      1650
4      1      2      1350
4      1      3      1250
4      1      4      1150
4      2      1      1350
4      2      2      1150
4      2      3      1250
4      2      4      1350
4      3      1      1300
4      3      2      1200
4      3      3      1300
4      3      4      1250
4      4      1      1250
4      4      2      1450
4      4      3      1450
4      4      4      1400
;
run;

proc print;
title 'Dados para verificação';
run;

proc glm;
title 'Análise de variância';
 class machos femeas descend;
 model PESO = machos femeas(machos);
 test h=machos e=femeas(machos);
run;