[求助]proc inbreed的输出结果如何调用

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

data inb;
input id sire dam;
cards;
1 . .
2 . .
3 . .
4 2 1
5 2 3
;
proc inbreed covar
   mattrix;
run;
/*------------------------------------*/
proc iml;
start main;
y={7,
   10,
   6};
x={1,
   1,
   1};
a={1  0  0 .5  0,
   0  1  0 .5 .5,
   0  0  1  0 .5,
  .5 .5  0  1 .25,
   0 .5 .5 .25 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;

麻烦哪位前辈帮我看看，前半个部分是用Proc inbreed得到近交系数矩阵，后半部分的a就是上半部分输出的矩阵，怎么直接调用proc inbreed的输出结果，不用自己手动输入阿？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：inbreed Breed 输出结果 Reed bree 结果 proc 输出 inbreed

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好像应该能够更简单些，我也在看。不知道能不能行。