[求助]如何用SAS估算育种值

麻烦哪位高手帮我看一下哪里出现问题了，GDATA的赋值总是不对，没到求育种值程序就终止了。谢谢了！liqiuing2000@yahoo.com.cn

data inb;
input id sire dam;
cards;
2 . .
3 . .
4 . .
5 1 3
6 1 3
7 2 4
8 2 4
9 2 4
;
proc inbreed covar mattrix outcov=g;
run;

data g;
set g;
keep col1-col8;
run;

data ex1;
input farm $ animal $ y;
cards;
1 2 225
1 3 200
1 4 255
2 5 250
2 6 198
2 7 245
2 8 260
2 9 235
;

proc mixed ;
class farm animal;
model y=farm/solution;
random animal/gdata=g solution;
run;