data inb;
input id sire dam;
cards;
1 . .
2 . .
3 . .
4 2 1
5 2 3
;
proc inbreed covar mattrix outcov=a;/*先把矩阵输入output到一个数据集a*/
run;

data a;/**只保留a中的矩阵部分前面的说明变量剔除掉，这个我没想到其他方法，做的有点笨，希望其他高手指点*/

set a;

keep col1-col5;

run;

proc iml;
start main;

use a；/*用数据集a*/

read all into a;/*读取为矩阵a*/
y={7,
   10,
   6};
x={1,
   1,
   1};

z={1 0 0 0 0,
   0 0 1 0 0,
   0 0 0 0 1};

lamda=1;

b=x`*x;
c=x`*z;
d=z`*x;
e=z`*z;
f=inv(a);
g=(e+(f#lamda));
h=(b||c)//(d||g);
l=x`*y;
m=z`*y;
o=(l)//(m);
p=inv(h);
q=p*o;

print b c d e f g h l m o p q;
finish main;
run;
quit;

多谢2楼，

好像代码不能执行哎！

proc inbreed里面有那个ODS有用么？看说明书好像行的。

我是菜鸟，呵呵

[此贴子已经被作者于2008-8-28 16:11:21编辑过]

proc inbreed种没有你的说的ODS啊？

你复制粘贴的问题吧？

找到错误了，有个分号是中文输入的。

多谢了！

在Proc inbreed帮助文件里面有关于PROC INBREED: ODS Table Names 的介绍，这个ODS好像应该能够更简单些，我也在看。不知道能不能行。